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Abstract 

Mismatched  epitaxial  layers  of  InaGai-^ajAs  and  InyAli_yAs  were  grown  on  (001) 
InP  by  molecular  beam  epitaxy.  The  layers  were  characterized  by  a  technique  we 
developed  known  as  variable  azimuthal- angle  elhpsometry.  It  reveals  large  optical 
anisotropy  for  many  strained  layers.  We  attribute  the  anisotropy  to  strain-induced 
surface  roughening  during  growth.  Samples  were  also  characterized  by  high-resolution 
x-ray  diffraction  (HRXRD)  to  assess  layer  quality  as  well  as  composition  and  strain. 
HRXRD  measurements  reveal  orthorhombic  distortion  of  partially  relaxed  layers  of 
InGaAs  and  InAlAs  in  tension  or  compression,  with  preferential  strain  relief  in  the 
[110]  direction.  We  show  that  HRXRD  epilayer  peak  width  and  interference  fringes 
are  sensitive,  non- destructive  criteria  to  judge  the  structural  quality  of  strained  het¬ 
erostructures.  For  layers  ranging  from  300  to  10,000  A,  with  lattice  mismatch  of 
±  1%  or  less,  the  crystalline  quality  consistently  remains  high  to  thicknesses  up  to 
3-9  times  the  Matthews-Blakeslee  critical  layer  thickness.  We  investigated  the  ther¬ 
mal  stability  of  these  layers,  using  HRXRD  to  measure  structural  changes  caused  by 
high-temperature  anneals.  We  also  compared  the  electron  mobility  of  modulation- 
doped  heterostructures  before  and  after  annealing.  Both  techniques  demonstrate  that 
our  high-quality  strained  layers  are  stable  to  temperatures  of  at  least  800-850® C.  We 
explain  this  result  by  the  limited  sources  available  for  the  nucleation  of  misfit  dis¬ 
locations.  The  findings  are  applied  to  the  design  and  growth  of  high-performance 
pseudomorphic  InAlAs/InGaAs/InP  heterostructure  field-effect  transistors  with  lay¬ 
ers  exceeding  the  Matthews-Blakeslee  limit. 
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Chapter  1 


Introduction 


1.1  Why  Ina.Gai_j.As/InP  and  InyAli_j/As/InP? 

Compound  semiconductors  are  of  interest  for  both  optical  and  electronic  devices. 
Solid-state  lasers  can  be  fabricated  from  III-V  semiconductors  such  as  gallium  arsenide 
(GaAs)  and  indium  phosphide  (InP)  because  these  materials,  unlike  silicon,  have 
fundamental  energy  gaps  which  are  direct.  These  materials  also  offer  the  possibility 
of  very  high-speed  electronic  devices.  If  only  binary  semiconductors  are  considered, 
however,  the  device  possibilities  are  somewhat  limited.  The  use  of  ternary  alloys  such 
as  Ala,Gai_j,As  and  In^Gai-seAs  allows  much  greater  freedom  to  choose  materials  with 
the  desired  characteristics  such  as  energy  gap  and  electron  mobihty. 

In  fig.  1.1,  we  plot  the  energy  gap  as  a  function  of  lattice  constant  for  selected  III- 
V  compound  semiconductors.  GaAs  and  AlAs  have  almost  identical  lattice  constants. 
Hence,  any  composition  of  AlxGai_xAs  will  be  nearly  lattice  matched  to  GaAs.  For 
this  reason,  much  research  has  focussed  on  the  AlxGai_xAs/GaAs  materials  system. 

In  contrast  to  GaAs,  nature  does  not  provide  a  ternary  alloy  system  which  is  al¬ 
ways  lattice  matched  to  InP.  As  shown  in  fig.  1.1,  InyAli_j,As  is  only  lattice  matched 
to  InP  when  y=0.52.  Similarly,  In^Gai-xAs  is  lattice  matched  if  and  only  if  x=0.53. 
If  the  quaternary  Ini-^-yGa^AlyAs  (with  x  +  y  K  0.47)  is  used,  alloys  with  energy 
gaps  covering  the  range  of  0.74  to  1.44  eV  can  be  lattice  matched  to  InP.  This  range 
includes  both  the  minimum  dispersion  wavelength  (1.3  fim)  and  the  minimum  at- 
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Figure  1-1:  Energy  gap  versus  lattice  constant  at  300K  for  selected  III-V  compound 
semiconductors  (courtesy  of  J.C.  Vlcek). 

tenuation  wavelength  (1.55  fim)  of  siUca  optical  fibers.  Primarily  for  this  reason, 
InP-based  compounds  are  the  materials  of  choice  for  lasers  and  detectors  in  optical 
telecommunications . 

Electronic  devices  based  upon  InP  are  also  of  interest  because  of  the  potential  for 
integration  with  optical  components  to  build  optoelectronic  integrated  circuits.  In 
addition,  the  InP  materials  system  has  several  intrinsic  advantages  for  electronic  de¬ 
vices.  For  example,  Ino.63Gao.47 As  has  a  room-temperature  mobility  of  11,000  cmVV- 
sec  and  a  peak  electron  velocity  of  3.0  x  10^  cm/sec}  These  values  are  considerably 
higher  than  those  for  GaAs  or  Si,  and  make  Ino.B3Gao.47As  a  promising  material  for 
high-frequency  microwave  transistors  and  high-speed  digital  circuits.  The  conduc¬ 
tion  band  offset  between  Ino.52Alo.48As  and  Ino.53Gao.47As  is  0.5  eV,  compared  to 
about  0.3  eV  for  the  AlGaAs/GaAs  system.^  Hence,  Ino.52Alo.48As  can  be  used  as  a 
pseudo-insulator  to  confine  electrons  in  the  Ino.53Gao.47 As  channel  of  heterostructure 
field-effect  transistors  (HFETs)  lattice-matched  to  InP. 

If  we  eliminate  the  requirement  of  lattice  matching,  the  range  of  available  material 
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properties  is  expanded.  For  example,  Inj,Gai_seAs  with  z  >  0.53  has  an  even  higher 
mobility  and  peak  velocity  than  Ino.B3Gao.47As,  and  can  be  used  as  the  channel  in 
HFETs.  The  lattice  mismatch  will  result  in  strained  epitaxial  layers.  In  the  case 
of  certain  optical  devices,  the  strain  itself  is  beneficial.  Strain  splits  the  light-  and 
heavy-hole  valence  band  degeneracy,  resulting  in  lower  threshold  current  for  solid- 
state  lasers.  Clearly,  the  use  of  mismatched  layers  in  the  InGaAs/InAlAs/InP  mate¬ 
rials  system  is  desirable  for  many  device  applications. 

1.2  Critical  Layer  Thickness:  Theory 

There  is  a  limit  to  the  amount  of  strain  which  can  be  accomodated  by  epitaxial  layers. 
If  a  layer  is  sufficiently  thick  and  mismatched,  the  formation  of  misfit  dislocations  at 
the  substrate/epilayer  interface  becomes  energetically  favorable.  These  dislocations 
reheve  strain,  but  also  degrade  the  structural,  electrical,  and  optical  quality  of  epi¬ 
taxial  layers,  often  making  them  unsuitable  for  device  applications. 

In  fig.  1.2,  we  show  schematics  of  epitaxial  growth  of  strained  layers.  We  assume 
the  epitaxial  layer  has  a  larger  lattice  constant  than  the  substrate.  If  the  lattice  mis¬ 
match  between  the  epilayer  and  substrate  is  small  and  the  layer  is  thin,  the  mismatch 
will  be  accomodated  entirely  by  strain  in  the  layer  and  atomic  bonds  will  be  continu¬ 
ous  across  the  interface  (A).  In  this  case,  the  symmetry  of  the  epilayer  lattice  distorts 
from  cubic  to  tetragonal,  and  the  layer  is  said  to  be  coherent  or  pseudomorphic.  As 
the  layer  thickness  is  increased,  misfit  dislocations  will  form,  reheving  a  portion  of  the 
strain  and  resulting  in  a  partially-relaxed  layer  (B).  Lattice  constants  perpendicular 
to  the  plane  of  growth,  a±,  and  parallel  to  the  plane  of  growth,  a||,  are  shown.  In 
case  A,  a||  is  equal  in  the  substrate  and  epilayer;  in  case  B,  the  values  are  different. 
In  both  A  and  B,  o_l  is  not  equal  in  the  substrate  and  epilayer.  Growth  modes  C  and 
D  will  be  discussed  later  in  this  chapter. 

For  a  fixed  lattice  mismatch,  the  thickness  at  which  misfit  dislocations  begin  to 
form  is  known  as  the  critical  layer  thickness.  Several  theories  have  been  proposed  to 
predict  the  critical  layer  thickness.  The  most  widely  accepted  theory  was  published 
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Figure  1-2:  Proposed  growth  modes  for  stained  epitaxial  layers.  A)  2-D  growth  of  a 
coherent  (dislocation-free),  tetragonally  distorted  layer.  B)  2-D  growth  of  a  partially 
relaxed  layer  in  which  misfit  dislocations  at  the  layer/substrate  interface  have  reUeved 
strain.  C)  2-D  growth  with  composition  modulation."  D)  3-D  growth  of  a  coherent 
island;  strain  is  relieved  by  the  deformation  of  the  substrate. 
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by  Matthews  and  Blakeslee  (M-B)  in  1974.*  M-B  considered  the  formation  of  misfit 
dislocations  at  the  layer/substrate  interface  by  the  elongation  of  threading  disloca¬ 
tions  from  the  substrate.  They  equated  the  force  exerted  on  the  dislocation  line  by 
misfit  stress  with  the  tension  in  the  dislocation  Une.  The  result  was  a  transcendental 
equation  relating  the  critical  layer  thickness  to  the  mismatch  and  elastic  constants. 
In  appendix  C,  we  give  the  equation  and  an  extension  of  the  theory  to  account  for 
anisotropic  elastic  properties.  To  first  order,  tc,MB>  the  M-B  critical  layer  thickness, 
is  inversely  proportional  to  the  lattice  mismatch,  /,  which  is  given  by: 

(1.1) 

o-s 

where  as  and  Ug  are  the  lattice  constants  of  the  substrate  and  epilayer,  respectively. 

To  get  some  idea  of  the  limitations  that  the  M-B  theory  imposes  on  strained-layer 
design  and  band-gap  engineering,  we  give  the  values  of  both  the  energy  gap®  and 
tc,MB  for  InAs-rich  InxGai_aAs  on  InP  in  table  1.1.  Complete  plots  of  tc,MB  for  both 
InxGai_xAs  and  In^Ali-yAs  are  given  in  appendix  C.  From  the  table,  we  see  that 
if  a  band-gap  of  0.07  eV  less  than  the  lattice  matched  value  is  required,  the  layer 
thickness  is  limited  to  300  A.  For  a  lattice  mismatch  of  1%,  the  critical  thickness  is 
only  about  130  A.  We  also  consider  a  layer  of  Ino.64oGao.46oAs  which  is  about  as  close 
to  lattice  matching  as  can  be  routinely  achieved  in  epitaxial  growth.  For  this  layer, 
ic,MB  is  only  4400  A;  thicker  layers  are  often  required  in  optical  devices.  It  is  clear 
that  the  M-B  limit  places  severe  constraints  on  the  design  of  strained-layer  and  even 
nominally  lattice-matched  heterostructures. 

An  epilayer  with  a  small  number  of  misfit  dislocations  may  stiU  be  suitable  for  cer¬ 
tain  device  applications.  In  fact,  InP  substrates  always  contain  threading  dislocations 
so  a  truly  defect-free  heterostructure  is  not  possible.  Elman  suggested  a  more  prac¬ 
tical  approach,  proposing  an  “effective  critical  layer  thickness,  defined  as  the  thickest 
layer  still  useful  for  high  quality  optical  devices.”^  This  effective  thickness  will  depend 
upon  the  specific  optical  (or  electronic)  device  in  question. 
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Table  1.1:  Lattice  mismatch,  energy  gap,  and  anisotropic  Matthews-Blakelsee  critical 
layer  thickness  for  InxGai-aAs/InP  at  300K. 


X 

f 

E,  (eV) 

tc,MB  (A) 

0.532 

0.00000 

0.736 

OO 

0.54 

0.00054 

0.728 

4400 

0.56 

0.00192 

0.707 

960 

0.60 

0.00469 

0.668 

300 

0.65 

0.00814 

0.620 

170 

0.70 

0.0116 

0.574 

110 

0.80 

0.0185 

0.490 

60 

1.00 

0.0323 

0.350 

30 

1.3  Critical  Layer  Thickness:  Experiments 

Over  the  past  decade,  a  large  number  of  experimental  investigations  of  critical  layer 
thickness  have  been  reported.®’®  Experimental  techniques  have  ranged  from  imaging 
dislocations  with  electron  microscopy  to  fabricating  heterostructure  devices  and  in¬ 
ferring  the  presence  of  dislocations  from  device  characteristics.  The  results  have  often 
been  contradictory,  with  critical  layer  thicknesses  ranging  from  about  tc,MB  to  more 
than  an  order  of  magnitude  greater  than  We  will  review  experimental  results 

in  this  and  the  following  sections.* 

Orders  and  Ushers^  used  high-resolution  x-ray  diffraction  to  measure  the  strain 
rehef  in  InGaAs  layers  on  (001)  GaAs  substrates.  They  found  that  layers  remained 
coherent  to  thicknesses  about  an  order  of  magnitude  larger  than  Gal  et  al} 

examined  the  photoluminescence  (PL)  peak  position  and  width  in  the  InGaAs/GaAs 
system,  and  obtained  results  similar  to  those  of  Orders  and  Usher.  In  the  same 
materials  system,  however,  Fritz  et  al^  apphed  both  PL  and  low-temperature  Hall 
measurements  and  obtained  fc’s  approximately  equal  to  U,mb-  Wang^®  fabricated 
AlGaAs/InGaAs/GaAs  HFET’s  with  the  InGaAs  channels  a  factor  of  two  and  four 

‘Most  of  the  experimental  investigations  of  critical  layer  thickness  have  been  in  the  SiGe/Si  and 
InGaAs/GaAs  systems,  with  only  a  few  reports  for  InGaAs/InP  and  InAlAs/InP  heterostructures. 
We  expect  that  InGaAs/InP  and  InAlAs/InP  may  behave  similarly  to  InGaAs/GaAs,  and  hence 
include  the  latter  in  our  literature  review. 
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Figure  1-3:  Critical  layer  thickness,  as  determined  by  Hall  mobility,  versus  InAs 
mole-fraction  in  Ino.52Alo.48As/Ina,Gai_a,As/InP  modulation-doped  structures  (from 
Tacano  et  al}^). 

thicker  than  tc,MB,  and  observed  no  signs  of  degradation  in  the  device  characteristics. 

In  the  InGaAs/InP  system,  conflicting  results  have  also  been  reported.  Tacano 
et  al}^  measured  the  lOK,  77K,  and  300K  Hall  mobility  of  modulation-doped  Ino.62* 
Alo.48As/Ina,Gai_j,As/InP  heterostructures.  As  shown  in  fig.  1.3,  they  obtained  criti¬ 
cal  thicknesses  which  exceeded  tc,MB  by  a  factor  of  three  to  five.  Temkin^^  also  studied 
the  InxGai_a,As/InP  system,  but  looked  for  signs  of  misfit  dislocations  by  monitoring 
the  reverse-bias  leakage  current  in  photodiodes.  His  results  were  approximately  in 
agreement  with  the  Matthews-Blakeslee  theory,  as  shown  in  fig.  1.4. 

It  is  now  clear  that  one  reason  for  the  apparent  discrepancies  in  the  literature  is  the 
different  experimental  resolutions  of  the  various  techniques.^^’^^  In  the  InGaAs/InP 
examples  listed  above,  it  seems  reasonable  that  a  significant  number  of  misfit  dislo¬ 
cations  are  required  before  electron  mobility  is  affected,  but  each  dislocation  could 
make  a  substantial  contribution  to  the  leakage  current.  It  is  apparent  that  any  serious 
study  of  critical  layer  thickness  requires  a  thorough  analysis  of  the  characterication 
techniques. 
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Figure  1-4:  Strained-layer  thickness  versus  InAs  mole  fraction  for  ^*Gai_xAs/InP 
p-i-n  photodiodes.  Solid  symbols  indicate  samples  which  contain  misfit  dislocations, 
based  upon  the  reverse-leakage  current.  Samples  which  exhibit  no  signs  of  disloca¬ 
tions  are  plotted  as  open  symbols.  The  line  is  the  Matthews-Blakeslee  critical  layer 
thickness  (from  Temkin  et 

1.4  Metastability 

The  Matthews-Blakeslee  theory  of  the  critical  layer  thickness  is  based  upon  equilib¬ 
rium  considerations.  A  strained  layer  exceeds  tc,MB  if  the  lowest  energy  configuration 
includes  misfit  dislocations.  In  techniques  such  as  molecular  beam  epitaxy  (MBE) 
and  metal-organic  chemical  vapor  deposition  (MOCVD),  growth  occurs  under  condi¬ 
tions  which  are  far  from  equilibrium. “  Hence,  it  is  possible  that  kinetic  bamers  to 
the  nucleation  and  growth  of  misfit  dislocations  could  result  in  metastable  layers} 

The  degree  of  relaxation  in  such  layers  can  be  a  function  of  the  growth  temperature, 
as  demonstrated  by  Whaley  and  Cohen  in  the  InGaAs/GaAs  system.^*  They  used 
reflection  high-energy  electron-diffraction  (RHEED)  to  monitor  the  value  of  a||  in 
the  epilayer  during  growth.  The  InGaAs  relaxed  much  sooner  when  a  high  growth 
temperature  (510"C)  was  used  compared  to  a  low  temperature  (470“C).  Elman  et 
al}  applied  PL  and  transmission  electron  microscopy  (TEM)  to  layers  of  InGaAs  on 
GaAs  and  observed  an  increase  in  the  critical  thickness  by  a  factor  of  seven  when  low 
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growth  temperatures  were  used. 

A  potential  problem  with  the  use  of  metastable  layers  in  devices  is  the  possibil¬ 
ity  of  layer  relaxation  during  high-temperature  processing  steps  or  device  operation. 
Peercyet  a/.^°  tested  this  hypothesis  by  annealing  InGaAs/GaAs  heterostructures  in 
which  the  epilayer  exceeded  tc.MB-  They  observed  a  dramatic  decrease  in  photolumi¬ 
nescence  intensity  after  the  anneal,  suggesting  lattice  relaxation.  A  similar  photolumi¬ 
nescence  experiment  by  Bertolet  et  yielded  the  opposite  results:  InGaAs/GaAs 
heterostructures  were  thermally  stable,  despite  exceeding 

Metastability  was  also  investigated  by  Fitzgerald  tt  who  compared  the  growth 
of  strained  layers  of  InGaAs  on  GaAs  mesas  and  planar  substrates.  They  found  that 
the  dislocation  density  was  lower  for  the  growth  on  mesas.  The  results  are  explained 
by  the  hmited  number  of  threading  dislocations  in  the  mesas  and  the  lack  of  suffi¬ 
cient  energy  for  the  homogeneous  nucleation  of  additional  dislocations.  Fitzgerald 
concludes  that  diiferent  apparent  critical  layer  thicknesses  may  be  observed,  depend¬ 
ing  on  the  active  dislocation  sources. 


1.5  Growth  Modes 

In  fig.  1.2,  cases  A  and  B,  we  illustrated  a  simple  model  in  which  growth  proceeds  in 
a  layer-by-layer  fashion,  with  the  introduction  of  misfit  dislocations  to  reheve  strain. 
The  Matthews-Blakeslee  theory  is  based  upon  this  type  of  model.  In  recent  years,  it 
has  become  clear  that  growth  of  strained  layers  is  often  more  complex.  For  example, 
both  phase  separation  and  three-dimensional  growth  have  been  observed  in  strained 
InGaAs  and  InAlAs  epilayers. 

Epitaxial  growth  has  been  classified  into  three  distinct  modes. We  will  briefly 
review  each  mode  and  then  discuss  experimental  results  for  InGaAs  and  InAlAs.  In 
the  Frank-van  der  Merwe  mode,  epitaxial  growth  takes  place  in  a  two-dimensional 
manner,  with  the  completion  of  one  monolayer  before  the  next  begins. f  In  the  Vb/mer- 

^  Recent  work  on  the  AlGaAs/GaAs  system  suggests  that  the  growth  front  may  actually  extend 
over  a  few  monolayers  in  so-called  layer-by-layer  growth.^^ 
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Weber  mode,  epitaxial  growth  occurs  in  a  three-dimensional  manner  by  the  formation 
and  coalescence  of  islands.  The  third  mode,  Stransh-Krastanov,  is  a  hybrid  with 
initial  layer-by-layer  growth  followed  by  the  formation  of  islands.  The  growth  mode 
of  any  system  will  be  determined  by  the  surface  energies  of  the  layer  and  substrate 
and  the  interfacial  energy.  Strain  can  be  an  important  factor  because  it  results  m  a 
higher  surface  energy  for  the  layer.**’ 

Tanaka  and  Ohkouchi**  investigated  the  growth  of  GaAs  on  InP  by  scannmg 
tunneling  microscopy.  They  observed  a  Stranski-Krastanov  mode,  with  2  D  growth 
for  the  first  two  monolayers,  followed  by  island  formation.  Hopkinson  et  o/.*®  grew 
quantum  wells  of  InAs  on  InP.  Based  upon  cross-sectional  TEM,  the  growth  remamed 
planar  to  at  least  30  A.  Similarly,  de  Miguel**  observed  2-D  growth  of  up  to  30  A 
of  InAs  on  Ino.82Alo.48As/InP.  Since  the  Stranski-Krastanov  mode  occurs  in  these 
extreme  cases  of  mismatch,  we  believe  that  growth  of  InGaAs  or  InAlAs  on  InP  will 
not  occur  in  the  Volmer-Weber  mode  under  normal  conditions.  Throughout  this 
thesis,  we  will  often  use  the  term  3-D  growth  when  actually  referring  to  Stranksi- 


Krastanov  growth. 

We  define  a  threshold  thickness,  ttfc,  as  the  thickness  at  which  the  growth  mode 
changes  from  2-D  to  3-D.  In  general,  we  are  interested  in  both  Uh  and  U,  the  critical 
thickness  for  formation  of  misfit  dislocations.  We  wish  to  know  which  thickness  is 
smaller  and  the  degree  of  layer  degradation  when  each  thickness  is  exceeded. 

One  view  of  strained  layer  growth  is:  a  layer  initiaUy  grows  in  a  2-D  mode;  af¬ 
ter  the  critical  layer  thickness  is  exceeded,  misfit  dislocations  form  to  reheve  strain; 
the  dislocations  serve  as  sites  for  the  nucleation  of  clusters  and  the  growth  mode 
becomes  3-D.*^’**  An  alternative  view,  supported  by  experimental  work  in  both  the 
SiGe/Si**  and  InGaAs/GaAs**’*®  systems,  suggests  that  strained  layers  may  switch  to 
the  3-D  growth  mode  in  the  absence  of  misfit  dislocations.*  This  so-called  “coherent 
Stranski-Krastanov”  mode  is  explained  by  the  accomodation  of  mismatch  by  elastic 
deformation  around  islands,  as  illustrated  in  fig.  1.2D. 


1  We  note  that  the  authori“~®®  u^d  electron  microscopy  to  image  idands  in  which  no  diriocations 
were  visible.  One  cannot,  however,  prove  the  absence  of  dislocations  in  t^  way.  ^nce,  it  is  more 
precise  to  say  that  3-D  growth  occurred  in  the  absence  of  significant  numbers  of  misfit  dislocataons. 
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Snyder  et  examined  InGaAs  layers  on  GaAs  by  both  RHEED  and  TEM. 
Based  upon  the  experimental  results,  they  postiUate  that  the  strain  initially  induces 
a  roughening  of  the  layer.  Eventually,  as  strain  energy  builds  up  with  each  successive 
mismatched  monolayer,  islands  form  in  a  “kinetically  controlled  coarsening  process.” 
This  scenario  allows  for  a  transition  to  3-D  growth  without  the  formation  of  misfit 
dislocations. 

Based  upon  a  theoretical  analysis  of  strain  energies  in  2-D  and  3-D  growth,  Berger 
et  al.^^  demonstrate  that  the  lowest  energy  surface  of  an  epilayer  may  be  3-D  in  the 
presence  of  strain.  They  apply  their  results  to  the  InGaAs/ GaAs  system  and  estimate 
that  3-D  growth  should  occur  when  the  mismatch  exceeds  2%.  Several  experimental 
studies  have  confirmed  3-D  growth  for  mismatches  greater  than  about  2%.®^“^*  Only 
Chang  et  al.^^  characterized  relatively  thick  layers  with  mismatches  less  than  2%. 
They  saw  no  signs  of  3-D  growth  for  a  1.0  ^m  layer  with  a  1%  mismatch. 

Very  recent  work  by  Gendry  et  examined  the  growth  mode  for  compressively- 

strained  layers  of  InGaAs  on  InP.  Using  RHEED  oscillations,  they  observed  a  tran¬ 
sition  from  2-D  to  3-D  growth  for  mismatches  of  2%  and  3%.  The  value  of  Uh 
increased  with  decreasing  growth  temperature.  Apparently,  growth  at  higher  tem¬ 
peratures  increases  the  mobility  of  surface  atoms,  allowing  the  formation  of  islands. 
For  a  mismatch  of  1%,  the  layer  remained  in  a  2-D  growth  mode  up  to  2200  A;  thicker 
layers  were  not  grown.  The  Gendry  results  support  the  idea  that  coherent  islands 
may  form  in  some  situations. 

Lievin  and  Fonstad^®  measured  RHEED  oscillations  during  the  growth  of  compres- 
sively-strained  InAlAs  on  InP  with  mismatches  from  1  to  3%.  They  observed  the 
apparent  onset  of  3-D  growth  at  thicknesses  ranging  from  10  A  (3%)  to  100  A  (1%). 
3-D  growth  began  earlier  for  layers  in  tension  with  comparable  mismatch.  They 
attribute  this  to  the  lower  In  composition  and  resulting  decrease  in  average  cation 
mobility  at  a  growing  surface. 

It  is  clear  that  strained  InGaAs  and  InAlAs  layers  frequently  grow  in  a  3-D  mode. 
If  the  islands  are  coherent,  the  material  quahty  may  remain  relatively  high.  Even  in 
this  case,  however,  the  variation  in  thickness  of  a  quantum  well  could  broaden  the 
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luminescence  peak  width.  Coherent  islands  can  only  accomodate  a  relatively  small 
amount  of  strain.  Eventually,  islands  will  coalesce,  forming  undesirable  defects  such 
as  edge  dislocations,  stacking  faults,  and  threading  dislocations. 

Another  possible  occurrence  during  the  growth  of  epitaxial  layers  is  a  modulation 
of  the  composition  in  a  direction  parallel  to  the  sample  surface.  This  phenomenon 
is  often  referred  to  as  phase  separation,  and  is  analogous  to  spinodal  decomposition 
in  bulk  alloys. Theoretical  work^®’'*®  has  shown  that  layers  with  composition  mod¬ 
ulation  may  have  a  lower  energy  than  homogeneous  layers.  We  show  a  schematic  of 
the  modulation  in  fig.  1.2C.  Note  the  change  in  lattice  parameter,  corresponding  to 
a  change  in  composition,  for  layers  near  the  surface. 

Composition  modulation  has  been  observed  in  ternary  and  quaternary  alloys  of 
III-V  compounds.®®"®®  In  the  case  of  InGaAs/InP,  the  period  of  oscillation  varied 
from  2000  to  4000  A,  and  decreased  with  increasing  thickness  for  a  fixed  mismatch.®^ 
The  composition  modulation  was  present  for  layers  with  thicknesses  of  1  —  2  x 
but  no  lattice  defects  were  observed.  For  thicker  InGaAs  layers,  both  composition 
modulation  and  stacking  faults  were  found.  These  results  suggest  that  composition 
modulation  can  accomodate  a  limited  amount  of  layer  strain,  but  eventually  other 
defects  will  be  introduced. 


1.6  Summary 

In  summary,  the  theory  of  Matthews  and  Blakeslee  predicts  an  equilibrium  critical 
layer  thickness  for  strained  epitaxial  layers.  Experimental  work  suggests  that  layers 
thinner  than  tc,MB  are  always  free  of  misfit  dislocations.  The  M-B  hmits,  however, 
place  severe  constraints  on  the  use  strained  layers  in  electronic  and  optical  devices. 
Layers  exceeding  tc,MB  are  clearly  desirable  in  many  applications. 

Several  important  questions  remain  for  the  In*Gai_*As/InP  and  InyAli_j,As/InP 
material  systems:  What  are  the  limits  on  layer  mismatch  and  thickness  for  high- 
quahty  heterostructures?  What  mechanisms  (misfit  dislocations,  three-dimensional 
growth,  composition  modulation,  etc.)  degrade  the  mismatched  layers?  Are  strained 
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layers  thermally  stable?  Why  or  why  not?  This  thesis  will  address  these  questions. 

1.7  Outline  of  Thesis 

As  mentioned  above,  any  serious  study  of  strained  layer  epitaxy  requires  a  careful 
analysis  of  the  characterication  techniques.  In  chapter  2,  we  analyze  the  use  of 
high-resolution  x-ray  diffraction  (HRXRD)  to  assess  the  composition,  relaxation,  and 
structural  quality  of  mismatched  epilayers.  We  demonstrate  that  the  crystalline  qual¬ 
ity  of  InGaAs  and  InAlAs  layers  consistently  remains  high  to  thicknesses  up  to  3-9 
times  the  Matthews-Blakeslee  critical  layer  thickness.  For  thicker  layers,  lattice  re¬ 
laxation  occurs  in  an  asymmetric  manner,  with  a  change  in  crystal  symmetry  from 
tetragonal  to  orthorhombic. 

Chapter  3  describes  a  new  characterization  technique  we  developed  known  as 
variable  azimuthal  angle  eUipsometry.  It  reveals  large  optical  anisotropies  in  mis¬ 
matched  InGaAs  and  InAlAs  layers.  We  attribute  the  anisotropy  to  strain-induced 
surface  roughening. 

In  chapter  4,  we  assess  the  thermal  stability  of  strained  InGaAs  and  InAlAs 
layers.  Using  HRXRD  and  HaU  mobility  measurements,  we  demonstrate  that  high- 
quality  layers  beyond  the  Matthews-Blakeslee  critical  thickness  are  relatively  stable 
to  temperatures  of  at  least  800-850“ C.  We  explain  this  result  by  the  limited  source 
of  threading  dislocations  and  the  lack  of  sufficient  thermal  energy  to  homogeneously 
nucleate  half-loop  dislocations  during  annealing. 

As  mentioned  earlier,  the  motivation  for  the  study  of  strained  layers  is  their  po¬ 
tential  application  in  devices.  We  applied  our  findings  to  the  HFET.  As  described 
in  chapter  5,  we  grew  HFETs  with  strained  Ina,Gai_j,As  channels  and  InyAlj.yAs 
pseudo-insulators.  The  InyAli_yAs  layers  were  about  twice  the  Matthews-Blakeslee 
critical  thickness.  HRXRD  showed  the  heterostructures  to  be  of  good  crystalline 
quality.  The  device  results  were  excellent,  suggesting  that  if  misfit  dislocations  were 
present,  they  did  not  affect  the  device  performance  in  any  appreciable  way. 

In  chapter  6,  we  discuss  our  results  and  the  relevance  of  theories  for  lattice 
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relaxation  and  strained  layer  growth.  We  concur  with  the  idea  that  a  source  of 
lattice  imperfections  is  sometimes  necessary  to  provide  the  nucleation  sites  for  misfit 
dislocations.  This  model  can  explain  the  thermal  stability  of  our  high-quality  layers 
as  well  as  the  instability  of  other  layers  reported  in  the  literature.  If  the  surface  of 
a  strained  layer  roughens  or  the  layer  grows  in  a  3-D  mode,  the  resulting  defects  can 
serve  as  a  source  of  misfit  dislocations  for  strain  relaxation  during  subsequent  growth 
or  anneaUng. 

The  key  conclusions  are  summarized  in  chapter  T .  We  also  discuss  promising 
directions  for  future  research  including  experiments  which  could  verify  our  model  for 
thermal  stability  and  lattice  relaxation. 

The  growth  of  over  200  InGaAs/InP  and  InAlAs/InP  heterostructures  by  solid- 
source  MBE  was  an  important  part  of  this  work.  We  generally  used  well-estabhshed 
growth  procedures  which  are  documented  in  the  literature.  We  describe  our  epilayer 
growth  in  appendix  A.  In  appendix  B,  we  compile  the  important  materials  prop¬ 
erties  of  InP,  InAs,  GaAs,  and  AlAs.  We  frequently  use  interpolations  from  these 
values  throughout  this  thesis.  Appendix  C  gives  the  equations  for  the  Matthews- 
Blakeslee  critical  layer  thickness  and  its  extension  to  the  case  of  elastically  anisotropic 
materials. 
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Chapter  2 


Characterization  by 
High- Resolution  X-Ray 
Diffraction 


High- resolution  x-ray  diffraction  (HRXRD)  is  routinely  used  to  measure  the  compo¬ 
sition  and  strain  of  epitaxial  layers  of  semiconductor  alloys  such  as  InxGai_*As.®^"®® 
We  apply  it  for  this  purpose  throughout  the  thesis  and  investigate  the  accuracy  of 
such  measurements.  We  focus  our  efforts  on  the  use  of  HRXRD  to  determine  the 
crystalline  quality  of  epitaxial  layers.  Compared  to  characterization  techniques  such 
as  transmission  electron  microscopy  and  the  fabrication  and  testing  of  heterostructure 
devices,  HRXRD  is  relatively  fast  and  simple.  Hence,  we  are  able  to  characterize  a 
large  number  of  layers  covering  a  wide  range  of  thickness  and  mismatch.  The  results 
allow  us  to  determine  the  range  of  thickness  and  mismatch  for  MBE  growth  of  high- 
quality  InGaAs  and  InAlAs  layers  on  InP.  In  later  chapters,  we  apply  HRXRD  to 
measure  the  thermal  stabihty  of  strained  layers  and  correlate  HRXRD  measurements 
of  crystalline  quality  with  electron  mobility  and  device  performance. 

In  section  2.1,  we  briefly  review  the  theory  and  experimental  techniques  required 
to  measure  layer  composition.  We  also  discuss  the  role  of  simulations  based  upon  the 
dynamical  theory  of  x-ray  diffraction.  Then,  we  present  our  contributions  in  three 
areas:  Section  2.2  describes  our  measurements  of  layer  quality  and  its  relationship 
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to  lattice  mismatch  and  thickness.  We  demonstrate  that  high-quality  InGaAs  and 
InAlAs  layers  can  be  grown  to  thicknesses  beyond  the  Matthews-Blakeslee  critical 
thickness.  In  section  2.3,  we  show  that  partially  relaxed  layers  of  both  InGaAs  and 
InAlAs  exhibit  orthorhombic  distortion  as  a  result  of  an  asymmetry  in  dislocation 
density.  This  fact  must  be  taken  into  account  when  measuring  the  composition  of 
such  layers.  In  section  2.4  we  determine  the  range  of  thicknesses  for  which  layer  com¬ 
position  can  be  accurately  determined  from  a  single  HRXRD  scan.  This  information 
is  important  for  the  design  of  calibration  layers  and  can  also  be  applied  to  multilayer 
structures  including  devices.  The  work  is  summarized  in  section  2.5. 


2.1  Background 

When  electromagnetic  radiation  with  a  wavelength  comparable  to  the  atomic  spac¬ 
ing  strikes  a  crystal,  atomic  scattering  wiU  result  in  constructive  interference  and 
a  diffracted  beam  in  certain  directions.  The  fundamental  equation  describing  the 
diffraction  is  Bragg’s  law: 


nX  =  2dsin6B  (21) 

where  6b  is  the  Bragg  diffraction  angle,  X  is  the  wavelength  of  the  radiation,  n  is  the 
order  of  the  diffraction,  and  d  is  the  interatomic  spacing,  d  is  related  to  the  lattice 
constant,  a,  by: 


Vh^  +  k^  +  P 


(2.2) 


where  (hkl)  are  the  indices  of  the  diffracting  plane. 

Constructive  interference  results  when  (2.1)  is  satisfied.  If  a  single  crystal  is 
scanned  through  6  to  produce  a  plot  of  intensity  versus  angle,  a  peak  will  be  observed 
at  6b-  The  spectral  width  of  the  source  combined  with  beam  divergence  will  result  in 
a  peak  width  of  several  minutes  of  arc.  This  resolution  is  not  good  enough  to  measure 
the  small  lattice  constant  differences  typically  found  in  epitaxial  semiconductor  layers. 
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1st  CRYSTAL  (REFERENCE) 


Figure  2-lt  Schematic  of  double-crystal  x-ray  diffraction  system.  The  case  of  sym¬ 
metric  reflections  from  both  crystals  is  shown. 

If  instead  an  x-ray  source  is  diffracted  off  a  reference  (first)  crystal  and  the  diffracted 
beam  is  then  used  to  produce  a  rocking  curve  off  a  second  crystal  (the  sample),  peak 
widths  as  small  as  a  few  arc-seconds  can  be  achieved.  This  measurement  technique 
is  known  as  double-crystal  x-ray  diffraction  and  is  illustrated  in  fig.  2.1.  For  some 
applications,  the  addition  of  a  third  and  even  fourth  crystal  can  be  advantageous.  We 
use  the  term  high-resolution  x-ray  diffraction  to  refer  to  any  system  employing  two 
or  more  crystals. 

When  using  HRXRD  to  measure  layer  composition,  the  only  parameter  required 
from  the  experiment  is  A9b,  the  difference  in  Bragg  angle  of  the  substrate  and  epi- 
luyer.  To  derive  the  relationship  between  and  composition,  we  combine  (2.1) 
and  (2.2): 


nXVh^  +  k^  +  P 
2sineB 

We  consider  the  symmetric  case  in  which  the  diffraction  planes  are  parallel  to  the 
substrate  surface,  and  the  incidence  angle  equals  the  exit  angle.  Differentiating  (2.3) 
with  respect  to  6  and  dividing  the  resulting  equation  by  (2.3)  yields: 
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A6b 

tanSs 


(2.4) 


(“) 

\  a  J  ±  zanuB 

where  we  have  used  a  small-angle  approximation.  We  can  also  express  (2.4)  as: 


^-KMb  (2.5) 

where  A"  is  a  constant  determined  by  the  lattice  constant  of  the  substrate  and  the 
diffraction  plane,  and  is  expressed  in  arc-sec.  We  include  the  ±  subscript  in  (2.4) 
and  (2.5)  because  a  symmetric  measurement  is  only  sensitive  to  the  lattice  mismatch 
perpendicular  to  the  substrate/epilayer  interface,  (Ao/a)x  (see  fig.  1.2).  Measure¬ 
ment  of  the  lattice  mismatch  parallel  to  the  interface,  (Aa/a)||,  requires  asymmetric 
measurements  and  will  be  discussed  in  section  2.3. 

In  order  to  relate  composition  to  lattice  mismatch,  it  is  convenient  to  define  a 
parameter  known  as  the  relaxed  lattice  mismatch.  It  is  the  mismatch  a  layer  would 
have  if  it  were  totally  relaxed  or  in  bulk  form,  and  is  conventionally  expressed  as:®® 


(2.6) 

\  a  J  r  1  +  V  \  a  J  X  1  +  V  \  a  /  \\ 

where  t/  is  Poisson’s  ratio  for  the  layer.  For  III-V  semiconductors,  u  is  typically  about 
1/3.  If  a  layer  is  fully  relaxed,  (Ao/o)j.  =  (Ao/a)||  =  (Aa/a),.  For  partially  relaxed 
layers,  (Aa/a)||  <  (Aa/a),.  <  (Aa/a)j..  If  a  layer  is  coherent,  (Aa/a)||  =  0,  and: 


The  final  step  is  to  relate  the  relaxed  lattice  mismatch  to  the  layer  composi¬ 
tion.  We  use  Vegard’s  law  which  assumes  that  the  lattice  constant  of  an  alloy  {e.g. 
InxGai_j,As)  is  a  linear  function  of  the  lattice  constants  of  the  constituent  binaries 
{e.g.  InAs  and  GaAs).  For  Ina,Gai_j,As,  we  have: 

X  =  14.48  0.5322  (2.8) 

using  the  lattice  parameters  in  appendix  B.  Similarly,  for  InyAli_yAs:®° 
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(2.9) 


y  =  14.82  +  0.5210 

We  collected  x-ray  data  with  Bede  model  300  and  model  150  two-crystal  systems 
using  Cu-Ka  radiation.  The  first  crystal  was  always  InP  oriented  for  the  (004)  reflec¬ 
tion.  Rocking  curves  were  measured  for  symmetric  (004)  as  well  as  asymmetric  (115) 
and  (224)  reflections  off  the  second  crystal  (sample  being  characterized).  In  a  typical 
scan,  data  points  were  collected  every  2  arc-sec  with  a  count  time  of  2  seconds  per 
point.  Such  scans  required  30-120  minutes,  depending  on  the  sample  mismatch  and 
resulting  angular  separation  between  substrate  and  layer  peaks.  When  measuring 
very  thin  layers  (100-200  A),  we  often  used  count  times  of  5-10  seconds,  requiring 
as  long  as  8  hours  per  scan.  The  spot  size  of  the  x-ray  beam  on  the  sample  was 
approximately  1  mm  by  2  mm. 

A  typical  HRXRD  scan  is  shown  in  fig.  2.2a.  The  sample  consisted  of  a  single  layer 
of  InyAli_yAs  on  InP;  the  symmetric  (004)  reflection  geometry  was  used.  (We  do  not 
measure  absolute  Bragg  angles,  only  differences  between  Bragg  angles.  Throughout 
this  work,  we  set  the  InP  peak  to  be  at  0  arc-seconds  for  convenience.)  In  this  case, 
both  the  substrate  and  the  layer  produce  distinct  peaks.  The  separation  between 
the  peaks,  A^,  is  -152  arc-sec.  Since  the  layer  has  a  smaller  Bragg  angle  than  the 
substrate,  the  lattice  constant  of  the  layer  must  be  larger  (eq.  2.1),  making  it  InAs- 
rich  compared  to  the  lattice-matched  composition  of  Ino.B2iAlo.479As.  From  appendix 
B,  the  lattice  constant  of  InP  is  5.8688  A.  Using  (2.1)  and  (2.2)  with  a  wavelength  of 
1.54056  A  (CuKoii),  we  obtain  31.668°  for  the  (004)  InP  Bragg  angle.  The  K  in  (2.5) 
is  7.8595  X  10"®,  and  the  perpendicular  mismatch  is  -1-0.00119.  If  we  now  assume  the 
layer  is  coherent  ((Aa/a)||  =  0),  (2.6)  gives  a  relaxed  mismatch  of  -1-0.000618.  Using 
(2.9),  the  layer  composition  is  Ino.53oAlo.47oAs.  We  will  discuss  the  accuracy  of  such 
measurements  in  section  2.4. 

In  fig.  2.2b,  we  replot  the  data  of  fig.  2.2a  on  a  logarithmic  scale.  We  ob¬ 
serve  not  only  the  substrate  and  layer  peaks,  but  also  a  series  of  low-intensity  peaks 
with  constant  spacing.  These  peaks  result  from  interference  effects  between  the  sub- 
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Figure  2-2:  HRXRD  (004)  rocking  curve  for  sample  1683,  1760  A  of  Ino.53oAlo.47oAs 
on  InP  plotted  on  linear  (a)  and  semi-log  (b)  scales.  Pendellosung  fringes  are  clearly 
visible  in  (b). 


strate/layer  and  layer/air  interfaces  and  are  known  as  Pendellosung  or  interference 
fringes.  For  a  single  layer,  the  spacing  between  the  fringes,  A/,  is  related  to  the  layer 
thickness,  t,  by:®®’®® 


A/ 


(2.10) 


tsin{20B) 

where  7/,  is  the  cosine  of  the  angle  between  the  diffracted  beam  and  the  surface  normal. 
In  the  case  of  symmetric  reflections,  7^  =  cos(90-5b).  For  the  (004)  reflection,  (001) 
InP  substrates,  and  CuKa  radiation: 


t  = 


18.67 

A/ 


(2.11) 


where  t  is  in  microns  and  A/  is  in  arc-sec.  For  sample  1683  in  fig.  2.2,  A/  —  106 
arc-sec  and  t  =  1760  A. 

We  can  now  check  our  assumption  of  coherency  for  the  epilayer  in  sample  1683. 
From  appendix  C,  the  modified  Matthews-Blakeslee  critical  thickness,  tc,MB,  for 
Ino.B3oAlo.470 As  is  3000  A.  Since  t  <  tc,MB)  the  layer  is  expected  to  be  coherent, 
and  our  calculated  composition  should  be  correct. 

The  identification  of  diffraction  peaks  is  usually  straightforward  in  the  case  of  a 
single  layer  on  a  substrate  such  as  sample  1683.  In  many  cases,  however,  we  wish  to 
characterize  multi-layer  heterostructures.  In  such  cases,  an  x-ray  scan  often  reveals  a 
plethora  of  diffraction  peaks.  The  clear  distinction  between  primary  diffraction  peaks 
and  Pendellosung  fringes  no  longer  applies.  The  best  way  to  determine  layer  com¬ 
positions  and  thicknesses  is  to  compare  the  experimental  rocking  curve  to  simulated 
curves. 

Both  the  kinematical  and  dynamical  theories  of  x-ray  diffraction  have  been  used 
to  simulate  rocking  curves.  The  kinematical  theory  simply  sums  diffracted  intensity 
contributions  from  each  unit  cell,  neglecting  absorption  losses  and  secondary  reflec¬ 
tions.  In  contrast,  the  dynamical  theory  provides  an  exact  solution  to  Maxwell’s 
equations  throughout  the  entire  structure.®®’®®’®^  The  kinematical  theory  is  compu¬ 
tationally  much  simpler  and  has  been  shown  to  provide  adequate  approximations  in 
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Figure  2-3:  HRXRD  (004)  experimental  and  simulated  rocking  curves  for  sample 
1876;  heterostructure  cross-section  is  shown. 

some  cases.  With  the  continuing  improvements  in  the  speed  of  personal  computers, 
however,  it  is  now  reasonable  to  apply  the  dynamical  theory  in  all  cases.  We  have  used 
the  commercial  program  RADS  (Rocking  curve  Analysis  by  Dynamical  Simulation) 
by  Bede  Scientific  Instruments  Limited  throughout  this  thesis. 

We  give  an  example  of  the  use  of  simulations  in  fig.  2.3.  The  structure  examined 
consists  of  (bottom  to  top)  an  (001)  InP  substrate,  a  500  A  layer  of  InyAli_„As,  and  a 
3000  A  layer  of  InxGai_*As.  For  the  simulation,  we  kept  the  layer  thicknesses  fixed  at 
their  nominal  values  (which  we  beheve  to  be  accurate  to  within  10%),  and  adjusted 
the  layer  compositions  until  reasonable  agreement  with  the  experimental  data  was 
obtained.  The  resulting  compositions  are  x  =  0.550  and  y  =  0.50. 


2.2  Crystalline  Quality 

Misfit  dislocations  reheve  lattice  strain,  as  illustrated  in  fig.  1.2.  For  a  constant 
composition,  as  strain  is  reUeved,  the  amount  of  tetragonal  distortion  and  (Ao/a)x 
decrease.  At  the  same  time,  (Aa/a)||  increases  but  (Ao/a),  is  constant.  One  can 
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define  lattice  relaxation,  R,  as:®^ 


We  often  express  as  a  per-cent.  For  a  fully  strained  layer,  R  =  0%;  for  a  completely 
relaxed  layer,  R  =  100%. 

In  theory,  it  should  be  possible  to  determine  critical  layer  thickness  by  measuring 
the  onset  of  relaxation  by  HRXRD.  Multiple  HRXRD  scans  are  required  to  inde¬ 
pendently  measure  both  relaxation  and  composition,  as  will  be  desribed  in  the  next 
section.  If  the  layer  composition  is  already  known,  a  single  (004)  scan  will  give  R.  For 
example,  if  thin  (t<  te)  and  thick  (t>  tc)  layers  are  grown  under  the  same  conditions, 
the  composition  of  the  thin  layer  can  be  measured  by  a  single  (004)  scan  since  R 
wiD  be  zero.  If  we  assume  that  the  thick  layer  has  the  same  composition,  an  (004) 
scan  on  the  thick  layer  will  give  its  R.  Such  an  experiment  is  shown  in  fig.  2.4  for 
the  InAlAs/InP  system.  The  four  epilayers  were  grown  the  same  day  with  identical 
In,  Al,  and  As  cell  temperatures.  The  300,  1000,  and  3000  A  layers  have  nearly 
equal  InP-InAlAs  peak  separations,  suggesting  neghgible  relaxation  (eq.  2.7  applies) 
and  constant  composition  throughout  the  growth  day.  The  peak  separation  for  the 
10,000  A  layer,  however,  is  considerably  reduced. 

Using  the  (Aa/a),  measured  for  the  thinner  layers  and  (2.6),  we  can  calculate 
(Aa/o)||.  Equation  (2.12)  then  gives  R  =  32%.  Apparently,  misfit  dislocations  have 
formed,  relaxing  strain  and  hence  reducing  the  tetragonal  distortion  and  (Aa/a)j.. 
These  measurements  suggest  the  onset  of  relaxation  for  this  composition  occurs  be¬ 
tween  t/tc, MB  =  18  (3000  A  layer)  and  i/tc,MB  =  65  (10,000  A  layer). 

In  similar  studies  of  SiGe/Si®®  and  InGaAs/GaAs,^  the  authors  also  observed  the 
apparent  onset  of  relaxation  at  t*  >>  <c,mb  and  called  t*  the  critical  layer  thickness. 
As  pointed  out  by  Fritz, however,  direct  measurement  of  relaxation  by  HRXRD 
can  give  anomalously  large  values  of  critical  thickness.  A  substantial  number  of 
dislocations  are  required  before  the  parallel  lattice  constant  changes  appreciably  and 
HRXRD  can  detect  the  resulting  change  in  strain.  Hence,  as  will  be  confirmed  in  the 
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Figure  2-4:  HRXRD  (004)  rocking  curves  for  four  Ino.42Alo.B8As  heterostructures 
with  varying  epilayer  thickness.  The  InP-InAlAs  peak  separation  for  sample  1804 
is  reduced  due  to  lattice  relaxation. 


experiments  below,  results  such  as  those  of  fig.  2.4  do  not  yield  a  rehable  measure  of 
the  critical  layer  thickness. 

In  addition  to  peak  separation,  the  width  of  diffraction  peaks  in  HRXRD  contains 
important  information.  The  full-width  at  half-maximum  (FWHM)  of  a  layer  peak  is 
often  cited  as  a  figure  of  merit  for  the  epilayer.  The  FWHM  is  sensitive  to  crystalline 
imperfection.  For  example,  the  presence  of  misfit  dislocations  will  result  in  a  local 
tilting  of  lattice  planes  and  a  broadening  of  the  peak.®®  Even  for  perfect  crystals, 
however,  the  FWHM  is  a  function  of  layer  thickness,  with  thinner  layers  producing 
broader  peaks.  For  example,  in  fig.  2.3  the  500  A  InAlAs  peak  is  much  broader  than 
the  3000  InGaAs  A  peak. 

We  determined  the  theoretical  (004)  FWHM  for  layers  of  InGaAs  and  InAlAs 
by  simulations  using  dynamical  diffraction  theory.®®  We  found  that  the  theoretical 
FWHM  of  both  InGaAs  and  InAlAs  layers  can  be  fit  by  the  function: 


FWHM  =  5.1  + 


1.582  X  10® 
t 


+ 


2.598  X  10® 


(2.13) 
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where  t  is  the  layer  thickness  in  angstroms  and  FWHM  is  in  arc-seconds.  Hence, 
as  a  figure  of  merit  for  crystalline  quality,  we  take  the  ratio  of  the  experimental  to 
the  theoretical  FWHM  for  the  (004)  reflection.  An  additional  measure  of  structural 
quality  is  provided  by  the  presence  or  absence-of  Pendellosung  fringes  in  the  rocking 
curves.  These  fringes  indicate  a  coherent,  high-quality  layer.®® 

For  sets  of  samples  with  constant  thickness  and  changing  composition,  as  mis¬ 
match  increases  we  typically  observe  a  transition  at  which  the  (004)  FWHM  increases, 
presumably  due  to  a  degradation  of  crystalline  quality.*®"®*  In  fig.  2.5,  we  show  an 
example  with  sample  1876  of  fig.  2.3  and  four  other  samples  grown  the  same  day.  The 
five  heterostructures  are  nominally  identical  except  for  the  composition  of  a  3000  A 
In*Gai_,As  layer.  In  addition  to  the  epilayer  and  substrate  peaks,  we  observe  a  series 
of  Pendellosung  fringes  on  four  of  the  samples.  For  samples  1876-1879  (0.477  <  x 
<  0.550),  the  experimental-to-theoretical  FWHM  ratio  is  between  1.1  and  1.5  and 
interference  fringes  are  present.  For  sample  1880  (x  =  0.466),  however,  the  FWHM 
ratio  increases  to  2.4  and  the  fringes  disappear.  The  transition  occurs  when  the  layer 
thickness  exceeds  tc, mb  hy  about  a  factor  of  eight. 

We  observe  a  similar  transition  for  mismatched  layers  of  Ina:Gai_*As  in  compres¬ 
sion  (®  >  0.53).  In  fig.  2.6  we  show  the  rocking  curves  for  four  heterostructures, 
each  consisting  of  a  single  1000  A  layer  of  Ina,Gai_a:As  on  InP,  with  x  varying  from 
0.515  to  0.69.  For  samples  1630,  1632,  and  1633,  the  FWHM  ratios  are  less  than 
1.5  and  Pendellosung  fringes  are  present,  indicating  high  crystalline  quality.  (The 
fringes  are  less  prominent  than  in  fig.  2.5  because  the  InGaAs  in  fig.  2.6  is  thinner 
and  a  second  layer  (InAlAs)  is  not  present.)  For  1634,  however,  the  epilayer  peak  is 
substantially  broadened  and  the  fringes  are  absent,  indicating  a  severe  degradation 
in  layer  quality.  In  a  separate  experiment,  we  grew  a  multi-layer  structure  which 
included  a  1000  A  layer  of  Ino.653Gaio.347As  {i/ic, mb  =  6.1).  The  FWHM  ratio  was 
1.3  and  fringes  were  present,  yielding  a  transition  in  the  range  6  <  t/tc, mb  <  8  for 
1000  A  compressively-strained  layers  of  InGaAs  on  InP. 

We  also  investigated  sets  of  1000  A  Inj,Ali_yAs  epitaxial  layers.  In  fig.  2.7  we  show 
the  rocking  curves  for  five  heterostructures  under  tension  (AlAs-rich)  with  y  varying 
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Figure  2-5:  HRXRD  (004)  rocking  curves  for  five  heterostructures  shown  in  inset. 
The  largest  peaks  (at  0  arc-sec)  are  the  InP  substrate.  The  second-largest  peaks 
(varying  from  -300  to  -1-1200  arc-sec)  are  the  3000  A  InGaAs  layers.  The  broad  peaks 
at  about  -1-500  arc-sec  are  the  500  A  InAlAs  buffer  layers.  t/te.MB  is  the  ratio  of  the 
InGaAs  thickness  to  the  anisotropic  Matthews-Blakeslee  critical  layer  thickness. 


Figure  2-6:  HRXRD  (004)  rocking  curves  for  four  1000  A  layers  of  In*Gai_*As  in 
compression.  The  FWHM  ratios  for  samples  1630,  1632,  and  1633  are  less  than  1.5; 
Pendellosung  fringes  are  also  present  for  these  three  samples. 
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Figure  2-7:  HRXRD  (004)  rocking  curves  for  five  1000  A  layers  of  Inj,Ali_yAs  in 
tension.  Samples  1701,  1705,  and  1702  exhibit  Pendellosung  fringes  and  peak  widths 
within  30%  of  the  theoretical  value. 

from  0.506  to  0.29.  Samples  1701.  1705,  and  1702  exhibit  small  FWHM  ratios  (1.2- 
1.3)  and  Pendellosung  fringes,  indicating  high  crystalline  quality.  For  samples  1703 
and  1704,  the  epilayer  peaks  are  substantially  broadened  and  the  fringes  are  absent, 
indicating  poor  structural  quahty.  For  this  set  the  transition  is  in  the  range  6  < 
t/tc,AfB  <  14. 

We  grew  a  complementary  set  of  1000  A  Inj,Ali_yAs  layers  in  compression  (InAs- 
rich).  The  results  are  shown  in  fig.  2.8.  For  the  first  sample,  3098,  the  lattice 
mismatch  between  the  epilayer  and  substrate  was  too  small  to  produce  two  distinct 
peaks.  (The  thickness  required  for  resolvable  peaks  will  be  discussed  in  section  2.4.) 
For  sample  3099,  the  FWHM  ratio  is  1.3  and  fringes  are  present,  indicating  excellent 
structural  quality.  The  other  three  samples  exhibit  very  broad  epilayer  peaks  with  no 
fringes.  We  can  compare  these  results  to  the  InAlAs  in  tension.  Samples  3100  (fig. 
2.8)  and  1702  (fig.  2.7)  each  have  t/tc  ,mb  =  6,  but  the  crystalline  quality  of  1702 
is  clearly  much  higher.  Similarly,  samples  3102  and  1703  are  14-15  times  tc,KB,  but 
1703  has  a  superior  rocking  curve.  These  results  show  that  the  crystalline  quality 
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Figure  2-8:  HRXRD  (004)  rocking  curves  for  five  1000  A  layers  of  In^Ali-yAs  in 
compression.  Sample  3099  has  Pendellosung  fringes  and  a  peak  width  within  30%  of 
the  theoretical  value. 

of  1000  A  InAlAs  layers  remains  unperturbed  to  larger  mismatches  in  tension  than 
in  compression.  We  have  confirmed  this  finding  for  other  InAlAs  layers  as  well.®® 
The  reason  for  the  differences  in  tension  and  compression  is  not  clear,  but  we  note 
that  Lievin  and  Fonstad,  using  reflection  high-energy  electron  diffraction  (RHEED) 
oscillations,  also  observed  differences  in  the  growth  of  InAlAs  in  compression  and 
tension,  and  related  it  to  differences  in  cation  mobihties.^® 

The  previous  figures  illustrate  the  apparent  degradation  in  layer  quality  when  the 
lattice  mismatch  is  too  large.  In  some  cases,  however,  factors  other  than  mismatch 
result  in  poor  crystalline  quality.  In  fig.  2.9,  we  show  (004)  HRXRD  rocking  curves  for 
samples  2006  and  2010.  Each  sample  consists  of  a  single,  nominally  lattice-matched 
layer  of  InAlAs  on  InP.  The  layer  in  2010  gives  a  sharp  diffraction  peak  and  about  20 
visible  Pendellosung  fringes.  In  contrast,  2006  shows  only  the  InP  peak  with  broad 
shoulders  on  either  side.  We  believe  that  2006  was  grown  at  too  low  a  temperature, 
resulting  in  the  poor  crystalline  quality. 

HRXRD  can  also  be  applied  to  device  heterostructures.  We  wiU  discuss  apphca- 
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Figure  2-9:  HRXRD  (004)  rocking  curves  for  samples  2006  and  2010.  The  sam¬ 
ples  consist  of  single  layers  of  InAlAs  with  approximately  the  same  composition. 
^^2010  exhibits  multiple  fringes  and  a  narrow  epilayer  peak  width,  indicating  excel¬ 
lent  crystalline  quality.  The  quality  of  2006  is  poor,  apparently  because  the  growth 
temperature  was  too  low. 
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Figure  2-10:  HRXRD  (004)  rocking  curves  for  two  MODFETs  with  nominally  identi¬ 
cal  structures,  including  500  A  channels  of  Ino.68Gao.32As.  The  electron  mobilities  of 
#4084  are  9200  and  40,000  cm^/V-sec  at  300  and  77  K,  respectively.  The  mobilities 
of  #4102  are  1900  and  2500  cm^V-sec  at  300  and  77  K,  respectively. 

tions  in  more  detail  in  chapter  5,  but  it  is  worthwhile  to  consider  one  example  here.  In 
fig.  2.10,  we  show  the  rocking  curves  for  two  modulation-doped  field-effect  transistors 
(MODFETs).  The  structures  consist  of  several  layers  of  both  InGaAs  and  InAlAs 
(see  chapter  4  for  details).  All  the  layers  are  nominally  lattice-matched  except  for 
500  A  layers  of  Ino.68Gao.32As  which  serve  as  the  channel  of  the  device.  We  measured 
the  electron  mobibty  of  both  structures  at  300  K  and  77  K.  For  sample  4084,  we 
obtained:  fisooK  =  9200  cm^/V-sec  and  firrK  =  40,000  cm^/V-sec.  These  mobilities 
are  reasonable  for  this  structure.*®-'^"  The  FWHM  of  the  channel  peak  is  close  to  the 
theoretical  value,  and  fringes  are  present.  In  contrast,  sample  4102  has  low  mobility 
with  fi300K  =  1900  cmVV-sec  and  fmK  =  2500  cmVV-sec,  and  a  poor  rocking  curve. 
We  are  not  certain  what  went  wrong  in  the  growth  of  4102.  Possibihties  include  in¬ 
correct  growth  temperature  and  contamination  in  the  MBE.  The  correlation  between 
electrical  and  structural  characteristics  suggests  that,  in  some  cases,  x-ray  data  can 
be  used  to  decide  whether  to  process  a  device  heterostructure. 

The  growth  of  poor  quality  epilayers  due  to  incorrect  growth  temperature  or  other 
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Figure  2-11:  HRXRD  (004)  peak  width  versus  InGaAs  or  InAlAs  thickness  for  layers 
in  compression  and  tension.  Samples  are  coded  based  upon  the  ratio  of  layer  thickness 
to  anisotropic  Matthews-Blakeslee  critical  thickness.  Solid  line  is  from  simulations 
using  dynamical  diffraction  theory  (see  eq.  2.13).®® 

unknown  reasons  was  relatively  rare,  occurring  less  than  10%  of  the  time.  We  elim¬ 
inated  such  samples  from  our  data  set.  In  fig.  2.11,  we  show  the  (004)  HRXRD 
epitaxial  peak  width  as  a  function  of  thickness  for  the  remainder  of  our  sample  set 
which  includes  InGaAs  and  InAlAs  samples  in  both  compression  and  tension.  The 
samples  are  coded  based  upon  the  ratio  of  the  layer  thickness  to  The  theoret¬ 

ical  values  of  FWHM  for  perfect  layers  (represented  by  the  line)  are  a  lower  limit  to 
our  results,  with  experimental  FWHM’s  approximately  equal  to  the  theoretical  values 
for  layers  ranging  from  200  A  to  10,000  A.  We  observe  that  it  is  not  necessary  for 
the  layer  to  be  thinner  than  tc.MB  in  order  to  achieve  a  FWHM  ratio  close  to  one. 

This  result  is  more  clearly  illustrated  by  fig.  2.12  in  which  we  plot  the  FWHM 
ratio  versus  the  ratio  of  thickness  and  for  InGaAs  (a)  and  InAlAs  (b).  For 

<  3  the  FWHM  ratio  is  always  between  0.9  and  1.8,  independent  of 
We  observe  a  transition  region  for  3  <  iftc,MB  <  9,  with  the  FWHM  ratio  varying 
from  1.0  to  10.  For  tjtc,MB  >  9,  the  ratio  ranges  from  2  to  60.  In  fig.  2.12  the 
samples  are  coded  for  layers  in  tension  or  compression.  For  InGaAs,  there  is  no  clear 
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difference  between  the  two.  For  InAlAs,  we  observe  the  same  trend  illustrated  in  figs. 
2.7  and  2.8:  higher  crystalline  quality  for  layers  in  tension  than  in  compression  (at 
constant  t/tc,MB)- 

Our  results  for  InAlAs  are  in  contrast  to  a  recent  report  by  Tournie  et  al7^  They 
investigated  nearly  lattice-matched  1.3  fim  layers  of  InAlAs  on  InP.  For  substrate 
temperatures  of  480-520" C,  the  (004)  HRXRD  FWHM’s  were  240-440  arc-sec,  com¬ 
pared  to  a  theoretical  value  of  only  18  arc-sec.  They  achieved  FWHM’s  of  less  than 
50  arc-sec  only  when  a  high  substrate  temperature  (600°C)  or  a  superlattice  buffer 
was  used.  As  shown  in  fig.  2.12b,  however,  when  the  mismatch  is  not  too  large  we 
consistently  obtain  FWHM’s  close  to  the  theoretical  value  for  InAlAs  grown  directly 
on  InP  at  460-510°C.  A  possible  explanation  for  this  discrepancy  is  errors  in  the 
measurement  of  growth  temperature. 

In  fig.  2.13,  we  compare  the  relaxation  and  crystalline  quality  criteria  by  plotting 
thickness  versus  lattice  mismatch  for  InGaAs  (a)  and  InAlAs  (b).  We  separate  the 
samples  into  three  distinct  groups.  For  i/tc,MB  <  3,  the  FWHM  ratio  appears  to  be 
randomly  scattered  between  0.9  and  1.8  (see  fig.  2.12).  Hence,  we  choose  1.8  as  a 
cut-off  point  and  consider  layers  with  a  smaller  FWHM  ratio  to  be  of  high  crystalline 
quality.  For  all  samples  with  R  >  10%,  the  FWHM  ratio  was  greater  than  1.8.  For 
all  samples  with  a  FWHM  ratio  less  than  1.8,  the  layer  relaxation  was  less  than  10%. 
The  remaining  samples  fall  into  a  third  intermediate  group  with  R  <  10%  and  the 
FWHM  ratio  greater  than  1.8. 

In  fig.  2.13,  we  include  curves  R  separating  layers  based  upon  relaxation,  with 
a  cut-off  at  R  =  10%.  (We  use  R  to  refer  to  an  average  relaxation,  as  discussed  in 
section  2.3.)  The  minimum  detectable  relaxation  depends  upon  layer  thickness  (which 
influences  the  width  and  intensity  of  the  peak),  experimental  background  noise,  and 
lattice  mismatch.  For  most  samples,  the  smallest  relaxation  we  can  measure  is  5-10%. 
Our  empirical  curves  R,  obtained  with  this  R  =  10%  criterion,  are  actually  similar  to 
the  calculated  People-Bean  critical  layer  thickness^^  and  estimates  from  HRXRD. 

We  also  include  curves  C  in  fig.  2.13.  They  separate  the  samples  based  upon 
crystalline  quality  as  determined  by  the  FWHM  ratio  with  a  cut-off  of  1.8.  A  com- 
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Figure  2-12;  Ratio  of  experimental  to  theoretical  HRXRD  (004)  peak  width  versus 
ratio  of  thickness  to  anisotropic  Matthews-Blakeslee  critical  thickness  for  (a)  InGaAs 
and  (b)  InAlAs  layers.  Samples  are  coded  for  compression  and  tension. 
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Figure  2-13:  Thickness  versus  relaxed  lattice  mismatch  for  a)  InGaAs  and  b)  In- 
AlAs  in  both  compression  and  tension  on  InP.  Samples  are  coded  by  the  HRXRD 
peak-width  ratio  and  lattice  relaxation.  Curves  C  and  R  are  empirical,  showing  the 
degradation  of  crystalline  quality  and  significant  lattice  relaxation,  respectively. 
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parison  of  C  and  R  confirms  that  the  FWHM  ratio  is  more  sensitive  to  crystalline 
imperfection  than  direct  measurements  of  lattice  relaxation.  From  a  practical  point 
of  view,  we  also  note  that  measurement  of  the  FWHM  ratio  is  much  faster,  requiring 
only  a  single  (004)  scan  rather  than  a  complete  set  of  asymmetric  scans  (see  section 
2.3). 

The  FWHM  of  an  epitaxial  layer  can  be  broadened  by  a  variety  of  factors.  One,  as 
mentioned  earher,  is  the  finite  thickness  effect.  We  removed  this  effect  by  examining 
the  FWHM  ratio.  Another  potential  broadening  factor  is  the  warpage  of  lattice  planes 
resulting  from  strain.^®  To  test  for  warpage  effects,  we  performed  a  series  of  scans  with 
varying  aperture  size  on  the  x-ray  collimator.  We  used  a  sample  with  a  FWHM  ratio 
of  2.4.  The  results  showed  that  the  epitaxial  peak  width  was  independent  of  the  x-ray 
spot  size.  From  this,  we  conclude  that  lattice  warpage  does  not  have  a  significant 
effect  on  the  epitaxial  peak  width  for  our  samples.  Improper  instrument  alignment 
can  also  broaden  x-ray  peaks.  We  have  achieved  InP  substrate  peaks  as  narrow  as  9 
arc-sec,  the  theoretical  minimum,®^  implying  proper  alignment. 

Having  eliminated  the  above  sources  of  peak  broadening,  we  conclude  that  the 
peak  broadening  observed  for  samples  with  FWHM  ratios  greater  than  1.8  is  primar¬ 
ily  the  result  of  crystalline  imperfection.  Traditionally,  such  broadening  in  strained 
layers  has  been  attributed  to  misfit  dislocations.  Recent  work  has  suggested  the 
possibility  of  3D  growth  without  misfit  dislocations.®®’®^  Such  growth  could  lead  to 
imperfections  (e.y.  stacking  faults)  which  might  also  broaden  an  HRXRD  peak  and 
eliminate  Pendellosung  fringes.  Hence,  we  may  be  observing  the  effects  of  misfit 
dislocations,  3D  growth,  or  both.  We  will  return  to  this  topic  in  chapter  6. 

We  have  shown  that  the  epilayer  FWHM  ratio  and  Pendellosung  fringes  can  be  a 
sensitive,  non-destructive  means  to  assess  crystalline  quality  of  mismatched  epilayers. 
Our  results  do  not  necessarily  imply  the  total  absence  of  misfit  dislocations  for  samples 
with  Pendellosung  fringes  and  FWHM  ratios  near  unity.  Although  we  cannot  rule 
out  the  possibility  of  dislocations  in  layers  with,  for  example,  t/tc, mb  =  3,  such  layers 
do  exhibit  high  crystalline  quality  and  may  be  useful  for  many  device  applications. 
Knowledge  of  the  degree  of  strain  (relaxation)  is  also  important  for  both  active  and 
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buffer  layers  in  device  heterostructures.  Our  results  show  that  layers  remain  fully 
strained  well  beyond  tc.MB- 

According  to  some  reports,  growth  temperature  may  play  an  important  role  in 
the  relaxation  of  mismatched  epitaxial  layers.^®’ ^  The  temperature  range  over  which 
we  could  obtain  high-quality  InGaAs  and  InAlAs  was  only  about  50"C.  In  addition, 
the  uncertainty  in  our  temperature  measurements  was  about  30®C.  Hence,  we  were 
not  able  to  make  a  systematic  study  of  growth  temperature  effects  for  our  MBE 
layers.  We  did,  however,  examine  InAlAs/InP  samples  grown  by  MOCVD.'^  These 
showed  FWHM  ratios  close  to  1.0  and  Pendellosung  fringes  for  1000  A  layers  of 
InyAli-yAs,  y=0.41  and  0.42,  with  =  7.  The  InAlAs  was  grown  at  625°C, 

over  lOO^C  higher  than  our  MBE  samples.  This  result  suggests  that  high  crystaUine 
quality  beyond  tc,MB  is  not  restricted  to  our  growth  conditions.  We  also  performed 
post-growth  annealing  experiments  on  our  MBE  layers  as  described  in  chapter  4. 


2.3  Lattice  Relaxation  and  Orthorhombic 
Distortion 

In  the  previous  section,  we  showed  that  when  mismatch  is  too  large,  the  HRXRD 
epilayer  peak  broadens.  For  some  samples,  we  found  that  this  broadening  was  a 
function  of  azimuthal  angle. Fig.  2.14  illustrates  the  geometry  of  our  measurements. 
The  angle  of  incidence  is  fixed  at  the  Bragg  angle.  The  sample  is  rotated,  and 
measurements  are  taken  at  different  azimuthal  angles,  a.  At  a  =  0",  the  incident 
x-ray  beam  is  orthogonal  to  the  primary  flat  and  parallel  to  the  [110]  direction.  At 
a  =  90°,  the  beam  is  parallel  to  the  [110]  direction.  The  array  of  misfit  dislocations 
will  be  explained  below. 

In  fig.  2.15,  we  show  the  (004)  rocking  curves  at  a  =  0°  and  90°  for  sample  1503, 
a  1700  A  layer  of  Ino.6oGao.4oAs.  For  this  layer,  =  5.2  and  the  theoretical 

FWHM  is  99  arc-sec.  The  layer  peak  is  broadened  in  both  directions,  but  the  FWHM 

^MOCVD  samples  were  supplied  by  Noren  Pan  of  Raytheon  Corporation 
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Misfit  \/ 
Dislocations 


Figure  2-14:  Geometry  of  x-ray  measurements;  a  is  the  azimuthal  angle  (angle  of  ro¬ 
tation)  of  the  projection  of  the  beam  on  the  sample  with  respect  to  the  [110]  direction. 
An  asymmetric  array  of  misfit  dislocations  is  also  shown. 

is  substantially  larger  at  a  =  0"  (683  arc-sec)  than  at  a  =  90*’  (486  arc-sec).  To 
confirm  this  effect,  we  made  (004)  measurements  at  eight  different  azimuthal  angles 
for  sample  1503;  the  results  are  plotted  in  fig.  2.16.  The  data  approximately  follow  a 
cosine  law,  with  FWHM  minima  at  a  =  90  and  270®.  We  obtained  similar  results  for 
samples  of  InGaAs  in  tension  as  well  as  InAlAs  in  both  tension  and  compression.**’^^ 
In  aJl  cases  with  significant  FHWM  variations,  the  minima  occurred  when  the  incident 
beam  was  pairaillel  to  the  [110]  and  [110]  directions. 

In  zinc-blende  semiconductors  such  as  InGaAs,  the  [lIO]  and  [IlO]  directions  are 
equivalent.  The  [110]  and  [110]  directions  are  also  equivalent,  but  [110]  is  not  equiv¬ 
alent  to  [110].^*  In  fact,  several  papers  report  an  asymmetry  in  misfit  dislocation 
density  for  mismatched  III-V  semiconductors,  with  a  greater  density  in  one  <  110  > 
direction  than  in  the  orthogonal  direction  for  (001)  substrates.*®’^*"”  Misfit  disloca¬ 
tions  cause  a  local  tilting  of  lattice  planes.  These  tilted  planes  will  satisfy  the  Bragg 
condition  at  angles  that  differ  slightly  from  the  Bragg  smgle  for  a  dislocation-free 
layer.  Based  on  these  facts,  we  can  now  explain  the  data  of  figs.  2.15  and  2.16  as 
follows.  If  the  number  of  dislocations  along  the  [lIO]  direction  exceeds  the  number 


Figure  2-15:  HRXRD  (004)  rocking  curves  for  sample  1503,  showing  different  epilayer 
peak  widths  at  a  =  0  and  90”. 


[110]  [110]  [110]  [110]  [110] 


Figure  2-16:  HRXRD  (004)  epilayer  peak  width  as  a  function  of  azimuthal  angle  for 
sample  1503.  SoHd  line  is  a  least-squares  fit  to  a  cosine  function. 
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GLANCING-EXIT 


Figure  2-17:  Geometry  of  glancing-exit  and  glancing-incidence  asymmetric  HRXRD 
measurements.  9b  is  the  Bragg  angle  and  (f)  is  the  angle  between  the  surface  normal 
and  the  diffracting  planes. 

along  [110]  (as  schematically  illustrated  in  fig.  2.14),  the  peak  width  observed  along 
[110]  (a  =  0")  should  be  greater  than  that  observed  along  [110]  (a  =  90"),  as  our 
experiments  indicate.  Such  an  anisotropic  dislocation  density  distribution  was  ob¬ 
served  with  HRXRD  by  Grundmann  et  al.  for  InGaAs/GaAs  heterostructures  under 
compression.*®  Our  work  shows  that  similar  behavior  occurs  for  InGaAs/InP  and 
InAlAs/InP  under  compression  or  tension. 

Since  misfit  dislocations  relieve  epilayer  strain,  an  asymmetry  in  misfit  dislocation 
density  should  result  in  different  parallel  lattice  mismatches  in  the  orthogonal  [110] 
and  [110]  directions.  We  can  measure  the  parallel  mismatches  with  HRXRD  by 
using  asymmetric  reflections.  In  an  asymmetric  reflection,  the  diffracting  plane  is  not 
parallel  to  the  surface.  Hence,  to  obtain  a  Bragg  reflection,  the  angle  of  incidence 
must  be  equal  to  ±  where  <p  is  the  angle  between  the  surface  and  the  diffracting 
plane.*®  The  detector  is  at  26b,  as  in  the  case  of  symmetric  reflections.  The  two 
asymmetric  cases  are  illustrated  in  fig.  2.17  and  referred  to  as  glancing-incidence  and 
glancing-exit  reflections. 
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For  asymmetric  reflections,  the  measured  peak  separation  between  the  substrate 
and  layer  peaks  is  no  longer  simply  the  difference  in  Bragg  angles.  Instead,  it  is  given 
by: 


M  =  Mb  ±  A(t>  (2.14) 

where  is  the  tilt  of  the  substrate  with  respect  to  the  epilayer.*^  In  the  case  of 
tetragonal  distortion,  equivalent  asymmetric  planes  of  the  substrate  and  layer  are  not 
parallel.  Hence,  A(f>  is  non-zero  and  contributes  to  the  Bragg  peak  separation.  The 
+  and  -  signs  refer  to  glancing-incidence  and  glancing-exit  reflections,  respectively. 
The  terms  in  (2.14)  can  be  written  in  terms  of  lattice  constant  differences:®^ 


(2.15) 


(2.16) 


In  the  case  of  symmetric  reflections  such  as  the  (004)  reflection  from  (001)  substrates, 
(f>  =  0°  and  (2.14)-(2.16)  reduce  to  (2.4).  Hence,  the  perpendicular  mismatch  can 
be  calculated  directly  from  an  (004)  rocking  curve,  as  we  did  in  section  2.1.  For  the 
general  case,  we  combine  (2.14),  (2.15),  and  (2.16): 


'j  (— sin^  (^tan  ±  sin  ^cos  ^) -F  ^  (— cos^i^tan^B  T  sin^cos^) 

\a/||  \  a  /  ± 

(2.17) 

where  the  upper  signs  are  for  glancing-exit  and  the  lower  signs  are  for  glancing- 
incidence  reflections. 

We  made  measurements  of  both  (115)  and  (224)  asymmetric  reflections.  Since 
the  substrates  are  (001),  it  is  straightforward  to  calculate:  <f>n6  =  15. 79°  and  ^224  = 
35.26°.  Using  (2.3)  to  calculate  the  Bragg  angles,  (2.17)  reduces  to  the  following  four 
equations  (for  InP  substrates  only). 

(115)  glancing-incidence: 
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(115)  glancing-exit: 


=  -0.601  (-^)^  -  0.331 
(224)  glancing-incidence: 


A9  =  -1.031  +  0.191 


(2.19) 


(2.20) 


(224)  glancing-exit: 

A^= -0.089  ^^)  _  0.751  (2.21) 

In  the  case  of  simple  tetragonal  distortion,  two  measurements  are  required  to  ob¬ 
tain  both  (Aa/o)_L  and  (Ao/o)||.  For  example,  we  could  measure  the  (115)  glancing- 
exit  and  glancing-incidence  reflections.  If  instead  (Ac/a)x  is  known  from  an  (004) 
measurement,  then  only  one  asymmetric  measurement  is  required  to  determine  ( Aa/a)|| 
In  this  case,  it  is  preferable  to  measure  one  of  the  glancing-exit  reflections  because  A6 
will  be  more  heavily  weighted  by  the  parallel  mismatch  according  to  (2.18)-(2.21). 

In  fig.  2.18  we  show  (224)  rocking  curves  for  sample  1245,  8500  A  of  InGaAs 
under  compression.  Both  the  CuKai  and  CuKa2  components  are  visible  for  the  InP 
peak.  The  peak  separation,  A9,  is  different  for  the  curves  measured  at  a  =  0  and  90®. 
Instead  of  a  single  (Aa/o)||,  we  have  different  values  in  the  two  orthogonal  directions 
on  an  (001)  substrate,  (Aa/a)||[no]  and  (Aa/a)||[iio]. 

To  calculate  the  parallel  mismatches,  we  average  scans  separated  by  180® 
[(Aa/o)||[iio]  is  calculated  from  the  average  of  A9(^224)  (a  =  0®)  and  A9(224)  («  = 
180®);  (Aa/o)||[iio]  is  calculated  from  the  average  of  A9^224)  («  =90®)  and  A9(224)  (a 
=  270®)]  to  eliminate  differences  in  A^  caused  by  possible  epilayer  tilt.  We  also  note 
that  the  epilayer  peak  is  sharper  for  a  =  90®  than  for  0®,  presumably  due  to  the  same 
reason  as  the  (004)  peak  width  variation:  an  asymmetry  in  dislocation  density. 


RELATIVE  BRAGG  ANGLE  (arc-seconds) 


Figure  2-18:  HRXRD  rocking  curves  for  (224)  glancing-exit  reflections  of  sample  1245 
at  azimuthal  angles  of  0°  and  90".  The  difference  in  peak  separation  for  the  two  curves 
indicates  asymmetric  strain  in  the  epilayer. 

An  epilayer  lattice  is  now  described  by  three  mismatches:  (Ao/a)x,  (Aa/a)||[iio], 
and  (Aa/a)||[iio]-  The  conventional  expression  for  the  relaxed  lattice  mismatch,  (2.6), 
must  be  modified.  The  new  expression  is: 


where  u  is  Poisson’s  ratio.  If  (Aa/a)j.  /  (Aa/a)||[iio]  ^  (Aa/a)||[iio],  the  distortion 
of  the  lattice  is  orthorhombic.  For  example,  the  data  for  sample  1245  of  fig.  2.18 
indicates  that  (Aa/a)i  =  8.33  x  10-^  (Aa/a)||[no]  =  2.94  x  lO'^,  and  (Aa/a)||[aio] 
=  1.88  X  10"^.  Using  v  =  1/3,  we  find  (Ao/a)r  =  5.37  x  10~^.  From  (2.8),  the 
composition  is  Ino.eioGa^o.SQoAs. 

Whereas  the  degree  of  relaxation  is  defined  by  a  single  parameter  in  the  case  of 
purely  tetragonal  distortion,  two  parameters  are  needed  to  characterize  orthorhombic 
distortion.  One  can  define: 
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(2.23) 


and 


a  /liriinl 


We  also  define  the  average  relaxation,  R: 


(2.24) 


g  =  (2.25) 

For  symmetric  relaxation,  R[uo]  =  -^[110])  ®<1-  (2.25)  reduces  to  (2.12),  the  conventional 
definition  of  relaxation.  For  sample  1245  we  find:  iZ[iio]  =  55%,  i2[iio]  =  35%,  and  R 
=  45%. 

For  most  samples,  we  determined  R  from  an  (004)  and  four  (115)  or  (224)  HRXRD 
measurements,  as  discussed  above.  R  can  also  be  calculated  from  only  the  symmetric 
(004)  rocking  curve  if  the  composition  is  known.  The  composition  is  estimated  from 
other  samples  grown  the  same  day  and  the  measured  activation  energy  of  the  A1  or 
Ga  cell.  (The  In  cell  temperature  is  held  constant  in  our  work.)  We  use  this  method 
for  some  thin  epilayers  {t  <  lOOOA)  because  of  the  weak  and  broad  asymmetric  layer 
peaks  in  HRXRD.  We  show  each  partially  relaxed  layer  on  a  plot  of  thickness  versus 
mismatch  for  a)  InGaAs  and  b)  InAlAs  in  figure  2.19.  The  value  of  R  is  given  for  each 
sample.  The  uncertainties  in  R  are  5-10%.  We  do  not  include  samples  with  R  <  5%. 

In  general,  we  observe  the  expected  trends  in  fig.  2.19.  At  a  constant  mismatch, 
R  generally  increases  with  increasing  thickness.  For  a  fixed  thickness,  R  usually 
increases  with  increasing  mismatch  (positive  or  negative).  We  note  that  one  of  the 
InAlAs  layers  is  114%  relaxed.  In  this  case,  (Aa/a)|||jj,j  and  (Aa/a)|||jjgj  are  slightly 
larger  than  (Aa/a)j_.  A  similar  “overrelaxed”  InGaAs  layer  was  reported  by  Chu  et 
a/.**  Apparently,  more  dislocations  are  introduced  during  growth  than  are  required 
to  fully  relax  the  strain. 

In  section  2.2,  we  showed  that  the  crystaUine  quality  of  epilayers  begins  to  degrade 
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Figure  2-19:  Thickness  versus  lattice  mismatch  for  a)  InGaAs  and  b)  InAIAs.  The 
average  lattice  relaxation,  R,  is  given  for  each  sample.  Samples  with  R  <  5%  are  not 


before  significant  lattice  relaxation.  Hence,  all  the  partially  relaxed  samples  in  fig- 
2.19  have  a  high  density  of  structural  defects,  making  them  unsuitable  as  active  layers 
in  most  device  applications.  There  are  some  structures,  however,  in  which  the  active 
regions  of  a  device  are  grown  on  top  of  a  fuUy  relaxed  buffer  layer.®®’ Only  three 
of  our  samples  are  more  than  90%  relaxed.  In  each  case,  the  layer  is  over  50  times 
thicker  than  Hence,  very  thick  layers  are  generally  required  for  complete  strain 

relaxation. 

We  have  observed  orthorhombic  distortion  in  all  of  our  heterostructures  with  R  > 
10%.  In  fig.  2.20,  we  plot  the  difference  in  i2[iio]  and  R[iio]  versus  for  all 

samples  on  which  a  complete  set  of  asymmetric  x-ray  scans  were  collected.  The 
samples  are  coded  by  material  and  strain.  For  the  two  samples  with  <  15, 

the  difference  in  relaxation  is  zero  within  the  experimental  error  bars.  The  average 
relaxation  is  also  negligible  for  these  two  samples.  For  all  the  other  samples,  (i2[iio]- 
■R[iIo])  >  0.  We  conclude  that  these  partially  relaxed  layers  are  orthorhombically 
distorted,  with  maximum  strain  rehef  in  the  [110]  direction,  implying  a  greater  density 
of  misfit  dislocations  running  parallel  to  the  [110]  direction^  for  layers  of  both  InGaAs 
and  InAlAs,  in  either  tension  or  compression. 

The  magnitude  of  the  relaxation  asymmetry  cannot  be  predicted  from 
For  example,  at  tjtc,MB  =  20,  the  values  of  {R[uo]-R[iw])  are  5,  9,  19,  31,  and  37%. 
The  extreme  values  of  5  and  37%  were  both  measured  for  samples  of  InGaAs  in 
compression.  Since  the  relaxation  asymmetry  is  often  substantial  and  apparently 
cannot  be  predicted,  a  complete  set  of  asymmetric  x-ray  scans  must  be  measured  to 
determine  the  composition  of  a  partially  relaxed  epilayer. 

We  note  that  for  many  partially  relaxed  layers  on  which  asymmetric  rocking  curves 
revealed  orthorhombic  distortion,  the  (004)  FWHM  did  not  vary  as  a  function  of 
azimuthal  angle.  We  speculate  that  once  the  dislocation  density  reaches  a  certain 
value  in  both  directions,  the  FWHM  saturates  and  its  angular  dependence  disap¬ 
pears.  Hence,  the  absence  of  FWHM  variations  with  a  does  not  imply  the  absence 

tPor  the  Sumitomo  wafers  used  in  this  study,  the  majority  of  dislocations  are  parallel  to  the 
major  flat  and  orthogonal  to  the  oval  defects. 
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Figure  2-20:  Difference  in  relaxation  in  the  [110]  and  [110]  directions  as  a  function 
of  the  ratio  of  thickness  to  Matthews-Blakeslee  critical  layer  thickness.  Most  of  the 
samples  exhibit  substantial  orthorhombic  distortion.  The  error  bar  shown  is  typical 
of  the  estimated  error  for  all  samples. 


of  orthorhombic  distortion. 

In  fig.  2.21  we  show  a  schematic  view  of  lattice  distortion  suggested  by  our  data. 
We  assume  the  epilayer  has  a  larger  lattice  constant  than  the  substrate.  For  layers 
less  than  a  critical  thickness  (t  <  tc),  tetragonal  distortion  results  (fig.  2.21a).  For  t 
>  tc,  misfit  dislocations  form  in  an  asymmetric  pattern,  partially  relaxing  the  strain 
and  resulting  in  orthorhombic  distortion  (b,c).  For  t  >>  tc,  the  strain  is  fully  relieved 
by  dislocations  and  the  epilayer  symmetry  becomes  cubic  (d). 

Although  there  are  many  reports  showing  an  asymmetry  in  misfit  dislocation 
density  on  (001)  zinc-blende  substrates,  there  is  some  confusion  in  the  literature 
regarding  which  direction  ([110]  or  [110])  has  the  maximum  density.^^’’^®-®®  Using  two 
different  techniques,  we  have  shown  that  the  majority  of  dislocations  are  parallel  to 
the  [lIO]  direction  for  both  InGaAs  and  InAlAs  in  compression  or  tension.  In  chapter 
6,  we  will  discuss  the  reason  for  the  asymmetry  in  dislocation  density.  We  will  also 
show  that  relaxation  may  be  more  complex,  with  features  such  as  elongated  islands 
contributing  to  asymmetry. 
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Figure  2-21:  Schematic  view  of  lattice  distortion  as  a  function  of  degree  of  relaxation 
for  an  epilayer  with  a  larger  lattice  constant  than  the  substrate:  a)  t  <  b)  and  c) 
t  >  tc,  d)  t  >>  tc. 

2.4  Optimal  Thickness  for  Composition 
Measurements 

As  a  result  of  the  orthorhombic  distortion  described  in  the  previous  section,  at  least 
five  HRXRD  scans  are  required  to  determine  (Aa/a)r  for  partially  relaxed  InGaAs 
and  InAlAs  on  InP.  Since  a  single  scan  typically  takes  an  hour  or  more,  measuring 
the  composition  of  a  partially  relaxed  layer  is  very  time-consuming.  In  addition, 
some  commercial  HRXRD  systems  are  only  designed  to  measure  (004)  symmetric 
reflections.  Hence,  when  calibrating  composition,  the  growth  of  relaxed  epitaxial 
layers  should  be  avoided  to  eliminate  the  need  for  asymmetric  measurements. 

Strain  relaxation  is  illustrated  by  the  four  (004)  rocking  curves  shown  in  fig. 
2.4.  As  discussed  in  section  2.2,  misfit  dislocations  formed  in  siimple  1804,  relaxing 
strain  and  hence  reducing  the  tetragonal  distortion  and  A6.  Hence,  the  (004)  peak 
separation  and  the  assumption  of  coherency  would  give  the  correct  composition  for 
samples  1808,  1806,  and  1805,  but  an  incorrect  value  for  sample  1804.  For  this  sample 
set,  it  is  clear  that  significant  relaxation  does  not  occur  until  tltc,MB  >  IS* 
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Figure  2-22:  Thickness  versus  lattice  mismatch  for  150  samples  of  InGaAs  and  InAlAs 
on  InP.  Samples  are  coded  by  lattice  relaxation.  Both  composition  and  relaxation 
were  measured  by  HRXRD.  Note  that  the  lattice  relaxation  is  always  less  than  10% 
for  layers  less  than  five  times  the  Matthews-Blakeslee  critical  thickness. 

We  can  confirm  the  findings  from  fig.  2.4  by  examining  our  entire  sample  set, 
including  both  InGaAs  and  InAlAs  in  compression  and  tension.  We  plot  the  thickness 
of  each  sample  as  a  function  of  mismatch  in  fig.  2.22.  The  samples  are  coded  by 
average  lattice  relaxation,  R.  We  conclude  that  for  these  materials,  a  layer  can 
exceed  ic,MB  by  at  least  a  factor  of  five  with  relaxation  of  less  than  10%. 

The  possibility  of  layer  relaxation  puts  an  upper  limit  on  the  thickness  of  calibra¬ 
tion  layers.  If  layers  are  too  thin,  however,  they  cannot  be  detected  by  HRXRD.  We 
illustrate  this  in  fig.  2.23  which  includes  simulated  rocking  curves  for  single  layers  of 
Ino.6oGao.4oAs  on  InP.  The  layer  thicknesses  are  100,  200,  500,  and  1000  A.  Epilayer 
peaks  are  visible  for  the  200-1000  A  layers.  As  a  result  of  the  tail  of  the  InP  substrate 
peak,  however,  no  peak  is  visible  for  the  100  A  layer.  Hence,  for  this  mismatch,  the 
lower  limit  for  layer  thickness  is  between  100  and  200  A.  For  smaller  mismatches, 
thicker  layers  will  be  required  because  the  substrate  peak  will  be  closer  to  the  layer 
peak. 

Based  upon  simulations  such  as  fig.  2.23,  we  determined  the  minimum  layer 
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Figure  2-23:  Simulated  HRXRD  rocking  curves  for  Ino.6oGao.4oAs/InP  heterostruc¬ 
tures  with  different  epilayer  thicknesses.*® 

thickness  required  to  produce  a  resolvable  peak.  We  plot  the  results  for  Ina,Gai_j,As 
layers  on  InP  as  a  function  of  mismatch  in  fig.  2.24.  To  within  5%,  the  minimum 
thicknesses  are  the  same  for  InGaAs  and  InAlAs  in  tension  or  compression.  If  the 
maximum  thickness  were  limited  by  the  range  of  acceptable  thicknesses  would 

be  quite  small.  Fortunately,  as  shown  in  figs.  2.4  and  2.22,  layers  can  exceed  tc.MB 
without  substantial  relaxation.  We  include  our  conservative  upper  limit  of  htc,MB  in 
fig.  2.24.  We  use  <c,mb  for  InGaAs  in  compression,  but  the  values  are  similar  for 
the  other  three  cases  (see  appendix  C).  The  other  items  in  fig.  2.24  will  be  discussed 
below. 

The  procedure  used  to  obtain  the  minimum  resolvable  peak  (MRP)  simulations 
line  in  fig.  2.24  imphcitly  assumes  “perfect”  layers  and  a  measurement  system  with 
infinite  resolution  (dynamic  range).  Although  defect-free  single  crystals  may  be  grown 
in  some  cases,  there  is  always  an  experimental  background  count-rate  which  can  make 
it  more  difficult  to  identify  layer  peaks.  Hence,  it  is  important  to  experimentally  verify 
the  MRP  criteria.  We  did  so  by  searching  our  sample  set  and  selecting  those  layers 
which  were  reasonably  close  to  the  MRP  line.  We  plot  these  samples  as  data  points  on 


61 


Figure  2-24:  Epilayer  thickness  versus  mismatch  showing  the  minimum  resolvable 
peak  (MRP)  criteria,  the  thickness  at  which  the  error  in  composition  due  to  layer 
peak  shift  is  0.005,  and  the  10%  relaxation  threshold  of  5tc,MB-  Experimental  data 
points  are  also  shown.  The  estimated  error  in  composition,  Aajj,  is  less  than  0.01  for 
samples  in  region  A.  In  region  B,  Air  <  0.03. 
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fig.  2.24.  The  rocking  curve  for  each  sample  had  a  distinct  layer  peak.  As  an  example, 
we  consider  sample  2010  in  fig.  2.9.  This  sample  was  a  typical  single-layer  calibration 
growth.  We  observe  distinct  peaks  from  both  the  substrate  and  the  epilayer  as  well  as 
interference  fringes.  The  peak  separation,  A6,  is  -|-84  arc-sec.  Using  (2.4),  (2.6),  and 
(2.9),  we  obtain  a  composition  of  Ino.5i6Alo.484As.  From  simulations,  the  minimum 
thickness  is  2000  A  for  a  resolvable  peak  at  this  mismatch.  The  experimental  fringe 
spacing  gives  a  thickness  of  3700  A  for  sample  2010.  Considering  the  clear  peak 
separation  in  fig.  2.9,  we  believe  that  layers  of  the  same  composition  would  produce 
resolvable  peaks  for  thicknesses  less  than  3700  A.  Based  upon  our  data  points,  we 
include  a  dashed  MRP-expt  line  in  fig.  2.24.  Epilayer  peaks  should  by  resolvable  by 
HRXRD  for  points  above  this  Hne. 

Our  experimental  results  and  simulations  suggest  that  the  composition  of  layers 
lying  above  the  MRP-expt  line  and  below  the  5tc,MB  line  can  be  measured  by  a  single 
HRXRD  scan.  We  now  discuss  six  sources  of  error  and  the  resulting  accuracy  for 
such  measurements. 

1.  Lattice  constants  and  Poisson  ratios:  Estimates  of  composition  from 
HRXRD  measurements  require  knowledge  of  lattice  parameters  of  the  substrate  (InP) 
and  end-members  of  the  ternary  alloy  (In As,  GaAs,  AlAs),  as  well  as  Poisson’s  ratio 
for  the  alloy.  Until  recently,  the  values  of  a^/4,  and  uaia»  were  not  known  to  high 
precision.  Two  independent  studies  in  1991  established  definitive  values  for  these 
parameters.*®  For  aU  four  materials,  the  lattice  constants  are  now  known  to  within 
0.001  A,  and  the  Poisson  ratios  to  within  0.005.*®^*  We  estimate  that  the  resulting 
error  in  layer  composition,  Aii,  is  less  than  0.003. 

2.  Epilayer  peak  shift:  Another  potential  source  of  error  results  from  the  shift 
of  epilayer  peak  position  as  a  result  of  the  tail  of  the  substrate  peak.*®  This  effect 
is  illustrated  in  fig.  2.23.  Note  that  the  separation  between  the  InP  and  InGaAs 
peaks,  A9,  increases  with  increasing  InGaAs  thickness.  If  the  composition  is  obtained 
directly  from  via  Bragg’s  law,  the  results  for  the  1000  A,  500  A,  and  200  A  layers 
are  x  =  0.5995,  0.5988,  and  0.5932,  respectively.  Using  simulations,  we  calculated 
the  thickness,  as  a  function  of  mismatch,  at  which  the  error  in  composition,  Ax2,  is 
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equal  to  0.005.  For  mismatches  less  than  0.002,  the  resulting  thicknesses  were  nearly 
identical  to  the  MRP  criterion.  For  larger  mismatches,  however,  the  values  are  greater 
than  the  MRP  criterion,  as  shown  by  the  dotted  line  in  lig.  2.24.  By  comparison  to 
simulations,  one  can  always  compensate  for  the  epilayer  peak  shift;  fig.  2.24  can  be 
used  to  determine  when  such  simulations  are  necessary. 

3.  Relaxation:  We  consider  a  worst-case  scenario  in  which  a  layer  with  R  =  10% 
is  assumed  to  be  coherent.  Using  (2.6),  (2.7),  and  Vegard’s  law,  we  obtain  A®3,  the 
error  in  composition: 


Ala  = 


O'InP 


j(T),  =  ‘(T), 


20(ffl/njlj  ^GaAt(AlAs)')  \  d  J r  \  d  ) f 

where  k=0.72  for  InxGai_j,As/InP  and  0.74  for  InyAli_yAs/InP.  For  example,  a  layer 
of  Ino.6ooGao.4ooAs  on  InP  has  (Aa/a),.  =  0.0047  and  Axa  =  0.003;  the  apparent 
composition  is  Ino.597Gao.403As,  a  tolerable  error  in  most  applications. 

4.  Experimental  errors  in  A^  will  translate  into  errors  in  the  composition. 
Such  errors  are  a  function  of  several  factors  including  layer  thickness  (peak  width), 
mismatch,  background  countrate,  count  time  per  point,  and  the  angular  step  size 
between  points.  We  estimated  the  error,  A^^,  for  the  samples  shown  in  fig.  2.24.  For 
all  the  samples  with  mismatch  less  than  10“^,  A^^  <  90  arc-sec,  which  translates  to 
A®4  <  0.005.  Errors  are  generally  much  smaller  for  layers  thicker  than  500  A.  For  the 
100  and  200  A  layers  (mismatch  of  0.011),  AO^  <  180  arc-sec,  and  A®4  <  0.010.  We 
also  note  that  poor  crystalline  quality  can  broaden  epilayer  peaks,  possibly  increasing 
Ad^.  By  remaining  below  i/ic, mb  =  5,  we  will  generally  avoid  substantial  broadening 


due  to  lattice  mismatch.  Other  factors  such  as  improper  growth  temperature  can  also 
degrade  layer  quality  and  broaden  x-ray  peaks. 

5.  Differential  Form  of  Bragg’s  Law:  An  additional  error  can  result  from 
using  (2.4)  when  A9  is  not  “small.”  The  exact  expression  is:®® 


(t) 


sin(^B) 


1  = - — - 1 

1  sin{0B  +  A6) 


(2.27) 


The  error  introduced  by  using  (2.4),  Axs,  is  only  0.001  when  (Aa/o),  =  ±  0.0050, 
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but  increases  to  0.010  for  (^J^afa)r  —  i  0.017.  This  error  can  be  eliminated  in  all 
cases  by  using  (2.27). 

6.  Vegard’s  law:  The  final  potential  source  of  error  is  deviation  from  Veg- 
ard’s  law.  Based  upon  experimental  studies  of  InGaAs®^"®^  as  well  as  theoretical 
considerations,®^  we  believe  that  Vegard’s  law  is  a  good  approximation  for  the  lattice 
constant  of  InGaAs  and  InAlAs.  We  do  not,  however,  have  sufficient  information  to 
make  a  good  estimate  of  the  probable  error  in  composition  measurements. 

To  obtain  an  estimate  of  the  total  error  in  composition,  A®r,  we  use: 

A®r  =  y/Axl  +  Axl  +  Axl  +  Axl  (2.28) 

We  note  that  for  thinner  samples,  the  A0  and  layer  peak  shift  errors  may  be  near 
their  maxima,  but  the  relaxation  error  will  be  negligible.  For  thicker  samples,  the 
opposite  is  generally  true.  Using  (2.28),  we  find  that  Axt  <  0.01  for  the  shaded 
region  A  in  fig.  2.23.  In  many  cases,  Axt  wiU  be  much  smaller.  Our  data  is  limited 
for  mismatches  larger  than  0.01,  but  we  estimate  A®r  <  0,03  for  region  B. 

Figure  2.24  can  be  used  to  determine  the  optimal  thickness  for  calibration  layers. 
For  example,  if  the  mismatch  for  a  calibration  layer  is  expected  to  be  between  0.001 
and  0.004,  a  thickness  of  1000-2000  A  should  be  chosen  to  minimize  Axj. 

We  note  that  the  MRP  criterion  was  developed  for  single  layers  on  substrates, 
and  hence  is  not  strictly  valid  for  multi-layer  structures.  Based  upon  simulated  and 
experimental  results  on  multi-layer  structures  including  devices,  however,  we  have 
found  the  MRP  criterion  is  still  approximately  valid  in  most  situations.  One  could 
construct  heterostructures  for  which  this  would  not  be  the  case.  For  example,  a 
200  A  layer  of  Ino.70Gao.30As  would  not  give  a  distinct  peak  if  a  2000  A  layer  of 
Ino.6BGao.3i As  were  also  present.  Hence,  fig.  2.24  may  be  used  (with  caution)  to 
determine  whether  HRXRD  can  provide  accurate  measurements  of  layer  composition 
in  complex  heterostructures  including  devices.  In  some  cases,  complete  dynamical 
simulations  will  be  necessary. 
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2.5  Summary 


In  summary,  we  used  HRXRD  to  investigate  the  impact  of  lattice  mismatch  on  epi¬ 
taxial  layers  of  InGaAs  and  InAlAs  grown  by  MBE  on  InP.  Measurement  of  (004) 
layer  peak  width  and  interference  fringes  is  shown  to  be  a  sensitive,  non-destructive 
technique  to  assess  structural  imperfection.  We  found  that  the  crystalline  quality 
of  InGaAs  and  InAlAs  epitaxial  layers  consistently  remains  largely  unperturbed  to 
thicknesses  up  to  3-9  times  the  Matthews-Blakeslee  critical  layer  thickness  for  lattice 
mismatch  up  to  ±  1%. 

When  layers  begin  to  relax,  misfit  dislocations  form  in  an  asymmetric  pattern,  with 
the  majority  of  dislocations  parallel  to  the  [lIO]  direction  for  InGaAs  and  InAlAs  in 
compression  or  tension.  As  a  result,  the  crystal  symmetry  changes  from  tetragonal 
to  orthorhombic.  This  orthorhombic  distortion  must  be  taken  into  account  when 
measuring  the  composition  of  partially  relaxed  epilayers. 

Based  upon  our  measurements  of  both  coherent  and  relaxed  layers  as  well  as 
HRXRD  simulations,  we  have  determined  the  range  of  epilayer  thicknesses  over  which 
a  single  HRXRD  scan  yields  the  composition  of  Ina,Gai_a,As  and  InyAli_yAs  layers 
to  within  1%.  When  possible,  calibration  layers  should  be  grown  within  this  range  to 
allow  fast  and  accurate  characterization. 
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Chapter  3 


Optical  Anisotropy 


In  chapter  2,  we  used  HRXRD  to  determine  layer  composition,  relaxation,  and  crys¬ 
talline  quality.  We  did  not,  however,  obtain  any  information  on  the  mechanism  of 
structural  degradation.  In  this  chapter,  we  use  a  new  technique  to  measure  the  opti¬ 
cal  anisotropy  (OA)  of  epilayers.  We  demonstrate  that  OA  is  related  to  the  growth 
mode,  with  surface  roughness  (possibly  associated  with  3-D  growth)  causing  large 
anisotropies.  Based  upon  a  comparison  of  OA  values  and  HRXRD  peak  widths,  we 
conclude  that  a  large  optical  anisotropy  is  associated  with  a  degradation  in  layer 
quality.  Hence,  measurement  of  OA  is  a  potentially  powerful  technique  for  probing 
the  limits  of  strained  layer  epitaxy. 

In  section  3.1,  we  briefly  review  optical  properties  of  sobds  and  their  measurement 
by  conventional  ellipsometry.  We  developed  a  new  optical  characterization  technique 
known  as  variable  azimuthal- angle  ellipsometry  (VAAE).  We  describe  the  technique 
and  our  results  in  Section  3.2.  We  also  measured  selected  samples  by  reflectance 
difference  spectroscopy  (RDS),  a  characterization  tool  which  is  extremely  sensitive  to 
optical  anisotropy.  Section  3.3  includes  the  RDS  results  and  their  correlation  with 
VAAE.  We  consider  the  origin  of  the  optical  anisotropy  in  section  3.4  and  summarize 
our  findings  in  section  3.5. 
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3.1  Background 


The  application  of  an  electric  field,  E,  to  a  material  results  in  an  electric  displacement, 
D.  In  scalar  form,  the  relationship  between  these  two  quantities  is: 


I>  =  (ei  +  te2)J5  (3.1) 

where  (ei  + 162)  is  the  dielectric  constant  of  the  material.  The  dielectric  constant  is 
related  to  the  complex  refractive  index,  iV,  by: 

—  {n-\-  ikf  =  Cl  + 162  (3.2) 


where  n  is  the  real  refractive  index  and  k  is  the  extinction  coefficient  of  the  material. 
Another  useful  quantity  is  the  absorption  coefficient,  a,  which  can  be  expressed  as: 


47rfc 


a  = 


(3.3) 


where  A  is  the  wavelength  of  the  hght.  The  reciprocal  of  a.  is  the  penetration  depth, 
6]  the  intensity  of  light  drops  by  a  factor  of  1/e  for  each  8. 

EUipsometry  is  used  to  determine  optical  constants  of  materials  by  measuring  the 
state  of  polarization  of  light  after  reflection.®®’®®  The  transverse  electric  (TE)  and 
transverse  magnetic  (TM)  components  of  light  behave  differently  upon  reflection. 
In  eflipsometry,  the  ratio  of  these  two  components  is  measured.  The  fundamental 
equation  is: 


—  =  tan(V’)e‘^  (3-4) 

rs 

where  Vp  is  the  reflection  coefficient  for  p-polarized  (TM)  waves,  and  rs  is  the  reflection 
coefficient  for  s-polarized  (TE)  waves.  The  experimentally  measured  quantities,  -ip  and 
A,  are  a  function  of  the  optical  constants  and  layer  thicknesses  of  a  heterostructure, 
via  the  Fresnel  reflection  coefficients.  Hence,  given  measured  values  of  tp  and  A  one 
can  estimate  unknown  thicknesses  and/or  optical  constants. 

EUipsometry  is  commonly  used  to  measure  the  thickness  and  refractive  index  of 
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dielectric  films  on  semiconductors.  Applications  of  ellipsometry  to  III-V  epilayers, 
however,  have  been  hmited.  This  is  partly  due  to  the  fact  that  the  epilayers  are  usu¬ 
ally  opaque  at  the  wavelength  of  conventional  ellipsometers,  making  it  impossible  to 
measure  layer  thickness.  Furthermore,  the  presence  of  surface  oxides  makes  it  difficult 
to  accurately  measure  the  optical  constants  of  the  epilayers.®^  Recently,  spectroscopic 
ellipsometry  has  been  applied  to  measure  layer  thicknesses  in  multi-layer  heterostruc¬ 
tures  such  as  field-effect  transistors.®® 


3.2  Variable  Azimuthal-Angle  Ellipsometry 

We  characterized  our  InGaAs  and  InAlAs  epitaxial  layers  with  a  Gaertner  model 
L116B  automated  ellipsometer.  Both  circularly  and  Hnearly  polarized  incident  light 
from  a  He-Ne  laser  (A  =  6328A)  were  used  with  the  angle  of  incidence  fixed  at  70" 
from  vertical.  For  Ino.63Gao.47As,  the  absorption  coefficient  at  6328  A  is  9  x  10^ 
cm~^.®®  The  corresponding  penetration  depth  is  1100  A.  We  estimate  a  =  4  x  10^ 
cm  ^  and  S  =  2500  A  for  Ino.B2Alo.48As. 

The  ellipsometer  is  equipped  with  a  rotating  stage;  measurements  were  made  as 
a  function  of  the  azimuthal  angle  a  as  shown  in  fig.  3.1.  We  use  the  same  convention 
as  in  chapter  2:  a  =  0"  when  the  projection  of  the  incident  laser  (x-ray)  beam  on 
the  wafer  surface  is  in  the  [110]  direction.  We  refer  to  this  technique  as  variable 
azimuthal- angle  ellipsometry  (VAAE).®®’^®® 

We  typically  observe  a  systematic  variation  of  A  as  a  function  of  azimuthal  angle. 
As  an  example,  we  plot  A  versus  a  for  sample  1530  in  fig.  3.2.  This  sample  consists 
of  a  single  1000  A  layer  of  Ino.38Gao.62As  on  InP.  We  note  the  two-fold  symmetry  with 
maxima  at  a  =  90"  and  270",  corresponding  to  incident  light  parallel  to  the  [110]  and 
[110]  directions.  The  data  can  be  fit  with  a  cosine  function  of  the  form: 

A  =  -f-  B^co8[2(a  —  C^)]  (3'5) 

where  Aa,  Ba,  and  Ca  are  fitting  parameters.  The  best-fit  to  (3.5)  was  calculated 
by  least  squares  and  is  shown  as  a  solid  line  in  fig.  3.2.  The  parameters  are:  A  a  = 
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Figure  3-1:  Geometry  of  ellipsometry  measurements;  a  is  the  azimuthal  angle  (angle 
of  rotation)  of  the  projection  of  the  beam  on  the  sample  with  respect  to  the  [110] 
direction.  The  angle  of  incidence  is  fixed  at  20®  &om  horizontal. 

145.1®,  =  8.21®,  and  =  90.7®.  For  this  layer,  =  9-  As  a  result  of  the 

large  mismatch,  the  crystalline  quality  is  poor,  as  illustrated  by  an  HRXRD  FWHM 
ratio  of  12. 

For  layers  with  small  mismatch,  we  still  observe  cosine-shape  variations  of  A  with 
a,  but  the  amplitude  of  the  cosine  function,  is  always  small.  As  an  example, 
we  show  the  data  for  sample  1877  in  fig.  3.3.  For  this  InGaAs  layer,  l-l) 

the  HRXRD  scan  (fig.  2.5)  indicates  high  crystalline  quality,  and  Ba  =  0.31®.  (The 
InGaAs  layer  in  1877  is  too  thick  for  a  significant  portion  of  the  light  to  reach  the 
InAlAs.) 

For  a  few  samples  with  very  large  azimuthal  variations  in  A,  we  also  observe  a 
small  but  systematic  variations  in  with  a  cosine  shape  and  maxima  at  a  —  90 
and  270®.  For  most  samples,  however,  the  variation  in  ip  with  a  appears  random.  We 
conclude  that  A  is  more  sensitive  to  optical  anisotropy  than  ijf.  Hence,  we  take  Ba 
as  a  measure  of  optical  anisotropy. 

We  also  observe  variations  in  A  a  and  A^,  the  average  values  of  A  and  ip,  from 
sample  to  sample.  We  have  found,  however,  that  Aa  and  A^  vary  substantially  with 
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[110]  [110]  [110]  [110]  [110] 


Figure  3-2:  EUipsometric  parameter  A  as  a  function  of  azimuthal  angle  for  sample 
1530,  mismatched  InGaAs  in  tension.  Solid  line  is  a  least-squares  fit  to  eq.  (3.5). 


[110]  [110]  [110]  [110]  [110] 


Figure  3-3:  EUipsometric  parameter  A  as  a  function  of  azimuthal  angle  for  sample 
1877,  nearly  lattice-matched  InGaAs.  Solid  line  is  a  least-squares  fit  to  eq.  (3.5). 
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Figure  3-4:  Amplitude  of  the  cosine  function  fits  for  the  ellipsometric  parameter  A 
as  a  function  of  lattice  mismatch  for  a  set  of  InGaAs  layers  in  tension. 

position  on  a  single  sample.  In  addition,  sample  cleaning  changes  Aa  and  A^.  In 
contrast,  is  relatively  insensitive  to  sample  cleaning  and  position. 

We  now  examine  for  sets  of  InGaAs  layers  in  tension  and  compression.  In  fig. 
3.4,  we  plot  jBa  as  a  function  of  mismatch  for  samples  1876-1880,  3000  A  layers  of 
InGaAs  in  tension.  The  HRXRD  rocking  curves  for  this  set  are  shown  in  fig.  2.5. 
Within  the  experimental  error  bars,  the  OA  is  constant  {B^  =  0.2  —  0.4°)  for  1876- 
1879;  these  samples  have  HRXRD  FWHM  ratios  close  to  unity  and  Pendellosung 
fringes.  In  contrast,  sample  1880  exhibits  degraded  crystalline  quality  and  a  larger 
OA,  with  B^  =  2.8°. 

The  OA  follows  a  similar  pattern  for  samples  1630-1634,  1000  A  layers  of  InGaAs 
in  compression.  is  plotted  as  a  function  of  mismatch  in  fig.  3.5;  the  HRXRD 
scans  were  given  in  fig.  2.6.  The  OA  is  small  for  the  samples  with  narrow  x-ray 
peaks  (1630,  1632,  and  1633)  and  large  for  the  sample  (1634)  with  a  broadened  x-ray 
rocking  curve. 

Optical  anisotropy  was  also  observed  for  InAlAs  layers.*®  As  an  example  of  a 
structure  with  large  anisotropy,  in  fig.  3.6  we  plot  A  versus  a  for  sample  1843, 


72 


3 


CO 

liJ 

lU 

oc 

o 

HI 

o 


00 


X  1-x 

1634  . 

I 

- 

t  ■ 

#1630  1632  1633 

■5  5 . J 

o.doo  0.005 

MISMATCH 

0.010 

Figure  3-5:  Amplitude  of  the  cosine  function  fits  for  the  elhpsometric  parameter  A 
as  a  function  of  lattice  mismatch  for  a  set  of  InGaAs  layers  in  compression. 

a  3000  A  layer  of  Ino.esAlo.asAs.  As  with  InGaAs,  the  data  approximately  follow 
a  cosine  law  with  maxima  at  a  =  90°  and  270°.  The  OA  of  InAlAs  samples  also 
correlates  with  HRXRD  results.  For  1843,  the  HRXRD  FWHM  ratio  is  23,  Ba  = 
7.1°,  and  =  20.  In  fig.  3.7,  we  plot  Ba  versus  lattice  mismatch  for  samples 

3098-3102,  1000  A  layers  of  InAlAs  in  compression  (see  fig.  2.8  for  HRXRD  scans). 
The  anisotropy  is  largest  for  samples  3100,  3101,  and  3102.  These  three  samples  have 
inferior  HRXRD  rocking  curves. 

As  shown  in  the  figures  3. 2-3. 7,  we  have  observed  large  optical  anisotropies  for 
InGaAs  in  tension,  InGaAs  in  compression,  and  InAlAs  in  compression.  In  each 
case,  the  orientation  dependence  is  the  same.  We  usually  do  not,  however,  observe 
large  anisotropies  for  InAlAs  in  tension  (AlAs-rich).  For  the  1701-1705  sample  set, 
HRXRD  scans  (fig.  2.7)  revealed  structural  degradation  for  two  samples,  but  Ba  was 
about  0.3°  for  aU  five  samples.  In  chapter  6,  we  will  discuss  possible  reasons  for  the 
differences  in  InAlAs  layers  in  tension  and  compression. 

We  have  demonstrated  that  OA  is  small  for  small  lattice  mismatch  and  usually 
becomes  large  when  mismatch  increases  and  crystalline  quality  deteriorates.  For  sev- 
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Figure  3-6:  Ellipsometric  parameter  A  as  a  function  of  azimuthal  angle  for  sample 
1843,  InAlAs  in  compression.  Solid  line  is  a  least-squares  fit  to  eq.  (3.5). 
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Figure  3-7:  Amplitude  of  the  cosine  function  fit  for  the  ellipsometric  parameter  A  as 
a  function  of  lattice  mismatch  for  a  set  of  InAlAs  layers  in  compression. 
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Figure  3-8:  Amplitude  of  the  cosine  function  lit  for  the  ellipsometric  parameter  A  as 
a  function  of  thickness  for  a  set  of  InAlAs  layers  in  compression. 

eral  sets  of  samples,  we  have  also  observed  that  the  OA  decreases  when  the  mismatch 
or  thickness  becomes  too  large.  An  example  is  shown  for  samples  1842-1846  in  fig- 
3.8.  Each  sample  contains  a  single  layer  of  Ino.esAlo.sBAs.  The  OA  initially  increases 
with  layer  thickness  but  then  drops  to  a  small  value  =  0.7°)  at  10,000  A.  We 
have  also  observed  this  behavior  for  InGaAs  in  tension  and  compression. 

To  confirm  the  apparent  connection  between  OA  and  crystalline  quality,  we  plot 
versus  the  HRXRD  FWHM  ratio  for  our  entire  set  of  samples  in  fig.  3.9.^  For 
InGaAs  (fig.  3.9a)  it  is  obvious  that  much  larger  values  of  optical  anisotropy  occur 
in  tension  than  in  compression,  but  several  of  the  layers  in  compression  have  B^’s 
which  are  clearly  larger  than  the  values  for  nearly  lattice-matched  samples.  The  OA 
is  always  small  (Ba  <  0.8°)  when  the  FWHM  ratio  is  less  than  1.5.  In  chapter 
2,  we  used  a  similar  FWHM  ratio,  1.8,  as  a  cutoff  for  “high”  crystalline  quality.^ 
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^We  excluded  multilayer  samples  except  for  a  few  two-layer  structures  with  a  neatly  lattice- 
matched  buffer  layer  such  as  samples  1876-1880. 

^Three  of  the  four  layers  in  fig.  3.9a  with  5a  >  1.0°  and  FWHM  ratios  between  1.5  and  1.8 
were  thin  (t  <  lOOOA).  For  such  thin  layers,  we  expect  HRXRD  to  be  less  sensitive  to  crystalline 
imperfection  than  for  thicker  layers  because  the  peaks  are  intrinsically  broad.  Hence,  for  thin  layers 
VAAE  may  detect  degradation  in  crystalline  quality  before  HRXRD. 
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For  InGaAs,  we  conclude  that  large  OA  correlates  with  a  degradation  in  crystalhne 
quality.  Figure  3.9a  includes  many  samples  with  a  large  FWHM  ratio  and  small 
Hence,  a  small  value  of  OA  does  not  always  imply  good  crystalline  quality. 

For  InAlAs  (fig.  3.9b),  we  observe  several  compressively-strained  samples  with 
large  optical  anisotropies,  but  only  one  sample  in  tension  has  a  Ba  greater  than 
1“.  We  again  observe  that  the  OA  is  always  small  for  samples  with  high  crystalline 
quality.  We  also  note  that  for  both  InAlAs  and  InGaAs,  Ba  »s  always  positive,  even 
for  the  samples  with  small  OA.  If  the  true  OA  were  zero  for  nearly  lattice-matched 
samples,  we  would  expect  experimental  noise  to  result  in  positive  and  negative  values 
of  Ba-  We  conclude  that  even  the  small  optical  anisotropies  are  real. 

Morris  ei  a/.'®'  recently  applied  the  VAAE  technique  to  InGaAs  grown  by  MOCVD. 
They  measured  significant  anisotropies  for  3.0  fim  layers  of  Ino.i6Gao.84As  on  GaAs. 
Their  observed  maximum  values  of  A  were  apparently  in  the  same  direction  as  ours. 
This  result  demonstrates  that  optical  anisotropies  are  not  confined  to  MBE-grown 

layers. 

3.3  Reflectance  Difference  Spectroscopy 

Reflectance  difference  spectroscopy  (RDS)  was  invented  in  1985.^°^'^°®  In  this  tech¬ 
nique,  hght  is  reflected  from  a  sample  at  near-normal  incidence  and  the  reflectance 
is  measured  in  2  different  directions.  The  incident  light  is  linearly  polarized.  In  the 
case  of  an  (001)  substrate,  the  hght  is  polarized  at  45"  to  the  [110],  hence  illuminating 
the  [110]  and  [iTO]  directions  equally,  as  shown  in  fig.  3.10.  The  reflected  hght  is 

measured  in  the  [110]  and  [iTO]  directions.  We  define: 

mo  =  complex  reflectance  paraUel  to  the  [110]  axis  of  substrate 
rjio  =  complex  reflectance  paraUel  to  the  [iTO]  axis  of  substrate 
If  the  sample  is  isotropic,  mo  =  and  the  reflected  hght  simply  reconstructs  the 
hnear  polarization  of  the  incident  hght.  Any  anisotropy  in  the  sample  will  result  in 
riio  7^  iiYo-  RDS  can  measure  reflectance  differences  as  small  as  5  xl0~®.^°^ 

A  major  difference  between  RDS  and  elhpsometry  is  that  elhpsometric  measure- 
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Figure  3-9:  Amplitude  of  the  cosine  function  fit  for  the  ellipsometric  parameter  A  as 
a  function  of  the  ratio  of  the  HRXRD  experimental  and  theoretical  (004)  peak  widths 
for  (a)  InGaAs  and  (b)  InAlAs  layers. 
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Figure  3-10:  Polarization  of  incident  light  with  respect  to  the  crystallographic  direc¬ 
tions  of  the  sample  for  reflectance  difference  spectroscopy  on  an  (001)  substrate. 

ments  are  made  at  an  oblique  angle  of  incidence.  Hence,  even  isotropic  samples  cause 
a  change  in  polarization.  It  is  this  change  in  polarization  that  is  measured  and  used 
to  estimate  layer  thickness  and  refractive  index  in  ellipsometry.  K  the  sample  is  op¬ 
tically  anisotropic,  the  analysis  becomes  very  difficult.  With  RDS,  however,  the  only 
factor  that  changes  the  polarization  is  optical  anisotropy. 

We  investigated  a  variety  of  matched  and  mismatched  InGaAs/InP  and  InAlAs/InP 
heterostructures,  measuring  reflectance  differences  over  the  energy  range  of  1.5  to  5.5 
eV.  Typical  results  are  shown  in  figs.  3.11  and  3.12  where  we  plot  both  the  magnitude, 
|Ar/rl,  and  phase,  A/3,  of: 

^  ~  (3.6) 

^  ^lio  +  ’’110 

for  samples  1879  and  1880  (see  figs.  2.5  and  3.4  for  HRXRD  and  VAAE  data). 
The  vertical  scales  differ  by  a  factor  of  ten  on  the  two  plots.  The  magnitudes  of 
the  optical  anisotropies  for  1879  are  less  than  0.004  over  the  entire  spectral  range.® 

^|Ar/rl  and  A)9  can  be  determined  from  a  single  RDS  measurement.  The  resolution  u  improved 
if  |Ar/r|  and  A/9  are  determined  by  subtracting  two  measurements  which  are  taken  90®  apart.  We 
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Figure  3-11:  Magnitude  and  phase  of  RDS  signal  as  a  function  of  photon  energy  for 
1879,  a  sample  with  good  crystalline  quality. 

(|Ar/r|  is  dimensionless;  has  units  of  radians.)  In  contrast,  1880  exhibits  much 
larger  anisotropies,  with  [At'/t’I  and  A/?  reaching  0.062  and  -0.027,  respectively.  The 
relative  magnitudes  are  in  qualitative  agreement  with  our  ellipsometry  measurements 
which  yielded  Ba  values  of  0.4"  and  2.8"  for  1879  and  1880,  respectively. 

For  all  25  samples  measured  by  RDS,  we  observed  a  strong  correlation  with  the 
optical  anisotropy  measured  by  ellipsometry.  Our  ellipsometry  measurements  were 
made  at  a  photon  energy  of  1.96  eV.  For  comparison,  we  plot  the  values  of  |Ar/r| 
and  A^  at  1.96  eV  versus  Ba  for  both  InGaAs/InP  and  InAlAs/InP  structures  in 
fig.  3.13.  There  is  clearly  a  correlation  between  |  Ar/rj  and  Ba,  and  an  even  stronger 
correlation  between  A/S  and  Ba- 

3.4  Physical  Origin  of  Optical  Anisotropy 

We  first  consider  the  cause  of  the  relatively  small  optical  anisotropies  observed  for  lay¬ 
ers  of  high  crystalline  quality.  The  RDS  signal  for  samples  such  as  1879  is  comparable 

used  this  latter  procedure;  it  accounts  for  the  smalls  negative  values  of  |Ar/r|  in  fig.  3.11. 
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Figure  3-12:  Magnitude  and  phase  of  RDS  signal  as  a  function  of  photon  energy  for 
1880,  a  sample  with  degraded  crystalline  quahty. 

in  magnitude  to  RDS  signals  from  bulk  GaAs.  Aspnes  attributed  the  GaAs  results 
to  surface  chemistry:  for  example,  the  absorption  associated  with  Ga-Ga  surface 
dimer  bonds. Similar  anisotropies  are  expected  in  nearly  lattice-matched  epilayers. 
Hence,  for  layers  of  high  crystalhne  quality,  we  attribute  the  relatively  small  optical 
anisotropy  (measured  by  RDS  and  VAAE)  to  surface  chemistry. 

We  have  observed  large  optical  anisotropies  in  layers  ranging  from  a  few  hundred 
angstroms  to  one  micron.  Since  the  penetration  depth  in  our  epilayers  is  1000  to 
3000  A,  virtually  none  of  the  hght  can  penetrate  through  a  one  micron  layer  to  the 
layer /substrate  interface.  We  conclude  that  the  origin  of  the  large  OA  is  in  the  top 
few  thousand  angstroms  of  epilayer,  not  the  substrate  or  the  interface. 

One  possible  explanation  for  the  large  OA  in  mismatched  layers  is  the  piezo-optical 
effect,  a  strain-induced  change  in  refractive  index.^®’^°®  If  piezo-optical  effects  were 
dominant,  we  would  expect  the  largest  anisotropies  in  samples  with  a  large  difference 
in  strain  in  the  [110]  and  [110]  directions.  The  strain  in  the  [110]  direction  is  equal 
to  (Ao/a)r  —  (Ao/a)||[iio])  and  similarly  for  the  [lIO]  direction.  Hence,  the  difference 
in  strain  is  (Ao/a)||[iio]  —  (Aa/o)||[iio]).  In  fig.  3.14,  we  plot  as  a  function  of  the 
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Figure  3-13:  Comparison  of  a)  magnitude,  and  b)  phase,  from  reflectance  difference 
spectroscopy  with  anisotropy  parameter  from  variable  azimuthal  angle  eUipsometry. 
The  photon  energy  is  1.96  eV  for  both  techniques. 
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Figure  3-14:  Optical  anisotropy  as  a  function  of  the  difference  in  strain  in  the  [110] 
and  [liO]  directions  for  partially  relaxed  InGaAs  layers. 

strain  difference  for  partially  relaxed  InGaAs  layers.  With  one  exception,  the  samples 
do  not  have  large  OA,  despite  the  large  asymmetries  in  strain.  Many  samples  (see 
fig.  3.9)  exhibit  large  OA,  but  are  essentially  unrelaxed  and  hence  do  not  have 
strain  asymmetries.  Furthermore,  piezo-opticjil  effects  should  produce  a  reversal  in 
the  direction  of  the  OA  when  changing  from  tension  to  compression.  Our  epilayers 
exhibited  the  same  pattern  for  both  tension  and  compression.  In  addition,  piezo- 
optical  effects  are  expected  to  exhibit  large  changes  near  the  critical  points  in  energy, 
but  we  observe  no  such  features  in  our  RDS  data.  Hence,  we  conclude  that  piezo- 
optical  effects  are  not  the  dominant  mechanism  for  optical  anisotropy. 

As  mentioned  in  chapter  1,  composition  modulation  may  occur  in  InGaAs  and  In- 
AlAs  layers.  We  considered  it  as  a  potential  cause  of  optical  anisotropy.  We  obtained 
plan- view  transmission  electron  microscopy  (TEM)  images  of  eight  samples.^  One  of 
the  InAlAs  samples,  3099,  exhibited  contrast  modulation  as  shown  in  fig.  3.15.  The 
contrast  modulation  has  a  period  of  365  A.  We  beheve  it  is  associated  with  variations 
in  composition.  X-ray  measurements  on  this  sample  (fig.  2.8)  reveal  good  crystaUine 

^TEM  by  F.  Peiro  and  Prof.  A.  Cornet  at  the  Univ.  of  Barcelona 
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Figure  3-15;  Plan-view  transmission  electron  micrograph  for  sample  3099,  1000  A  of 
Ino.579Alo.421  As.  A  bright-held,  two-beam  image  with  g=[220]  was  used.  A  contrast 
modulation  with  a  period  of  365  A  is  visible  along  the  <  100  >  directions. 

quality.  The  OA  (fig.  3.7)  was  small.  The  lack  of  OA  in  a  single  sample  with  com¬ 
position  modulation  is  not  definitive.  We  note,  however,  that  the  modulation  is  in 

a  <  100  >  direction,  This  is  in  agreement  with  other  reports.®^’ The  <  100  > 

directions  correspond  to  a  =  45,  135,  225,  and  315"  on  figures  such  as  3.2.  Hence, 

composition  modulation  could  only  produce  optical  anisotropies  in  these  directions, 
not  along  the  [110]  and  [lIO]  as  observed. 

Sun  et  al.  reported  large  asymmetries  in  the  electron  mobility  for  MOCVD-grown 
InGaAs/GaAs  heterostructures. They  suggest  that  the  formation  of  islands  is  re¬ 
sponsible.  Islands  are  known  to  extend  faster  in  the  [110]  direction  than  in  the  [110] 
direction  in  MOCVD.^°®  Such  a  mechanism  might  also  explain  our  results.  In  MBE, 
the  preferential  growth  direction  is  perpendicular  to  that  in  Thus, 

we  could  expect  islands  elongated  along  the  [lIO]  direction  when  three-dimensional 
growth  occurs.  Such  islands  were  observed  in  the  InGaAs/GaAs  system.®®  It  is  plau¬ 
sible  that  elongated  islands  would  affect  a  layer’s  optical  properties. 

Further  support  to  the  theory  of  elongated  islands  is  found  in  the  work  of  Acher 
et  They  performed  in  situ  RDS  measurements  during  the  growth  of  InAs 
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Figure  3-16:  Plan-view  transmission  electron  micrograph  for  sample  1530.  The  image 
was  taken  under  bright-field,  two-beam  conditions  with  g  =  [220].  Agglomerations  of 
threading  dislocations  are  observed  (A).  Islands  aligned  along  [lIO]  are  clearly  visible 
(B).  Misfit  dislocations  appear  to  have  formed  (C)  but  are  partially  masked  by  other 
features. 

on  InP  and  InP  on  GaAs  by  MOCVD.  Optical  anisotropies  up  to  5%  for  InAs/InP 
and  over  100%  for  InP/GaAs  were  observed.  They  attempted  to  fit  their  data  with 
models  based  upon  roughness  anisotropy  at  both  the  epilayer /substrate  interface  and 
the  epilayer  surface.  The  models  (with  fitting  parameters)  gave  reasonable  agreement 
to  the  observed  optical  anisotropies  as  a  function  of  growth  time. 

To  confirm  the  existence  of  elongated  islands  or  roughness  anisotropy,  we  obtained 
TEM  micrographs  of  sample  1530  which  has  a  large  OA  (fig.  3.2).  In  fig.  3.16,  we 
show  the  plan- view  image.  The  TEM  micrograph  reveals  three-dimensional  growth 
with  islands  along  the  [110]  direction.  Cross-sectional  TEM  (fig.  3.17)  shows  that  the 
film  is  continuous,  but  has  a  very  irregular  morphology,  with  layer  thickness  ranging 
from  830  to  1400  A.  Stacking  faults  are  also  observed. 

We  also  obtained  plan-view  TEM  micrographs  of  samples  1879  and  1880,  figs. 
3.18  and  3.19.  For  1879,  the  surface  morphology  was  smooth,  indicative  of  two- 
dimensional  growth.  Only  a  few  misfit  dislocations  were  found.  The  intersection 
of  two  dislocations  is  shown  in  the  figure.  These  dislocations  were  confirmed  to  be 
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Figure  3-17:  Cross-section  transmission  electron  micrograph  for  sample  1530.  The 
image  was  taken  in  the  [110]  direction  under  a  g  =  [220]  bright-held  condition.  Note 
the  variations  in  epilayer  thickness  and  stacking  faults  (SF). 

of  the  60"-type.^^^  This  sample  exhibited  excellent  crystalline  quality  (fig.  2.5)  and 
small  OA  (figs.  3.4  and  3.11).  In  contrast,  TEM  reveals  a  rough  surface  morphology 
with  hillocks  for  sample  1880.  Dislocations  are  also  visible,  but  in  the  [lIO]  direction 
only,  with  an  average  spacing  of  about  1.5  fim.  The  sample  has  a  large  OA  (figs. 
3.4  and  3.12)  and  degraded  structural  quality  (fig.  2.5).  Elongated  islands  are  not 
clearly  visible  on  sample  1880  as  they  were  on  sample  1530.  The  OA  for  sample  1880, 
although  large  compared  to  lattice- matched  samples,  is  much  smaller  than  sample 
1530. 

As  mentioned  earlier,  we  often  observe  that  the  OA  of  mismatched  layers  returns 
to  a  small  value  when  the  layer  thickness  becomes  too  large.  This  usually  occurs  when 
the  layer  begins  to  relax  substantially  {R  >  10%).  A  TEM  image  of  such  a  sample 
(not  shown)  revealed  poor  crystalline  quality,  with  >  10*”  cm“*  threading  dislocations 
and  stacking  faults.  We  do  not  believe  that  the  layer  has  returned  to  smooth,  two- 
dimensional  growth.  Instead,  the  growth  is  thought  to  be  rough  or  three-dimensional, 
but  relatively  isotropic.  We  note  that  Acher  et  al.  also  observed  the  OA  returning  to 
a  small  value  when  the  thickness  of  InAs  on  InP  reached  3000  A.**^ 
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Figure  3-18:  Plan-view  transmission  electron  micrograph  for  sample  1879  with  g  = 
[220].  Two  intersecting  60°  misfit  dislocations  are  visible.  The  morphology  is  smooth, 
indicating  two-dimensional  growth. 


Figure  3-19:  Plan-view  transmission  electron  micrograph  for  sample  1880  with  g 
=  [220].  Misfit  dislocations  are  visible  along  the  [110]  direction  only.  The  surface 
morphology  is  rough,  with  hillocks  elongated  toward  [110]  and  [110]. 
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3.5  Summary 


We  developed  a  new  technique,  variable  azimuthal-angle  ellipsometry,  for  measuring 
optical  anisotropy.  We  observed  anisotropy  in  InGaAs/InP  and  InAlAs/InP  het¬ 
erostructures  using  both  VAAE  and  Reflectance  Difference  Spectroscopy.  For  nearly 
lattice-matched  epilayers  with  high  crystalline  quality  (measured  by  HRXRD),  the 
optical  anisotropy  is  always  small  but  non-zero.  We  attribute  this  small  anisotropy  to 
surface  chemistry.  For  samples  with  larger  mismatch  and  degraded  crystalline  qual¬ 
ity,  the  opticsJ  anisotropy  is  often  large.  Based  upon  a  comparison  to  TEM  images 
and  other  arguments,  we  conclude  that  large  anisotropies  result  from  strain-mduced 
surface  roughness,  possibly  associated  with  three-dimensional  growth. 
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Chapter  4 


Relaxation  Kinetics  and  Thermal 
Stability 

In  chapter  2,  we  demonstrated  that  complete  relaxation  of  epilayer  strain  does  not 
occur  when  the  equilibrium  critical  layer  thickness  is  reached.  Similar  results  in  other 
material  systems  have  been  attributed  to  metastability.  The  assumption  is  that  misfit 
dislocations  do  not  form  or  propagate  during  growth  due  to  kinetic  barriers.  This 
suggests  that  layers  beyond  tc,MB  with  high  crytahine  quality  may  relax  during  high- 
temperature  processing  steps  or  device  operation.  If  so,  they  would  not  be  suitable  for 
most  device  applications.  Our  primary  goals  in  this  chapter  are  to  assess  the  thermal 
stabihty  of  high-quahty  strained  InGaAs  and  InAlAs  layers  and  to  understand  the 

kinetics  of  relaxation  during  growth  and  anneahng. 

We  assessed  the  thermal  stability  of  our  epilayers  by  two  techniques.  In  section 
4.1,  we  describe  the  use  of  HRXRD  to  measure  structural  changes  resulting  from  high- 
temperature  anneals.  In  section  4.2,  we  study  the  effect  of  anneahng  on  the  electron 
mobihty  in  MODFETs  with  mismatched  InGaAs  channels.  We  examine  relaxation 
kinetics  in  section  4.3,  comparing  our  as-grown  relaxation  data  to  a  pubhshed  model. 
We  then  explain  the  thermal  stabihty  of  our  mismatched  layers  by  considering  the 
energetics  for  nucleation  of  dislocation  half-loops.  Our  findings  are  summarized  in 

section  4.4. 
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4.1  Annealing:  Structural  Stability 

We  assessed  the  thermal  stability  of  mismatched  InGaAs  and  InAlAs  layers  by  HRXRD 
measurements  before  and  after  rapid  thermal  annealing  (RTA).  The  anneals  were  per¬ 
formed  in  an  AG  Associates  Heatpulse  210  system  with  a  nitrogen  ambient.  Samples 
were  placed  face-down  on  a  GaAs  wafer  during  the  anneal  to  minimize  thermal  de¬ 
composition.  The  temperature  was  measured  by  a  thermocouple  attached  to  a  Si 
wafer  and  is  beUeved  to  be  accurate  to  within  10"C.  For  MBE-grown  In,Gai_,As 
and  InyAli-yAs,  we  typically  observe  variations  in  ®  or  y  of  0.005  across  a  sample. 
Such  variations  could  be  mistaken  for  relaxation  during  annealing.  To  avoid  this 
problem,  we  performed  sequential  anneals  on  small  (~5  x  5  mm)  pieces,  measuring 
HRXRD  rocking  curves  after  each  anneal. 

We  begin  by  considering  a  layer  which  relaxed  substantially  during  growth.  Sam¬ 
ple  1441  consists  of  1.0  /im  of  Ino.6iGao.39 As  on  InP.  The  (004)  HRXRD  rocking  curve 
is  shown  in  fig.  4.1.  Based  upon  this  and  (115)  scans,  we  determined  that  R  was  52% 
before  annealing.  We  annealed  the  sample  for  60  seconds  at  850  C.  As  shown  in  fig. 
4.2,  annealing  caused  the  layer  peak  to  move  closer  to  the  substrate  peak.  This  shift 
results  from  additional  relaxation,  with  (Aa/a)]]  increasing  and  (Ao/a)j.  decreasing. 
After  the  first  anneal,  R  =  67%.  We  performed  additional  2-minute  and  7-minute 
850'’C  anneals  (not  shown)  on  the  same  piece  of  1441  and  observed  further  relaxation, 
with  R  =  70%  and  75%,  respectively.  We  note  that  after  the  second  and  third  anneals 
at  850‘’C,  severe  degradation  of  the  surface  was  obvious  to  the  naked  eye,  despite  the 
use  of  a  GaAs  cap.  Hence,  it  may  not  be  feasible  to  achieve  complete  relaxation  by 
annealing  samples  such  as  1441.  Based  upon  these  results,  we  use  850  C  as  an  upper 
limit  for  annealing  temperature. 

The  InGaAs  layer  of  sample  1441  was  33  times  and,  as  expected,  exhibited 

poor  crystalline  quality  even  before  annealing.  For  most  device  applications,  layers 
of  high  crystalline  quality  are  required.  We  devote  the  rest  of  this  and  the  following 
section  to  the  thermal  stability  of  high-quality  strained  layers. 

We  are  primarily  interested  in  lattice  relaxation,  but  interdiffusion  may  also  occur 
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Figure  4-1;  HRXRD  scans  of  sample  1441,  a  partially  relaxed  layer  of  InGaAs  wRh 
tjtcMB  =  33)  before  and  after  annealing.  The  shift  in  the  layer  peak  indicates  lattice 
relaxation  during  the  anneal. 

during  high-temperature  anneals.  Using  an  interdiffusion  coefficient  of  Dm  =  10"^'* 
cmVsec  for  InGaAs/InP  at  we  calculate  =  77  A  for  a  60  second 

anneal.  This  suggests  that  interdiffusion  may  degrade  the  interface.  To  determine  the 
effect  of  interdiffusion  on  HRXRD  rocking  curves,  we  performed  a  series  of  anneals 
on  a  nearly  lattice-matched  epilayer.  Sample  1331  consists  of  a  single  4000  A  layer  of 
Ino.B37Gao.463As,  with  tItcMB  =  0.5.  This  layer  should  be  thermodynamically  stable 
with  respect  to  misfit  dislocations.  In  fig.  4.2,  we  show  the  rocking  curves  before 
anneals  and  after  60  seconds  at  650,  750,  and  850°C.  The  as-grown  curve  exhibits 
Pendellosung  fringes  and  an  epilayer  FWHM  ratio  of  1.1,  indicating  high  crystalline 
quality.  Interdiffusion  across  the  epilayer/substrate  interface  could  cause  a  loss  of 
fringes  and  a  broadening  of  the  layer  peak.  After  annealing,  there  is  a  reduction  in 
the  fringe  intensity  on  the  low-angle  side,  but  no  change  on  the  high-angle  side.  The 
FWHM  ratio  remained  between  1.1  and  1.3.  We  conclude  that  if  interdiffusion  is 
occurring,  its  effect  on  the  HRXRD  rocking  curves  is  small. 

We  now  examine  the  thermal  stability  of  high-quality  layers  beyond  the  Matthews- 
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Figure  4-2:  HRXRD  scans  of  sample  1331,  a  coherent  layer  of  InGaAs  with  — 

0.5,  before  and  after  annealing. 

Blakeslee  critical  layer  thickness.  Before  discussing  our  results  for  InGaAs  and  InAlAs, 
it  is  instructive  to  consider  a  published  result  for  the  SiGe/Si  system.  A  1700  A 
layer  of  Sio.82Geo.i8  was  grown  at  500" C  and  annealed  at  temperatures  from  550°  C  to 
850°C  for  60  seconds. The  rocking  curves  are  shown  in  fig.  4.3.  For  temperatures 
up  to  600°C,  almost  no  changes  are  observed.  After  a  650°C  anneal,  however,  the 
Pendellosung  fringes  disappear.  Higher  temperatures  cause  a  broadening  of  the  layer 
peak  and  a  decrease  in  the  epilayer-substrate  peak  separation,  indicating  substantial 
lattice  relaxation.  We  note  that  =  8  for  this  layer. 

In  fig.  4.4,  we  show  the  effects  of  anneahng  on  sample  1879  which  included  a 
3000  A  layer  of  Ino.477Gao,523As.  For  the  InGaAs  layer,  =  7.3.  The  results 

are  quite  different  than  for  the  SiGe  layer,  although  the  growth  temperature,  layer 
thickness,  and  mismatch  were  similar.  The  only  effect  of  annealing  is  a  slight  loss  in 
fringe  intensity,  similar  to  1331,  the  nearly  lattice-matched  “control.”  We  conclude 
that  the  strained  InGaAs  layer  is  not  relaxing  during  the  anneals.  As  shown  in  fig. 
4.5,  InAlAs  behaves  in  a  similar  way.  In  this  case,  the  1800  A  layer  of  Ino.44Alo.B6As  is 
8.6  times  We  will  discuss  the  reason  for  the  different  behavior  of  InGa(Al)As 
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Figure  4-3:  HRXRD  scans  of  a  1700  A  layer  of  Sio.82Geo.i8  on  Si  with  =  8- 

The  layer  was  grown  at  500®C  and  annealed  for  60  seconds  at  550-850"C. 

and  SiGe  in  chapter  6. 

We  also  used  HRXRD  to  examine  the  structural  stability  of  a  device  heterostruc¬ 
ture.  The  cross-section  of  the  modulation-doped  field-effect  transistor  (MODFET)  is 
shown  in  fig.  4.6.  All  the  layers  are  nominally  lattice-matched  to  the  InP  substrate 
except  the  500  A  channel.  We  performed  anneals  on  MODFET  4089  which  has  an 
Ino.MGao.3sAs  channel,  with  =  2.6.  HRXRD  scans  for  as-grown  and  60”  @ 

800"C  are  shown  in  fig.  4.7.  In  this  case,  different  pieces  of  the  wafer  were  used  for  the 
two  measurements.  Hence,  the  slight  difference  in  the  position  of  the  channel  peak 
is  probably  a  result  of  wafer  nonuniformity.  The  important  finding  is  that  fringes 
remain  after  annealing,  demonstrating  the  thermal  stability  of  this  pseudomorphic 
device  structure. 


4.2  Annealing:  Electronic  Stability 

As  discussed  in  chapter  2,  we  have  improved  the  sensitivity  of  HRXRD  to  crystalline 
imperfection  by  examining  the  FWHM  ratio  and  interference  fringes  in  addition  to 
relaxation.  We  cannot,  however,  detect  the  first  misfit  dislocations.  For  example. 
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Figure  4-6:  Cross-section  of  modulation- doped  structures  used  in  this  study.  With 
the  exception  of  the  InGaAs  channel,  all  epilayers  are  nominally  lattice-matched  to 
the  InP  substrate. 
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Figure  4-7:  HRXRD  scans  of  sample  4089,  a  MODFET  heterostructure,  before  and 
after  annealing.  The  InGaAs  channel  is  a  factor  of  2.6  thicker  than 

sample  1879  has  a  small  FWHM  ratio  and  fringes  (fig.  2.5),  but  TEM  reveals  misfit 
dislocations  (fig.  3.18);  the  density  is  too  low  to  obtain  a  reliable  density.  Hence, 
it  is  possible  that  degradation  is  occuring  during  our  high-temperature  anneals,  but 
HRXRD  cannot  detect  it. 

The  mobility  of  modulation-doped  structures  is  very  sensitive  to  layer  and  interface 
quality.^^®  To  obtain  a  measure  of  thermal  stability  which  is  more  sensitive  than 
HRXRD,  we  grew  MODFET  structures  with  mismatched  channels  and  measured  the 
electron  mobihty  before  and  after  anneals.  The  basic  structure  was  given  in  fig.  4.6. 
The  Si  dopant  atoms  are  placed  in  a  high  band-gap  InAlAs  layer.  Most  electrons 
donated  by  the  Si  diffuse  into  the  low  band-gap  InGaAs  layer  (channel)  which  is 
free  of  dopant  atoms  and  hence  has  a  high  mobihty.  In  this  way,  both  high  carrier 
concentration  and  high  mobihty  are  achieved  simultaneously. 

The  mobihty  of  MODFETs  is  a  function  of  both  the  background  carrier  concen¬ 
tration  in  the  channel  and  the  abruptness  of  interfaces.  Interdiffusion  of  In,  Ga,  and 
A1  during  anneahng  will  make  the  interfaces  less  abrupt.  In  addition.  Si  might  diffuse 
through  the  100  A  spacer  and  into  the  channel,  lowering  the  mobihty.  The  diffusion 
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of  Si  in  GaAs  has  been  shown  to  occur  via  a  Si^o-Si^,  complex.^^’'  The  diffusion 
coefficient  for  the  complex  is  2  x  10~^^  cm*/sec  at  800"C.  For  a  5”  anneal,  — 

100  A.  The  diffusion  coefficient  of  Si  in  InAlAs  is  probably  different,  so  this  is  at 
best  an  order-of-magnitude  calculation.  It  does  suggest,  however,  that  Si  diffusion 
may  be  a  problem.  We  again  need  a  control  sample  to  let  us  separate  the  effects  of 
interdiffusion  and  relaxation.^ 

We  investigated  three  MODFET  structures,  each  with  a  500  A  In*Gai_g.A6  chan¬ 
nel.  Sample  4168  was  nominally  lattice  matched.  As  expected,  the  channel  did 
not  give  a  distinct  HRXRD  peak,  but  we  are  confident  that  x  =  0.53  ±  0.01  and 
tltc,MB  <  1-  For  MODFET  4089,  HRXRD  gives  x  =  0.64  (fig.  4.7),  and  = 

2.6.  For  MODFET  4084,  ®  =  0.68  (fig.  2.10),  and  t/tc,MB  =  3.8. 

Initially,  we  tried  to  perform  successive  anneals  on  a  single  piece  of  a  sample. 
We  observed  anomalous  results,  apparently  due  to  diffusion  of  the  In  used  for  ohmic 
contacts  during  the  high-temperature  anneals  (600-800‘’C).  Instead,  we  scribed  each 
sample  into  several  ~  5  x  5  mm  pieces.  Some  pieces  were  annealed  at  high  tempera¬ 
tures  before  the  formation  of  the  ohmic  contacts  and  the  HaU/van  der  Pauw  (HVDP) 
measurements.  Others  were  only  subjected  to  a  375‘’C  contact  anneal,  presumably 
yielding  the  as-grown  mobihty,  fi,  and  sheet  carrier  concentration,  ns-  We  observed 
nonuniformities  in  fi  and  ng  of  5  to  15%  across  each  sample. 

We  plot  the  electron  mobility  and  sheet  carrier  concentration  as  a  function  of 
anneal  temperature  and  time  in  fig.  4.8  (room-temperature  HVDP  measurements) 
and  fig.  4.9  (77K).  We  measured  3-4  pieces  without  high-temperature  anneals  on  each 
wafer.  The  error  bars  shown  are  the  standard  deviations  of  these  measurements. 

For  sample  4168,  the  control,  the  changes  in  fisooK  with  annealing  are  relatively 
small  and  may  reflect  as-grown  nonuniformities.  The  changes  in  ng^sooK,  however, 
are  much  larger  than  the  error  bar.  Annealing  causes  a  large  reduction  in  ng^sooK, 
with  a  decrease  of  a  factor  of  5  after  60”  at  800"  C.  It  is  well  known  that  Si  can  be  an 
amphoteric  dopant  in  III-V  compounds.  We  speculate  that  Sij„_^j-Si]i,  complexes  are 

complete  separation  of  the  effects  may  not  be  possible  since  strain  is  known  to  enhance 
diffusion. 
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Figure  4-8:  Electron  mobility  (a)  and  sheet  carrier  concentration  (b)  at  300  K  as  a 
function  of  annealing  for  MODFET  heterostructures. 
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Figure  4-9:  Electron  mobility  (a)  and  sheet  carrier  concentration  (b)  at  77  K  as  a 
function  of  annealing  for  MODFET  heterostructures. 


forming  during  the  anneals,  hence  reducing  the  measured  electron  concentration. 

For  sample  4084,  ^isooK  was  essentially  unchanged  by  5”  anneals  at  700  and  SOO^C. 
For  the  same  anneals,  sample  4089  exhibited  decreases  in  fisooK-  After  60  at  800®C, 
however,  fisooK  =  10,300  cm^V-sec,  about  the  same  as  the  pre-anneal  values.  Also, 
we  note  that  ns, book  was  largest  (2.3  xl0”/cm^)  for  the  4089  sample  with  the  lowest 
mobility,  5”  at  800“ C.  In  fact,  for  the  three  MODFET  structures  described  here  as 
well  as  all  others  that  we  grew,  /iaoojc  was  always  less  than  7500  cm^/V-sec  when 
ns,3ooK  was  greater  than  2.1  xlO^V^m^  This,  suggests  that  the  variations  in  room- 
temperature  mobility  for  4089  are  the  result  of  wafer  nonuniformities,  not  annealing. 

The  77K  mobility  of  both  4084  and  4089  dropped  after  5”  at  800“  C,  while  the 
lattice-matched  sample,  4168,  was  not  affected.  After  60”  at  800“C  the  77K  mobil¬ 
ity  of  4089  dropped  from  52,000  to  36,000  cm^V-sec,  whereas  the  300K  mobility 
was  apparently  unchanged.  We  believe  that  the  77K  drop  is  real,  not  a  result  of 
nonuniformity. ^  After  60”  at  800“C,  fijTK  for  4168  dropped  from  40,000  to  13,500 
cm^/V-sec,  a  much  larger  decrease  than  4089  exhibited.  A  likely  cause  of  the  drop 
in  mobility  is  Si  diffusion  into  the  channel.  It  is  not  clear  why  the  effect  is  larger  for 
the  lattice  matched  structure. 

In  summary,  the  room-temperature  mobihty  of  our  pseudomorphic  MODFETs  is 
not  affected  by  anneals  up  to  at  least  60”  at  800“C.  The  77K  mobility  is  more  sensitive 
and  decreases  after  800“C  anneals.  Decreases  in  (Ittk  sdso  occur  for  a  lattice-matched 
MODFET,  suggesting  that  interdiffusion  may  be  responsible.  We  conclude  that  sig¬ 
nificant  lattice  relaxation  is  not  occuring  during  annealing.  These  results  suggest 
that  our  high-quality  epilayers  beyond  the  Matthews-Blakeslee  critical  thickness  are 
thermally  stable  and  suitable  for  use  in  devices  such  as  FETs. 

study  of  pseudomorphic  InGaAs/InAlAs  MODFETs  (without  annealing)  also  found  that 
degradation  in  mobility  appeared  first  at  77K  and  then  at  300K.®® 
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4.3  Relaxation  Kinetics 

Research  on  a  variety  of  semiconductors  including  III-V  compounds  has  demon¬ 
strated  that  relaxation  processes  are  governed  by  both  thermodynamic  and  kinetic 
constraints.  The  Matthews-Blakeslee  theory  discussed  in  chapter  1  and  appendix  C 
only  considers  thermodynamics:  if  the  lowest  (equilibrium)  energy  state  of  a  het¬ 
erostructure  includes  misfit  dislocations,  the  layer  is  said  to  exceed  the  critical  layer 
thickness.  In  MBE,  however,  layer  growth  occurs  under  conditions  which  are  far  from 
equilibrium,  resulting  in  metastable  structures.  The  degree  of  relaxation  may  be  a 
function  of  the  nucleation  and  growth  of  misfit  dislocations. 

Several  theories  have  been  proposed  to  explain  the  relaxation  kinetics  of  strained 
layers.  Intuitively,  relaxation  should  increase  with  increasing  mismatch  or  thickness. 
In  addition,  since  the  kinetics  of  dislocation  nucleation  and  growth  are  involved,  both 
time  and  temperature  will  be  factors.  If  one  is  modeling  the  relaxation  during  growth, 
the  time  required  to  grow  the  layer  and  the  growth  temperature  are  the  relevant 
parameters.  One  can  also  model  the  relaxation  during  post-growth  anneahng.  In  that 
case,  the  temperature  and  duration  of  the  anneal  are  the  appropriate  parameters. 

At  this  point,  it  is  probably  fair  to  say  that  none  of  the  relaxation  kinetics  theories 
is  widely  accepted.^^-^®-”®’"®  The  theory  by  Dodson  and  Tsao^®'^^“  has  shown  some 
success  in  applications  to  SiGe/Si®®’^^^  and  InGaAs/GaAs®^  heterostructures.  We 
will  briefly  outline  their  model  and  apply  it  to  our  data. 

Dodson  and  Tsao’s  model,  in  part  based  upon  the  earher  work  of  Alexander  and 
Haasen,^^^  considers  both  the  kinetics  of  dislocation  formation  and  the  increase  in 
dislocation  density  caused  by  dislocation  multiplication.  The  result  is  a  differential 
equation  for  7(t)>  f^e  strain  relief: 

^  =  -7(0-KM]*[7(0  +  7o]  (4-1) 

where  D  is  a  phenomenological  parameter  describing  dislocation  mobility  and  multi¬ 
plication  and  7o  represents  a  source  term  for  misfit  dislocations;  both  D  and  7o  are 
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“adjustable  parameters.”®  The  function  r(h)  is  the  homogeneous  strain  retained  by 
the  overlayer,  given  by: 

_  b(l  —  i/cos^a)ln[Ahfb) 
h47r(l  +  u) 

where  h  is  the  layer  thickness,  i/  is  Poisson’s  ratio,  a  is  an 
cation  orientation  (see  appendix  C),  and  b  is  the  Burger’s 
The  strain  rehef  in  the  Dodson-Tsao  model,  7(t),  is 
mismatch.  Since  orthorhombic  distortion  may  be  present, 

=  V  °  >'ll[no]  V  °  /llfiio] 

2 

We  can  also  relate  strain  rehef  to  the  dislocation  density  by:®^ 

7(t)  =  p{t)b 

where  p{t)  is  the  hnear  density  of  misfit  dislocations  lying  in  the  interface, 
that  (4.4)  is  strictly  vahd  only  if  misfit  dislocations  are  the  only  defects 
layer  strain.  Strain  may  also  be  reheved  by  defects  such  as  stacking  faults  and  edge 
dislocations  which  often  result  from  three-dimensional  growth  (see  chapters  1  and  6). 

Previous  work  in  the  InGaAs/GaAs  system  suggests  that  the  parameter  D  in 
the  Dodson-Tsao  model  may  be  a  function  of  lattice  mismatch.®'*  Thus,  we  apply  the 
model  to  sets  of  samples  with  constant  composition  (mismatch)  and  varying  thickness 
(see  fig.  2.19  for  our  complete  set  of  relaxation  data).  We  begin  with  a  set  of  four 
InGaAs  samples  with  (Aa/a)^  =  -f-0.010.  We  used  the  Runga-Kutta  method  of  order 
4  to  numerically  integrate  the  differential  equation,  (4.1).  We  adjusted  D  and  7o 
until  a  reasonable  fit  to  the  experimental  data  was  achieved.  The  fit  was  obtained 
with  D  =  70  see"*  and  70  =  1  x  10“®.®^  The  data  and  model  fit  are  shown  in  fig. 

*Our  D  is  equivalent  to  Dodson  and  Tsao’s  where  p  is  the  shear  modulus.  We  note  that 

the  correct  form  of  r(h)  is  given  in  Dodson  and  Tsao’s  erratum.  However,  the  equation  equivalent 
to  our  (4.1)  apparently  has  a  sign  error  in  the  erratum  but  is  correct  in  the  original  paper. 

’In  all  cases,  we  use  a  t  equal  to  the  growth  time  plus  100  seconds  to  represent  the  cool-down 
period  after  growth.  The  results  are  not  very  sensitive  to  the  t  that  is  used. 


(4.3) 

(4.4) 

We  note 
reheving 


(4.2) 

angle  describing  the  dislo- 
vector  of  the  dislocations, 
equivalent  to  the  parallel 
we  have: 
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Figure  4-10:  Parallel  lattice  mismatch  as  a  function  of  layer  thickness  for  two  sets  of 
InGaAs  layers  in  compression.  The  soUd  hnes  are  the  predicted  strain  reUef  from  the 
Dodson- Tsao  model,  eq.  (4.1). 

4.10.  The  fits  are  not  very  sensitive  to  the  value  of  7o.  We  note  that  for  a  fixed 
D,  relaxation  decreases  when  7o  increases,  contrary  to  what  one  might  expect  for  a 
dislocation  source  term.  A  second  set  of  InGaAs  in  compression  is  also  available  for 
comparison.  The  relaxed  lattice  mismatch,  0.013,  is  larger  than  the  previous  set.  We 
attempted  to  fit  the  data  using  the  same  values  of  D  and  7o.  The  predicted  curve  is 
a  reasonable  fit  to  the  data,  as  shown  in  fig.  4.10. 

We  also  used  the  Dodson-Tsao  model  to  fit  sets  of  InGaAs  in  tension  (fig.  4.11) 
and  InAlAs  in  compression  (fig.  4.12),  using  D’s  of  55  and  70  sec  respectively.  We 
note  that  these  values  of  D  are  within  a  factor  of  1.5  of  the  SiGe  results  of  Dodson 
and  Tsao.^2°  We  were  unable  to  fit  the  10,000  A  data  point  in  fig.  4.12.  This  is  the 
sample  with  R  =  114%.  The  model  will  never  predict  a  relaxation  greater  than  100%. 

In  the  case  of  InAlAs  in  tension,  a  10,000  A  layer  with  a  mismatch  of  -0.0066  is 
only  32%  relaxed.  In  order  to  obtain  a  reasonable  fit  to  this  data  set  (fig.  4.13),  we 
had  to  use  a  very  small  value  of  D,  0.05  sec-^  This  suggests  that  the  relaxation 
mechanism  is  much  different  than  the  other  three  cases.  We  wiU  return  to  this  issue 
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Figure  4-11:  Parallel  lattice  mismatch  as  a  function  of  layer  thickness  for  InGaAs 
layers  in  tension.  The  solid  line  is  the  predicted  strain  relief  from  the  Dodson-Tsao 
model,  eq.  (4.1). 


THICKNESS  (A) 

Figure  4-12:  Parallel  lattice  mismatch  as  a  function  of  layer  thickness  for  InAlAs 
layers  in  compression.  The  solid  Hne  is  the  predicted  strain  relief  from  the  Dodson- 
Tsao  model,  eq.  (4.1). 
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Figure  4-13:  Parallel  lattice  mismatch  as  a  function  of  layer  thickness  for  InAlAs 
layers  in  tension.  The  solid  line  is  the  predicted  strain  relief  from  the  Dodson-Tsao 
model,  eq.  (4.1). 

in  chapter  6. 

We  now  consider  the  results  of  our  annealing  experiments.  In  sections  4.1  and 
4.2,  we  demonstrated  the  thermal  stabihty  of  high-quality  InGaAs  and  InAlAs  layers 
beyond  the  Matthews-Blakeslee  (equilibrium)  critical  layer  thickness.  The  funda¬ 
mental  question  is:  If  the  layers  are  metastable,  why  is  no  relaxation  observed  for 
high-temperature  anneals? 

In  order  for  relaxation  to  occur,  a  source  of  misfit  dislocations  is  required.  One 
potential  source  is  threading  dislocations  from  the  substrate.  These  threading  dis¬ 
locations  can  bend  into  misfit  dislocations  at  the  epilayer/substrate  interface,  as  il¬ 
lustrated  in  fig.  4.14.  Our  InP  wafers  are  guaranteed  to  have  less  than  5  x  10^/cm^ 
threading  dislocations.  As  an  order-of-magnitude  calculation,  we  assume  2  x  10^/cm* 
threading  dislocations  and  consider  a  1  cm  by  1  cm  sample.  We  also  assume  all  the 
threading  dislocations  are  bent  into  misfit  dislocations  and  propagate  to  an  edge  of 
the  sample.  With  these  assumptions,  the  average  length  of  a  misfit  dislocation  is  0.5 
cm  and  the  linear  density,  p,  is  5000/cm  or  0.5//xm.  Using  (4.4),  we  obtain  (Aa/o)|| 
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Figure  4-14;  Behavior  of  substrate  threading  dislocations  in  strained-layer  epitaxy. 
For  layers  below  the  critical  thickness,  threading  dislocations  propagate  to  the  surface 
(a).  For  layers  above  the  critical  thickness,  a  misfit  dislocation  segment  may  form  at 
the  interface  (b). 

=  2  X  10"^  which  is  about  the  sensitivity  hmit  of  HRXRD,  but  is  a  factor  of  30-60 
smaller  than  the  mismatch  for  the  samples  we  annealed.  We  conclude  that  threading 
dislocations  do  not  provide  a  sufficient  source  of  misfit  dislocations  to  substantially 
relax  our  strained  layers. 

Misfit  dislocation  segments  can  also  be  formed  from  half-loops  propagating  from 
the  epilayer  surface.  These  half-loops  must  overcome  an  energy  barrier  to  nucleate. 
The  energetics  of  half-loop  nucleation  was  considered  by  Fitzgerald  et  for  the  case 
of  InGaAs  on  GaAs  and  Kamat  and  Hirth”®  for  GaAs  on  other  materials.  We  will 
apply  the  procedure  of  Fitzgerald  to  InGaAs/InP  and  InAlAs/InP  heterostructures. 
The  change  in  energy  due  to  the  formation  of  a  semicircular  loop  is:®^ 


E  = 


Ghr 


-r-  6(2  -  -  87rr(Aa/a),(l  -|-  v)co3Xcos4>  -  26(1  -  v)sina 

I/)  L  ' 


8(1-1/) 


(4.5) 


where  r  is  the  loop  radius,  Ce  is  the  core  energy  factor,  and  the  other  parameters  were 
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Figure  4-15:  Formation  of  a  misfit  dislocation  from  a  dislocation  half-loop.  The 
critical  radius  of  the  half-loop  is  r*. 

previously  defined.  We  maximize  this  energy  to  obtain  r*,  the  critical  loop  radius  for 
surface  nucleation.  The  result  is; 

.  _  6(2  -  u)  [ln(8(7sr*/e^6)  -f  1]  -  2fc(l  -  u)sina 
^  ~  167r(Ao/o),(l -t- i/)cosAcos(^ 

The  energy  of  the  system  will  be  at  a  maximum  if  r  =  r*.  Hence,  if  a  loop  has 
sufficient  energy  to  reach  r*,  it  can  then  spontaneously  grow  until  it  reaches  the 
interface  and  forms  a  misfit  dislocation  as  illustrated  in  fig.  4.15. 

We  solved  (4.6)  using  Newton’s  method,  assuming  Ce  —  4,^^®  A  =a  =  60°,  and  4>  = 
35.26°.^*  Other  values  were  interpolated  from  appendix  B.  The  result  for  InGaAs/InP 
is  plotted  in  fig.  4.16a.  As  the  magnitude  of  the  mismatch  increases,  the  size  of  the 
critical  radius  decreases.  Using  the  calculated  values  of  r*,  we  can  determine  the 
critical  energy,  E* ,  from  (4.5).  The  result  is  shown  in  fig.  4.16b.  As  lattice  mismatch 
increases,  E*  decreases.  In  other  words,  homogeneous  nucleation  is  more  likely  in 
highly  strained  systems.  This  makes  physical  sense  because  the  amount  of  elastic 
energy  released  by  a  half-loop  is  proportional  to  the  strain.  The  results  for  InAlAs, 
fig.  4.17,  are  similar.  The  critical  energy  for  AlAs  is  higher  than  for  InAs  or  GaAs 
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primarily  because  AlAs  has  a  larger  shear  modulus. 

We  now  consider  the  energy  available  for  the  nucleation  of  surface  loops.  Estimates 
for  the  energy  available  during  growth  are  50kT^*^  and  88kT.^^®  We  assume  that 
energy  available  during  annealing  is  given  by  the  same  expressions.  Using  88  kT  and 
our  maximum  annealing  temperature  of  1123K,  we  obtain  8.6  eV.  From  figs.  4.16 
and  4.17,  the  energy  required  for  loop  nucleation  is  always  at  least  25  eV.  For  the 
compositions  used  in  our  annealing  experiments,  E*  is  100-500  eV.  We  conclude  that 
sufficient  energy  is  not  available  for  the  homogeneous  nucleation^l  of  surface  half-loops 
during  anneals. 

Our  findings  are  consistent  with  those  of  Fitzgerald  et  and  Kamat  and 
Hirth.^^^  Fitzgerald  estimated  that  a  mismatch  of  about  5%  is  required  for  half¬ 
loop  nucleation  during  the  growth  of  InGaAs  on  GaAs  at  550^0.  Kamat  and  Hirth 
calculated  that  a  mismatch  of  at  least  9%  is  required  for  homogeneous  nucleation  on 
a  GaAs  surface  at  TSO^C. 


4.4  Summary 

In  summary,  we  have  shown  that  significant  relaxation  of  InGaAs  and  InAlAs  layers 
occurs  during  growth  when  the  mismatch  and  thickness  are  relatively  large.  Layer 
thicknesses  of  50-100  times  tc.MB  are  required  to  achieve  almost  complete  strain  re¬ 
laxation. 

We  examined  the  thermal  stabihty  and  relaxation  kinetics  of  strained  layers.  For  a 
fixed  mismatch,  the  degree  of  relaxation  during  growth  can  be  fit  by  the  kinetic  model 
of  Dodson  and  Tsao,  although  adjustable  parameters  are  required.  Using  HRXRD, 
we  observe  additional  relaxation  when  partially  relaxed  InGaAs  layers  are  annealed. 
High-quality  (unrelaxed)  InGaAs  and  InAlAs  layers  with  tftc,MB  =  3-9  do  not  exhibit 
relaxation  during  anneals  up  to  850°C.  In  addition,  the  room-temperature  mobility 
of  InGaAs/InAlAs  MODFET  heterostructures  with  pseudomorphic  channels  is  not 

II Following  the  convention  in  the  literature,  we  use  the  term  homogeneous  nucleation  to  refer  to 
the  nucleation  of  a  dislocation  half-loop  on  a  defect-  and  particulate-free  epilayer  surface. 
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Figure  4-16:  Critical  radius  (a)  and  energy  (b)  for  homogeneous  nucleation  of  dislo¬ 
cation  half-loops  in  InGaAs. 
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Figure  4-17:  Critical  radius  (a)  and  energy  (b)  for  homogeneous  nucleation  of  dislo¬ 
cation  half- loops  in  InAlAs. 
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affected  by  anneals  up  to  SOO'C.  Degradation  in  the  77K  mobility  is  observed,  but 
is  also  present  in  a  lattice  matched  MODFET.  This  suggests  that  inter  diffusion,  not 
relaxation,  is  the  primary  cause.  We  explain  the  thermal  stabiHty  of  our  epilayers  by 
the  lack  of  sufficient  energy  to  homogeneously  nucleate  half-loop  dislocations  and  the 
absence  of  a  sufficient  number  of  threading  dislocations  in  the  substrate. 


no 


Chapter  5 


Device  Application: 
Heterostructure  Field-Effect 
Transistors 


The  critical  layer  thickness  of  Matthews  and  Blakeslee,  is  often  used  as  an 

upper  hmit  in  the  design  of  pseudomorphic  device  heterostructures.*  As  we  showed 
in  chapter  2,  however,  complete  lattice  relaxation  does  not  abruptly  occur  when  an 
InGaAs  or  InAlAs  layer  reaches  tc,MB-  In  fact,  layers  often  retain  high  crystaUine 
quality  up  to  thicknesses  several  times  tc,MB-  Such  layers  are  thermally  stable  to 
temperatures  of  at  least  800°C  (chapter  4).  These  results  suggest  that  layers  thicker 
than  tc,MB  may  be  suitable  for  some  device  applications. 

We  fabricated  heterostructure  field  effect  transistors  (HFETs)  in  the  InAlAs/ 
InGaAs/lnP  materials  system.  In  this  chapter,  we  present  both  the  HRXRD  and 
device  results  for  an  HFET  in  which  one  layer  exceeds  tc.MB-  We  also  illustrate  how 
HRXRD  can  be  used  to  determine  which  samples  should  be  processed  into  devices. 
Finally,  we  compile  data  from  the  hterature  for  HFETs  with  strained  InGaAs  and 
InAlAs  layers,  and  compare  it  to  <c,AfB  and  our  HRXRD  results. 
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Figure  5-1;  Cross-section  of  doped- channel  InAlAs/InGaAs  heterostructure  field- 
effect  transistor. 

5.1  HFET  Characterization  by  HRXRD 

Our  typical  device  structure  is  shown  in  fig.  5.1  and  consists  of  (bottom  to  top): 
semi-insulating  (001)  InP(Fe)  substrate,  1000  A  InAlAs  buffer,  75  A  InGaAs  sub¬ 
channel,  100  A  Si-doped  n+-InGaAs  channel,  300  A  InAlAs  pseudo-insulator,  and 
50  A  InGaAs  cap.  We  note  that  unlike  a  MODFET,  the  dopants  are  located  in  the 
channel  in  our  device.  This  results  in  lower  channel  mobility.  However,  these  devices 
typically  operate  in  the  velocity  saturation  regime,  and  saturation  velocity  is  not 
very  sensitive  to  doping.^^®  The  doped-channel  approach  offers  several  advantages 
over  MODFETs  including  higher  breakdown  voltage  and  lack  of  transconductance 
collapse. 

In  the  past,  doped-channel  HFETs  were  fabricated  using  lattice-matched  layers  of 
InGaAs  and  InAlAs  on  InP.^^®  There  are,  however,  advantages  to  using  strained  layers. 
Drain  current,  transconductance,  and  cutoff  frequency  can  be  increased  by  employing 
InAs-rich  InGaAs  channels. If  AlAs-rich  insulating  layers  are  used,  the  band-gap 
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of  the  InAlAs  and  the  conduction  band  offset  between  the  InGaAs  and  InAlAs  are 
increased.  This  reduces  hot-electron  effects,  improves  electron  confinement  in  the 
channel,  and  increases  the  breakdown  voltage  of  the  device.*®’^*® 

Both  the  channel  and  insulator  were  intentionally  strained  in  HFET  1870.  The 
channel  was  nominally  doped  to  6  x  10^®/cm®.  Hall/van  der  Pauw  measurements  yield 
a  sheet  carrier  density  of  2.1  x  10^^/cm^  and  an  electron  mobility  of  3660  cm  /V-sec 
at  room  temperature. 

The  (004)  HRXRD  scan  of  1870  is  shown  in  fig.  5.2.  We  simulated  the  structure 
using  the  nominal  layer  thicknesses  (which  we  beheve  are  accurate  to  within  10%) 
and  adjusting  the  compositions  to  obtain  a  good  match  with  the  experimental  data. 
The  best-fit  simulated  rocking  curve  is  shown  as  curve  b)  in  fig.  5.2.  The  buffer  is  the 
thickest  layer  and  hence  gives  the  peak  with  the  largest  intensity  and  narrowest  width 
(except  for  the  InP  substrate).  We  see  that  the  composition  of  the  buffer,  j/i,  =  0.506, 
is  close  to  the  lattice-matched  value,  0.521.  The  300  A  insulator  layer  also  gives  a 
distinct  but  broader  peak,  with  j/i  =  0.411.  The  100  A  channel,  75  A  subchannel,  and 
50  A  cap  are  too  thin  to  produce  distinct  peaks  (see  section  2.4).  The  structure  in  the 
A0  =  -2000  to  -200  arc-sec  range  is  primarily  caused  by  the  nominally  x=0.65  channel, 
but  we  cannot  obtain  an  accurate  measure  of  x  from  this  sample.  To  demonstrate 
the  sensitivity  of  HRXRD  to  the  buffer  and  insulator  compositions,  we  simulated  a 
structure  in  which  the  buffer  and  insulator  have  In  As  mole  fractions  which  are  0.005 
and  0.010  less  than  the  best-fit  values.  The  result  is  shown  as  curve  c)  in  fig.  5.2. 
Clearly,  changes  in  composition  of  this  size  can  be  detected  by  HRXRD. 

The  InAlAs  pseudoinsulator  layer  in  sample  1870  is  a  factor  of  2.0  thicker  than 
tc, mb-  The  experimental  FWHM  of  this  layer  is  520  arc-sec,  compared  to  500  arc-sec 
for  the  simulation.®®  This  good  agreement,  along  with  the  presence  of  Pendellosung 
fringes  in  the  experimental  curve,  indicates  a  coherent  heterostructure.  The  100  A 
Ino.65Gao.3s As  channel  is  thinner  than  icMB-  Based  upon  our  studies  and  the  recent 
work  of  Gendry  et  al.,^®’'*'^  we  believe  that  higher  In  As  mole  fractions  (or  thicker 
layers)  could  be  used  in  the  channel  without  degradation  from  lattice  relaxation  or 
three-dimensional  growth.  The  resulting  devices  would  be  expected  to  have  higher 
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Figure  5-2:  Experimental  HRXRD  rocking  curve  for  HFET  1870;  peaks  from  the 
InAlAs  buffer  and  pseudo-insulator  layers  are  visible,  b)  Best-fit  simulation  using 
structure  of  fig.  5.1.  c)  Simulation  of  heterostructure  with  different  buffer  and  pseudo¬ 
insulator  compositions. 

transconductance  but  lower  breakdown  voltage  than  1870. 


5.2  Device  Characteristics 

The  ultimate  test  of  a  strained-layer  heterostructure  is  its  device  characteristics.  In 
summary,  devices^  on  sample  1870  with  a  1.9  fim  gate  length  {Lg)  and  a  200  ^m  gate 
width  have  a  peak  transconductance  of  190  mS/mm,  a  maximum  drain  current  of  302 
mA/mm,  and  a  power  density  of  0.40  W/mm.^^®’"”  The  output  I-V  characteristics 
are  shown  in  fig.  5.3.  The  device  exhibits  complete  channel  pinch-off  and  a  breakdown 
voltage  of  14  V. 

The  unity  current  gain  cut-off  frequency,  ft,  and  unity  power  gam  cut-off  frequency, 
are  plotted  as  a  function  of  gate-source  voltage  in  fig.  5.4.  The  maximum  values 
are;  ft  =  15  GHz  and  fmax  =  101  GHz.  For  this  device,  we  calculate  an  electon 
velocity  of  1.8  x  10^  cm/sec  and  an  ft  x  Lg  product  of  28  GHz-^m.  These  values  are 

t  Device  processing  and  testing  was  performed  by  Sandeep  Bahl  at  MIT. 
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Figure  5-3:  Output  I-V  characteristics  for  HFET  1870.  The  gate  length  was  1.9  /in 
and  the  gate  width  was  200  fim. 

comparable  to  those  obtained  on  lattice- matched  MODFETs  of  similar  gate  length, 
but  the  power  density  and  breaikdown  voltage  of  our  device  are  2-3  times  higher. 

We  also  measured  device  characteristics  as  a  function  of  angular  orientation  on  the 
surface  and  observed  no  asymmetries.  The  presence  of  asymmetries  is  an  indication 
that  misfit  dislocations  have  formed.*®  We  conclude  that  if  dislocations  are  present 
in  HFET  1870,  their  density  is  insufficient  to  degrade  device  performance  in  any 
appreciable  way. 

5.3  Heterostructure  Evaluation 

In  addition  to  measuring  layer  composition,  HRXRD  can  be  used  to  determine  which 
HFET  structures  should  be  processed  into  devices.  As  an  example,  we  show  the 
HRXRD  rocking  curves  for  two  samples  in  fig.  5.5.  Both  samples  have  a  structure 
similar  to  1870  (fig.  5.1).  For  HFET  4164  (a),  the  nominal  channel  composition  has 
been  increased  from  x  ==  0.65  to  x  =  0.70.  The  channel  is  now  clearly  visible  in  the 
HRXRD  scan.  By  comparison  to  simulations,  we  confirmed  that  x  =  0.70  i  0.02.  In 
contrast,  we  do  not  observe  distinct  peaks  for  the  channel,  insulator,  or  buffer  layers 
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Figure  5-4:  Unity  current  and  power  gain  cut-off  frequencies  for  HFET  1870.  The 
drain-source  voltage  is  fixed  at  4.5  V. 

of  HFET  4143  (fig.  5.5b).  The  unintentional  mismatch  of  the  buffer  layer  normally 
produces  a  distinct  peak.  Since  it  is  the  first  layer,  the  absence  of  a  buffer  layer 
peak  for  sample  4143  suggests  that  something  went  wrong  with  the  growth  from  the 
beginning.  One  possibility  is  that  the  growth  temperature  was  too  low  and  the  native 
oxide  was  not  completely  desorbed.  Based  upon  the  HRXRD  results,  4164  should  be 
processed  but  4143  should  not.  This  conclusion  is  supported  by  mobility  measure¬ 
ments.  The  room-temperature  electron  mobility  of  sample  4164  is  4110  cm’/V-sec, 
about  what  one  would  expect  for  an  InGaAs  layer  doped  at  6  x  10^  /cm  .  The 
mobihty  for  sample  4143  is  only  230  cm’/V-sec,  confirming  that  the  sample  is  poor. 

We  have  not  observed  clear  differences  in  the  HRXRD  scans  between  “good 
HFET  (/1300K  ^  4000  cm*/V-sec)  and  “mediocre”  HFET  {fiaooK  ^  2000  cm*/V-sec) 
heterostructures.  We  also  note  that  our  biggest  problem  in  device  growth  has  been 
achieving  the  desired  doping  level.  Hence,  we  routinely  perform  both  HRXRD  and 
Hall/ van  der  Pauw  measurements  on  all  device  heterostructures.  Each  measurement 
can  be  completed  in  less  than  one  hour.  We  use  the  results  to  determine  which  samples 
to  process  and  to  make  appropriate  adjustments  to  Ga,  Al,  and  Si  cell  temperatures 
for  the  next  growths. 
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Figure  5-5:  HRXRD  rocking  curves  for  HFETs  with  structures  similar  to  fig.  5.1. 
Sample  4164  (a)  exhibits  distinct  peaks  for  the  buffer,  insulator,  and  channel  layers, 
but  4143  (b)  does  not. 
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5.4  Survey  of  Pseudomorphic 

In*Gai_a,As/Iiij,Ali_yAs/InP  HFETs 

Our  results  suggest  that  high-performance  HFETs  can  include  layers  exceeding  the 
Matthews-Blakeslee  critical  layer  thickness.  To  confirm  this  finding,  we  surveyed  the 
hterature  for  HFETs  using  mismatched  layers  of  either  InAs-rich  In*Gai_*  As  or  AlAs- 
rich  In„Ali_yAs  on  In  most  reports,  the  channel  thickness 

does  not  exceed  tc,MB,  and,  as  expected,  the  devices  show  no  sign  of  misfit  dislocations. 
In  some  studies,  however,  devices  were  fabricated  on  a  series  of  heterostructures 
with  increasing  mismatch,  and  a  degradation  in  performance  was  observed  when  the 
mismatch  was  too  large.  In  fig.  5.6  we  plot  layer  thickness  versus  mismatch  for 
InP-based  HFETs  with  InGaAs  (a)  and  InAlAs  (b)  strained  layers.  Closed  symbols 
indicate  device  structures  showing  no  signs  of  dislocations  and  open  symbols  represent 
degraded  devices.  Fig.  5.6  includes  the  calculated  values  of  tc.MB  (see  Appendix  C). 

All  the  device  results  in  fig.  5.6  suggest  that  HFETs  can  accomodate  strained 
layers  up  to  2-4  times  te,MB  without  degradation.  We  also  include  hnes  showing  esti¬ 
mates  for  the  limit  of  high  crystalline  quality  based  upon  our  experimental  HRXRD 
criteria  from  chapter  2  (FWHM  ratio  <  1.8).  We  observe  reasonable  agreement  be¬ 
tween  the  HRXRD  results  and  the  hmited  device  data  available.  For  mismatches 
greater  than  ±1%  we  use  a  dotted  line  because  of  our  limited  HRXRD  data  in  that 
region. 

In  the  InGaAs/GaAs  system,  a  recent  paper^°  reported  high-performance  metal- 
semiconductor  field-effect  transistors  in  which  the  InGaAs  layers  were  2  and  4  times 
tcMB-  The  author  concluded  that  the  mismatched  layers  are  strain-relaxed  because, 
“As  the  InGaAs  layer  grows  thicker  than  the  critical  thickness,  the  surface  lattice 
constant  increases  abruptly  and  becomes  close  to  the  bulk  InGaAs  lattice  constant.” 
Our  works  suggests  an  alternative  interpretation:  the  InGaAs  layers  are  primarily 
strained  with  a  low  density  of  misfit  dislocations  which  is  insufficient  to  appreciably 
degrade  device  performance. 

As  mentioned  in  chapter  1,  Temkin  et  alP  examined  the  reverse-bias  leakage 
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Figure  5-6:  Layer  thickness  versus  mismatch  for  InGaAs  (a)  and  InAlAs  (b).  Symbols 
represent  InAlAs/InGaAs  HFET’s  with  (open)  and  without  (filled)  signs  of  degrada¬ 
tion  due  to  mismatch.  Estimates  of  critical  thickness  from  the  anisotropic  extension 
of  the  Matthews-Blakeslee  theory  are  shown  as  solid  lines.  We  also  include  lines 
showing  our  HRXRD-based  limit  (FWHM  ratio  <  1.8)  for  high  crystalline  quality; 
the  lines  are  dotted  for  mismatches  greater  than  ±  1%  because  of  our  limited  data  in 
that  region. 
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current  in  In»Gai_*AsInP  quantum-well  p-i-n  photodiodes  and  found  an  apparent 
critical  thickness  of  1.0-1.5  Because  of  their  larger  mesa  area,  such  structures 

may  be  more  sensitive  to  misfit  dislocations  than  HFETs.  In  general,  the  amount 
of  strain  which  can  be  tolerated  in  a  device  heterostructure  may  depend  upon  the 
type  of  device.  For  example,  a  small  density  of  misfit  dislocations  may  not  affect  an 
HFET,  but  the  same  dislocation  density  in  the  active  region  of  a  laser  could  severely 
degrade  its  performance. 

5.5  Summary 

In  summary,  HRXRD  experimental  rocking  curves  combined  with  simulations  can 
be  used  to  determine  the  composition  of  epitaxial  layers  in  device  heterostructures. 
HRXRD  also  gives  a  qualitative  measure  of  the  material  quality  in  a  device  structure. 
This  information  can  be  used  to  select  samples  for  device  processing.  In  chapter  2,  we 
found  that  high  crystalline  quality,  as  measured  by  HRXRD,  can  be  maintained  to 
thicknesses  beyond  In  this  chapter,  we  have  used  actual  device  measurements  to 

demonstrate  that  layers  thicker  than  Umb  can  be  used  in  high-performance  HFETs. 
Our  device  results  are  in  general  agreement  with  other  reports  for  HFETs  in  the 
InAlAs/InGaAs/InP  materials  system. 
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Chapter  6 


Discussion 


In  this  chapter,  we  first  consider  the  possibihty  that  composition  modulation  may 
affect  the  crystalline  quality  and  lattice  relaxation  in  strained  layers.  We  then  explore 
the  connection  between  heterogeneous  dislocation  nucleation,  strain  relief,  thermal 
stability,  and  growth  modes.  We  explain  the  thermal  stability  of  our  InGaAs/InP 
and  InAlAs/InP  pseudomorphic  heterostructures  as  weU  as  conflicting  reports  in  the 
literature  for  the  InGaAs/GaAs  and  SiGe/Si  systems.  Based  upon  our  findings,  we 
suggest  procedures  for  the  growth  of  high-quahty,  thermally-stable  strained  epilayers. 
In  addition,  we  consider  the  cause  of  asymmetry  in  misfit  dislocation  density  for  our 
samples. 


6.1  Composition  Modulation 

As  discussed  in  section  1.5,  others  have  observed  composition  modulation  in  both 
InGaAs  and  InAlAs  layers.  In  chapter  3,  we  showed  that  composition  modulation 
was  in  the  wrong  direction  to  account  for  the  large  optical  anisotropy  observed  for 
many  strained  layers.  Here,  we  will  address  the  role  of  composition  modulation  in 
the  crystalline  quality  and  lattice  relaxation  of  strained  layers. 

The  period  of  composition  modulation  is  typically  a  few  hundred  to  a  few  thousand 
angstroms  compared  to  our  HRXRD  spot  size  of  about  1x2  mm.  Hence,  we  expect 
the  x-ray  beam  to  sample  many  periods  of  oscillation.  Substantial  variations  in  lattice 


121 


parameter  would  clearly  broaden  the  HRXRD  layer  peak. 

We  do  not  know  of  any  measurements  of  the  modulation  amplitude  for  InAlAs  or 
InGaAs.  In  the  case  of  In,Gai_,As„Pi_„,  variations  in  ®  and  y  are  as  large  as  0.1.®« 

In  such  a  quaternary  system,  the  lattice  parameter  can  remain  constant  as  both  x 
and  y  change.  In  our  ternary  layers,  however,  any  change  in  composition  also  changes 
the  lattice  constant.  Contrast  modulation  with  a  period  of  365  A  was  observed  for 
sample  3099,  1000  A  of  Ino.579Alo.42iAs  (fig.  3.15).  The  HRXRD  scan  for  3099  (fig. 
2.8)  exhibited  interference  fringes  and  a  FWHM  ratio  of  1.3.  We  speculate  that  the 
amplitude  of  the  modulation  in  sample  3099  is  very  smaU  (Ay  «  0.01)  and  hence 

does  not  affect  the  HRXRD  measurements. 

In  section  2.3,  we  observed  that  asymmetric  broadening  of  the  (004)  HRXRD 
layer  peaks  occurs  in  some  cases.  The  peak  widths  were  a  maximum  in  the  [110]  and 
[liO]  directions  and  a  minimum  in  the  [liO]  and  [110]  directions.  The  observed  com¬ 
position  modulations  are  in  the  <  100  >-type  directions.^^’®^'^"®  Hence,  composition 

modulation  cannot  explain  the  asymmetric  peak  broadening. 

We  have  TEM  micrographs  for  three  samples  with  degraded  crystalline  quality: 
1530  (fig.  3.16),  1880  (fig.  3.19)  and  1486  (not  shown).  In  each  case,  composition 
modulation  is  not  visible.  It  is  possible,  however,  that  composition  modulation  is 
present  but  is  obscured  by  the  extreme  surface  roughness.  Hence,  we  cannot  com¬ 
pletely  rule  out  the  possibility  that  composition  modulation  is  a  factor  in  HRXRD 

peak  broadening  and  fringe  loss  for  some  samples. 

Composition  modulation  can  lower  the  total  energy  of  a  strained  epitaxial  layer. 
Unlike  crystalline  defects  such  as  misfit  dislocations  and  stacking  faults,  however,  it 
cannot  reUeve  a  significant  amount  of  the  lattice  strain.  Peiro  et  investigated 
layers  of  Ino.54Gao.46As  on  InP  by  TEM  and  HRXRD.  For  samples  which  were  1- 
2  X  they  observed  composition  modulation,  but  no  crystalline  defects  or  lattice 

relaxation.  Samples  3-7  x  tc,MB  exhibited  composition  modulation,  high  densities  of 
stacking  faults,  and  substantial  lattice  relaxation.  Based  upon  these  results  and  our 
measurements  on  sample  3099,  we  believe  that  composition  modulation  alone  is  not 
necessarily  harmful  to  epitaxial  layers.  Composition  modulation  might  influence  the 
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point  at  which  misfit  dislocations,  surface  roughness,  or  3-D  growth  begins. 

Finally,  we  note  that  Yao  et  used  TEM  to  observe  the  transition  from  2-D 
to  3-D  growth  in  the  InGaAs/GaAs  system.  They  did  not  see  any  evidence  of  com¬ 
position  modulation.  All  the  layers  they  investigated  had  a  mismatch  of  at  least  1%, 
compared  to  less  than  0.1%  for  Peiro’s  InGaAs  layers.  This  suggests  that  composition 
modulation  may  be  more  likely  in  systems  with  small  mismatch  because  3-D  growth 
does  not  occur. 


6.2  Heterogeneous  Nucleation  of  Dislocations 

In  chapter  4,  we  demonstrated  the  thermal  stability  of  high-quality  InGaAs  and 
InAlAs  layers  which  exceeded  the  Matthews-Blakeslee  critical  layer  thickness.  We 
suggested  that  the  layers  were  stable  because  of  1)  an  insufficient  number  of  thread¬ 
ing  dislocations  to  bend  into  misfit  dislocations  and  relieve  strain  and  2)  insufficient 
thermal  energy  to  homogeneously  nucleate  half-loop  dislocations.  We  also  showed 
an  example  from  the  literature  in  which  a  comparable  SiGe  layer  relaxed  substan¬ 
tially  during  annealing.  Calculations  show  that  even  more  thermal  energy  is  required 
to  homogeneously  nucleate  half-loops  on  SiGe  than  on  InGaAs  or  InAlAs.  Hence, 
homogeneous  half-loop  nucleation  cannot  explain  the  difference  in  thermal  stability. 

A  survey  of  the  SiGe  literature  reveals  conflicting  results  for  thermal  stability  with 
some  layers  relaxing  during  anneals  of  600-800® as  illustrated  in  fig.  4.4, 
and  others  remaining  stable  to  at  least  800®C.^^®-^^®  A  few  annealing  experiments  have 
also  been  reported  in  the  InGaAs/GaAs  system.  Peercy  et  a/.^°  observed  a  drastic 
reduction  in  photoluminescence  intensity  when  device  heterostructures  with  InGaAs 
layers  thicker  than  ^c.Mn  were  annealed.  Bertolet  et  annealed  similar  samples 
and  found  no  changes  in  the  photoluminescence.  Bertolet  et  al.  point  out  that  their 
structures  were  grown  by  MOCVD  whereas  Peercy  et  al.  used  MBE.  They  suggest 
that  the  MOCVD  layers  are  more  stable  because  of  the  higher  growth  temperature. 
Watson  et  examined  the  stability  of  MBE-  and  MOCVD-grown  InGaAs  on 

patterned  GaAs  by  counting  dark-line  defects  imaged  by  cathodoluminescence.  They 
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found  the  MOCVD  layers  to  be  much  more  stable,  and  attributed  the  dilference  to 
heterogeneous  nucleation  sources,  with  more  sources  available  in  the  MBE  layers. 

Our  growth  temperatures  were  sUghtly  lower  than  the  temperatures  used  by 
Peercy  et  al.^°  (SbO^C)  and  Watson  et  (520‘’C)  for  MBE-growth  of  InGaAs 

on  GaAs.  Our  samples  do  not,  however,  show  signs  of  relaxation  during  annealing.^ 
Based  upon  this,  we  beheve  the  growth  temperature  is  not  the  primary  factor  deter- 
mining  thermal  stability. 

Differences  in  heterogeneous  nucleation  sources  can  explain  our  data  as  well  as  the 
InGaAs/GaAs  and  SiGe/Si  results.  Apparently,  our  samples  did  not  have  a  sufficient 
number  of  sources  to  nucleate  half  loops  to  relieve  strain  and  degrade  the  HRXRD 
rocking  curves  and  electron  mobility.  Since  our  samples  were  grown  by  MBE,  we 
conclude  there  is  nothing  inherent  in  the  MBE  process  which  causes  high  densities  of 
nucleation  sources. 

We  now  consider  the  origin  of  the  heterogeneous  dislocation  sources.  Dodson 
suggests  that  dislocation  half-loops  will  be  generated  to  relieve  any  local  stress  con¬ 
centration.  The  stress  concentrations  could  result  from  clusters  of  impurities  such 
as  oxygen,  nitrogen,  or  carbon.  Eaglesham  identified  “diamond  defects”  as  a  hetero¬ 
geneous,  regenerative  nucleation  source  in  SiGe  films.  These  defects  are  produced 
by  the  dissociation  of  stacking  faults.^^^  Another  possible  source  of  nucleation  sites  is 
oval  defects  which  are  thought  to  be  caused  by  Ga-spitting  (in  MBE)  and  particulates 
on  the  surface. This  suggests  that  sample  preparation,  transfer,  and  outgassing 
may  influence  the  stabihty  of  strained  layers.  Our  use  of  epi-ready  substrates  without 
additional  cleaning  (see  appendix  A)  could  be  important. 

Point  defects  might  also  be  a  significant  source  of  misfit  dislocations.  In  the  case 
of  intentionaUy  doped  layers,  the  dopant  atoms  themselves  are  probably  the  most 
abundant  source  of  point  defects.  Watson  et  al.^^^  noted  that  their  layers  were  n- 
type,  with  sulfur  as  the  dopant  in  the  MOCVD  material  and  sihcon  in  the  MBE 
layers.  They  suggest  that  the  Si  atoms  might  be  better  nucleation  sources.  In  the 

tWe  are  using  different  characterization  techniques  than  Watson  and  Peercy,  but  the  post¬ 
annealing  dislocation  densities  on  MBE-grown  samples  observed  by  Watson  are  much  larger  than 
the  detection  limit  for  our  HRXRD  measurements. 
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case  of  the  MODFETs  we  annealed,  the  Si  atoms  were  in  the  lattice-matched  InAlAs, 
not  the  mismatched  channel.  Hence,  Si  should  not  play  a  role  unless  it  diffused 
into  the  channel.  The  strained  layers  in  samples  1879  (fig.  4.5)  and  1682  (fig.  4.6), 
however,  were  intentionally  Si-doped,  with  nominal  concentrations  of  2xl0^'^/cm®  and 
6  X  10^^/cm®,  respectively.  These  samples  did  not  relax  during  anneals  up  to  850  C. 
We  conclude  that  substitutional  Si  atoms  are  not  an  efficient  source  of  nucleation 
sites  at  these  concentrations. 

6.3  Asymmetries  in  Misfit  Dislocation  Density 

In  chapter  2,  we  used  two  HRXRD  techniques  to  demonstrate  asymmetric  lattice 
relaxation  in  both  InGaAs  and  InAlAs  in  compression  or  tension.  We  are  now  in  a 
position  to  explain  the  asymmetries.  Fox  and  Jesser  performed  a  detailed  analysis  of 
dislocation  asymmetries  in  MOCVD-grown  GaAso.95Po.05  on  GaAs.  They  consid¬ 
ered  four  possible  sources:  substrate  misorientation,  thickness  gradients,  differences 
in  nucleation  energies,  and  differences  in  Peierls  barriers.  We  wiU  evaluate  each  of 
these  sources  for  our  samples. 

1)  Substrate  Misorientation:  All  of  our  substrates  were  guaranteed  to  be 
within  0.5"  of  the  (001)  orientation.  It  seems  unlikely  that  such  small  misorienta- 
tions  could  cause  substantial  asymmetries.  Furthermore,  we  have  observed  the  same 
asymmetry  pattern  on  substrates  from  many  different  InP  boules.  If  we  assume  that 
the  misorientations  are  randomly  distributed  in  direction,  the  asymmetry  patterns 
should  vary  from  one  substrate  to  another.  We  conclude  that  substrate  misorienta¬ 
tion  is  not  the  cause  of  the  observed  asymmetries. 

2)  Thickness  Gradient:  Unlike  the  MOCVD  layers  of  Fox  and  Jesser,  our  MBE 
layers  do  not  have  significant  thickness  gradients.  Hence,  we  rule  this  out  as  a  cause 
of  asymmetries. 

3)  Nucleation  Energy:  Dislocations  along  the  [110]  and  [110]  directions  in  zinc- 
blende  semiconductors  have  different  dislocation  cores.  Abrahams  et  al.  illustrated 
the  difference  and  we  show  their  drawing  as  fig.  6.1.^^  In  case  b)  of  fig.  6.1,  the  core 
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Figure  6-1:  a)  Atom  positions  in  three  unit  cells  of  the  zinc-blende  structure.  Atoms 
labelled  A  represent  In,  Ga,  or  A1  in  our  materials.  Atoms  labelled  B  represent  As^ 
b)  Core  of  60“  dislocation  parallel  to  [Oil];  the  slip  plane  is  (111),  c)  Core  of^O 
dislocation  parallel  to  [OU];  the  slip  plane  is  (Ill),  (figure  from  Abrahams  et  al.  ) 


of  the  60“  dislocation  is  paraUel  to  the  [Oil]  direction;  in  c),  it  is  parallel  to  [Oil].  The 
dislocation  motion  occurs  on  the  group  III  or  group  V  ghde  set  of  planes.  Smce  the 
dislocation  cores  are  different,  we  expect  different  energy  barriers  for  the  two  types 
of  dislocations. 

Fox  and  Jesser  argued  that  in  their  samples,  sufficient  threading  dislocations  were 
present  to  account  for  the  observed  misfit  dislocations.  Hence,  nucleation  of  additional 
dislocations  was  not  required.  As  we  showed  above  and  in  chapter  4,  nucleation  of 
half-loops  is  necessary  to  substantially  relax  our  layers.  Thus,  we  cannot  rule  out 
the  possibility  that  differences  in  nucleation  energy  contribute  to  the  orthorhombic 
distortion  shown  in  fig.  2.20. 

4)  Peierls  Barriers:  Based  upon  the  influence  of  doping  on  the  dislocation 
asymmetries,  observations  of  dislocation  motion  in  bulk  GaAs,  and  the  elimination  of 
the  other  three  sources.  Fox  and  Jesser  concluded  that  differences  in  Peierls  barriers 
were  responsible  for  the  asymmetries  in  misfit  dislocation  density.  We  have  TEM 
images  of  InGaAs  in  tension  (fig.  3.16)  and  InAlAs  in  tension  (not  shown)  which 
reveal  misfit  dislocations  along  the  [lIO]  direction  only,  with  an  average  spacing  of 
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about  1  fim.  These  relatively  low  densities  of  misfit  dislocations  can  be  attributed 
to  threading  dislocations.  Hence,  the  Fox-Jesser  analysis  apphes  for  samples  with 
low  dislocation  densities  and  negligible  relaxation.  Differences  in  the  Peierls  barriers 
could  be  responsible  for  the  observed  asymmetry  in  such  samples. 

The  preferred  direction  for  dislocations  is  consistent  in  our  TEM  and  HRXRD 
measurements.  Assuming  a  difference  in  Peierls  barriers  exists  to  account  for  the 
asymmetries  in  the  samples  with  a  low  dislocation  density,  the  difference  in  barriers 
should  also  be  present  in  the  layers  which  have  substantially  relaxed,  presumably 
from  heterogeneous  nucleation  of  dislocation  half-loops.  Hence,  the  Peierls  barrier 
differences  could  also  cause  asymmetries  at  high  dislocation  densities.  Differences  in 
nucleation  energy  may  also  contribute  to  the  asymmetry. 


6.4  Surface  Roughness  and  Relaxation 

We  have  proposed  that  heterogeneous  nucleation  sources  are  required  for  substantial 
lattice  relaxation.  We  explained  the  thermal  stabihty  of  our  high-crystalline-quality 
strained  layers  by  the  lack  of  these  sources  in  our  MBE  process.  We  also  observed 
significant  lattice  relaxation  in  many  as-grown  samples  (fig.  2.19).  This  raises  the 
question:  What  is  the  dislocation  source  for  the  relaxed  epilayers? 

In  chapter  3,  we  observed  large  optical  anisotropies  before  the  onset  of  significant 
relaxation  in  most  cases.  We  attributed  the  optical  anisotropies  to  surface  roughness, 
possibly  associated  with  3-D  growth.  HRXRD  and  TEM  results  suggest  that  this 
roughness  results  in  the  initial  loss  of  Pendellosung  fringes  and  broadening  of  the 
layer  peak,  although  misfit  dislocations  are  probably  a  factor  as  well.  We  propose 
that  the  roughness  provides  a  source  of  crystalline  imperfections  which  serve  as  sites 
for  subsequent  heterogeneous  nucleation  of  misfit  dislocations  during  growth  or  an¬ 
nealing.  This  model  is  supported  by  the  work  of  Guha  et  al.^^  who  found  that  large 
numbers  of  defects  are  introduced  at  the  boundaries  between  islands.  We  observed  a 
high  density  of  stacking  faults  in  plan-view  TEM  of  sample  1486  (not  shown),  a  6000 
A  layer  of  Ino.46oGao.s4o As  with  R  =  52%.  These  stacking  faults  reheve  strain  and 
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may  also  serve  as  a  source  of  diamond  defects,  as  suggested  by  Elaglesham  ct  dl. 

Based  upon  the  as-grown  characteristics  of  epilayers  and  their  thermal  stabihty, 
we  propose  classifying  layers  into  four  categories.  We  show  the  regions  on  a  plot  of 
thickness  versus  composition  for  InGaAs/InP  in  fig.  6.2.  In  region  I,  the  layers  grow 
in  a  2-D  mode  with  no  misfit  dislocations  and  hence  no  relaxation  (A  in  fig.  1.2).  The 
HRXRD  characteristics  are  good,  with  peak  widths  near  the  theoretical  value  and,  if 
the  layer  is  thick  enough,  Pendellosung  fringes.  The  optical  anisotropy  is  small  and 
the  layers  are  thermally  stable.  The  only  change  in  region  II  is  that  a  low  density  of 
misfit  dislocations  is  present,  probably  formed  from  threading  dislocations  from  the 
substrate  (B  in  fig.  1.2).  The  relaxation  is  no  longer  zero,  but  still  very  small.  The 
lines  dividing  regions  I  and  II  are  at  tc.MB-  This  boundary  is  difficult  to  determine 
experimentally;  we  use  tc,MB  based  upon  the  work  of  others. 

Going  from  region  II  to  region  III,  the  sample  morphology  and  x-ray  character¬ 
istics  become  poor.  In  addition,  the  optical  anisotropy  changes  from  small  to  large. 
We  expect  the  layers  to  be  thermally  unstable,  but  we  have  not  verified  this  exper¬ 
imentally.  We  separate  the  regions  by  lines  at  4  x  tc.MB-  We  use  this  value  based 
upon  our  HRXRD  and  VAAE  measurements.  It  is  the  same  as  the  C  lines  in  fig. 
2.13a.  We  note,  however,  that  based  upon  the  HRXRD  FWHM  ratio,  we  found  a 
transition  region  from  3  to  9  x  UmB'  Hence,  the  boundaries  between  regions  in  fig. 
6.2  are  only  approximate.  It  is  hkely  that  they  depend  on  variables  such  as  growth 
temperature. 

Region  IV  represents  samples  which  have  undergone  substantial  relaxation  (R  > 
10%).  The  strain  is  reheved  by  a  high  density  of  misfit  dislocations  and  possibly 
other  defects  such  as  stacking  faults.  The  layers  may  relax  further  if  annealed  at  high 
temperatures.  The  optical  anisotropy  is  small,  perhaps  because  the  surface  morphol¬ 
ogy  is  no  longer  dominated  by  elongated  islands.  We  used  a  value  of  10  x  tc.MB 
separate  regions  III  and  IV.  Again,  this  is  a  rough  approximation  (same  as  curves  R 
in  fig.  2.13),  based  upon  our  HRXRD  and  VAAE  data. 

For  clarity,  we  do  not  show  most  of  our  data  points  in  fig.  6.2;  the  reader  is 
referred  to  fig.  2.13.  We  do,  however,  include  the  samples  on  which  we  measured 
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THICKNESS 


Figure  6-2:  Proposed  regions  based  upon  material  quality  and  thermal  stability  of 
our  InGaAs  layers  on  InP.  The  soUd  hnes  are  at  tc.MB,  Stc.MB,  and  12tc,MB,  and  are 
only  approximate  boundaries.  The  descriptions  in  regions  II,  III,  and  IV  apply  to 
both  tension  and  compression;  they  are  separated  because  of  limited  space  on  the 
plot.  Samples  used  in  thermal  stabiUty  experiments  are  indicated.  MDs  =  misfit 
dislocations,  OA  =  optical  anisotropy.  Morph  =  morphology. 
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thermal  stability  in  chapter  4.  Samples  1331  and  4168  are  in  region  I  and,  based 
upon  HRXRD  and  Hall  mobihty  measurements,  respectively,  are  thermally  stable. 
Samples  4089  and  4084  in  region  II  are  also  thermally  stable.  In  our  model,  these 
samples  did  not  relax  because  their  surfaces  did  not  roughen  during  growth.  Sample 
1879  is  in  region  III,  but  was  thermally  stable  based  upon  HRXRD  measurements 
(fig.  4.4).  Its  morphology  and  optical  anisotropy  are  also  consistent  with  region  II 
samples.  We  note  that  this  layer  is  in  the  HRXRD  transition  region,  with  ijicMB  = 
7.3.  Sample  1441  is  well  into  region  IV  and  exhibits  all  the  expected  characteristics. 
The  additional  relaxation  during  annealing  of  this  layer  can  be  explained  by  the 
nucleation  of  additional  half-loop  dislocations.  Another  possibility  is  that  existing 

misfit  dislocations  lengthened  during  the  annealing  cycle. 

Much  of  fig.  6.2  may  also  apply  to  layers  of  InAlAs  on  InP.  As  shown  in  fig.  2.13b, 
HRXRD  measurements  can  be  used  to  separate  layers  into  regions  similar  to  II,  III, 
and  IV  in  fig.  6.2.  There  was  a  difference  in  tension  and  compression,  with  layers 
in  tension  retaining  high  crystalline  quality  longer.  InAlAs  layers  in  tension  did  not 
usually  exhibit  large  optical  anisotropy.  In  addition,  layers  in  tension  did  not  relax 
as  much  as  comparable  layers  in  compression  (figs.  4.12  and  4.13).  We  speculate 
that  the  evolution  of  surface  roughness  is  much  different  in  tension  and  compression, 
possibly  due  to  differences  in  cation  mobilities.  The  difference  could  result  in  fewer 
heterogeneous  nucleation  sites  and  hence  less  lattice  relaxation  in  tension. 

For  an  alternative  interpretation  of  the  InAlAs  relaxation  data,  we  return  to  the 
Dodson-Tsao  model  discussed  in  chapter  4.  For  InAlAs  in  tension,  the  value  of  the 
fitting  parameter  D  was  about  three  orders  of  magnitude  lower  than  D's  for  InAlAs 
in  compression  or  InGaAs  in  tension  or  compression.  D  is  a  function  of  dislocation 
mobility.  One  possibility  is  that  at  the  relatively  high  densities  of  dislocations  required 
for  significant  relaxation,  dislocations  are  pinning  each  other,  reducing  the  average 
dislocation  mobility.  This  phenomenon  is  also  known  as  work  hardening.^^®  It  is  not 
clear  why  dislocation  pinning  should  be  different  in  tension  and  compression,  but  it 
could  be  related  to  differences  in  surface  roughening  and  the  types  of  dislocations 

present. 
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In  section  6.1,  we  considered  the  role  of  heterogeneous  nucleation  in  lattice  relax¬ 
ation.  Based  upon  our  proposed  model,  we  now  discuss  techniques  which  can  be  used 
to  grow  high-quality,  thermaUy-stable  strained  layers. 

We  know  of  no  reports  of  misfit  dislocations  or  3-D  growth  in  InGaAs  or  InAlAs 
layers  thinner  than  tc,MB-  Hence,  one  is  always  safe  in  staying  below  the  Matthews- 
Blakeslee  limit.  In  fact,  for  device  or  circuit  applications  in  which  a  single  misfit 
dislocation  would  be  catastrophic,  tc,MB  should  probably  be  used  as  an  upper  limit. 
For  applications  in  which  a  relatively  low  density  of  misfit  dislocations  is  tolerable 
{e.g.  HFETs),  it  may  be  possible  to  exceed  the  severe  constraints  imposed  by  the 
Matthews-Blakeslee  limit.  Klem  has  shown  that  the  misfit  dislocation  density  is  a 
function  of  the  density  of  threading  dislocations  in  the  substrate  for  InGaAs  growth 
on  GaAs.^^’^  Hence,  one  should  begin  by  choosing  substrates  with  the  lowest  possible 
density  of  threading  dislocations.^  As  mentioned  earlier,  there  are  several  potential 
sources  of  heterogeneous  nucleation  including  oval  defects  and  particulates.  At  this 
point  in  time,  it  is  not  clear  which  sources  are  dominant.  We  will  not  attempt  to 
review  the  literature  on  substrate  preparation  and  elimination  of  oval  defects,  but  one 
should  take  steps  to  minimize  probable  nucleation  sources.  Finally,  surface  roughness 
and  3-D  growth  should  be  avoided. 

Considerable  progress  has  been  made  in  the  modelling  of  crystalline  growth  modes 
by  molecular  dynamics  simulations.^®' To  our  knowledge,  however,  there  is  not  yet 
a  simple  way  to  predict  the  onset  of  3-D  growth  or  surface  roughening  for  a  given 
material  and  set  of  growth  conditions.  Experimental  evidence  suggests  that  3-D 
growth  does  not  occur  until  thicknesses  well  beyond  tc,MB  fo^  both  InGaAs 
and  InAlAs.**®  Growth  temperature  plays  a  role,  with  3-D  growth  occurring  sooner  at 
higher  temperatures.®®’^**’*^®  Hence,  growth  temperatures  as  low  as  possible*  should  be 
used.  If  surface  roughness  and  3-D  growth  are  avoided,  and  the  density  of  nucleation 
sources  is  low,  epilayers  with  high  crystalline  quality  and  thermal  stability  should 
result. 

I  As  shown  by  Fitzgerald  ei  alP  and  Watson  ei  the  use  of  patterned  substrates  can  reduce 

the  dislocation  density  by  limiting  the  number  of  threading  dislocations  available. 

^If  the  growth  temperature  is  too  low,  material  quality  will  degrade,  independent  of  strain. 
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6.5  Summary 

In  summary,  we  have  shown  that  the  contradictory  results  from  studies  of  thermal 
stability  of  strained  layers  can  be  explained  by  the  heterogeneous  nucleation  of  dis¬ 
location  loops.  Our  high  quaUty  InGaAs  and  InAlAs  layers  were  thermally  stable 
because  of  a  relatively  low  density  of  nucleation  sources  in  our  MBE  process.  We  pro 
pose  that  if  a  layer  undergoes  surface  roughening  or  3-D  growth,  additional  nucleation 
sources  are  made  available,  allowing  substrate  relaxation  during  growth  or  subsequent 
annealing.  Based  upon  these  findings,  we  have  suggested  procedures  for  the  growth  of 
high-qi'ality,  thermally-stable  strained  layers  beyond  the  Matthews-Blakeslee  cntical 

layer  thickness . 
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Chapter  7 


Summary 


7.1  Conclusions 

We  used  high-resolution  x-ray  diffraction  (HRXRD)  extensively  in  this  work  and 
found  it  to  be  a  powerful  and  versatile  characterization  tool.  It  is  routinely  used  to 
measure  the  composition  of  single  epitaocial  layers.  We  showed  that  if  the  thickness  of 
a  calibration  layer  is  chosen  properly,  the  InAs  mole  fraction  in  InGaAs  or  InAlAs  on 
InP  can  be  determined  to  within  0.01  with  a  single  measurement.  If  HRXRD  scans 
are  combined  with  dynamical  simulations,  one  can  determine  the  composition  in  com¬ 
plex  multi-layer  heterostructures.  In  addition,  HRXRD  gives  information  about  the 
crystalline  quality  of  epitaxial  layers.  A  good  figure  of  merit  is  the  ratio  of  the  exper¬ 
imental  to  the  theoretical  epilayer  peak  width,  with  values  close  to  unity  representing 
high  crystalline  quality.  The  presence  of  interference  fringes  also  indicates  a  coher¬ 
ent,  high-quality  layer.  These  criteria  can  detect  degradation  in  crystalline  quality 
earlier  than  direct  measurements  of  lattice  relaxation.  HRXRD  measurements  are 
non-destructive  and  relatively  simple  and  fast  in  most  cases. 

Variable  azimuthal- angle  ellipsometry  (VAAE)  is  a  new  characterization  technique 
which  measures  the  optical  anisotropy  of  epitaxial  layers.  Nearly  lattice-matched 
layers  of  InGaAs  and  InAlAs  always  exhibit  a  small  anisotropy  which  is  caused  by 
surface  chemistry  effects.  Layers  with  a  large  mismatch  often  exhibit  large  optical 
anisotropies.  We  have  attributed  these  large  values  to  strain-induced  roughness.  This 
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roughness  is  correlated  with  a  degradation  in  crystalline  quality  measured  by  HRXRD. 
In  some  cases,  the  morphology  consists  of  islands  elongated  in  the  [110]  direction. 
VAAE  is  simple,  non-destructive,  and  extremely  fast.  It  emerges  as  a  potentially 
powerful  technique  for  probing  the  limits  of  strained  layer  epitaxy. 

Using  the  HRXRD  FWHM  ratio  and  interference  fringes,  we  found  that  the  crys- 
taUine  quality  of  InGaAs  and  InAlAs  epitaxial  layers  consistently  remains  unper¬ 
turbed  to  thicknesses  up  to  3-9  x  <c,MB  for  lattice  mismatch  up  to  ±  1%.  For 
InAlAs,  layers  in  tension  can  withstand  larger  mismatch  than  layers  in  compression 
(at  constant  thickness).  For  InGaAs,  we  observe  no  clear  differences  between  com¬ 
pression  and  tension.  Layers  which  are  beyond  tc.MJB  but  have  high  crystalline  quahty 
may  have  a  low  density  of  misfit  dislocations  at  the  epilayer /substrate  interface.  Both 
TEM  and  VAAE  show  that  such  layers  have  good  surface  morphology,  suggesting  2-D 
growth.  In  summary,  the  Matthews-Blakeslee  critical  layer  thickness  is  not  a  sharp 
limit  beyond  which  epilayer  quality  abruptly  degrades. 

Asymmetric  HRXRD  measurements  show  that  partially  relaxed  epitaxial  layers 
of  InGaAs  and  InAlAs  on  (001)  InP  exhibit  orthorhombic  distortion.  The  maximum 
strain  relief  is  in  the  [110]  direction  for  either  material  in  tension  or  compression. 
Symmetric  HRXRD  measurements  reveal  anisotropic  peak  broadening  for  certain 
samples,  with  a  maximum  peak  width  in  the  [110]  direction.  All  observations  are 
consistent  with  a  model  in  which  strain  is  reheved  by  60®  misfit  dislocations,  with 
a  majority  parallel  to  the  [lIO]  direction.  TEM  observations  on  selected  samples 
confirm  this  model.  For  other  samples,  the  situation  may  be  more  complex,  with 
surface  roughness  and  islands  elongated  in  the  [lIO]  direction.  The  asymmetry  is 
attributed  to  the  fact  that  dislocations  along  the  [110]  and  [lIO]  directions  have 
different  dislocation  cores  in  zinc-blende  structures.  Orthorhombic  distortion  must 
be  taken  into  account  when  measuring  the  composition  of  partially  relaxed  epitaxial 
layers  by  HRXRD. 

High-quahty  layers  of  InGaAs  and  InAlAs  with  thicknesses  up  to  9  x  <c,mb  were 
annealed  at  temperatures  from  600  to  850®C.  HRXRD  measurements  before  and  after 
the  anneals  showed  no  degradation  of  crystalhne  quahty.  In  addition,  MODFET 
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structures  with  strained  InGaAs  channels  were  annealed  and  characterized  by  Hall 
mobility  measurements.  The  300K  mobilities  were  not  affected  by  anneals.  Some 
degradation  in  the  77K  mobility  occurred,  but  it  is  apparently  the  result  of  diffusion 
of  Si  dopant  atoms  into  the  channel  rather  than  lattice  relaxation.  We  attribute  the 
thermal  stability  of  our  strained  layers  to:  1)  insufficient  energy  to  homogeneously 
nucleate  half-loop  dislocations,  2)  insufficient  number  of  threading  dislocations  in 
the  substrate,  and  3)  insufficient  number  of  sites  (possibly  impurity  clusters  or  oval 
defects)  for  heterogenerous  nucleation  of  dislocations. 

Based  upon  TEM,  VAAE,  and  HRXRD  measurements,  a  roughening  transition 
occurs  when  the  layer  thickness  and  mismatch  become  too  large.  This  transition  may 
be  equivalent  to  the  so-called  2-D  to  3-D  transition  other  authors  have  observed  using 
RHEED.  We  cannot,  however,  confirm  that  true  3-D  growth  results  in  all  cases.  One 
possibility  is  that  surface  roughness  is  sometimes  an  intermediate  stage  between  2-D 
growth  with  a  smooth  morphology  and  true  3-D  or  island  growth.  We  propose  that 
surface  roughening  provides  defect  sites  for  the  heterogeneous  nucleation  of  dislocation 
loops.  Layers  which  reach  the  roughening  transition  can  then  relax  substantially.  This 
model  is  also  consistent  with  the  thermal  stability  of  high-quality  strained  layers:  such 
layers  did  not  reach  the  roughening  transition  during  growth. 

The  fact  that  thermally  stable  layers  with  high  crystalline  quality  can  be  grown 
to  thicknesses  greater  than  suggests  that  such  layers  may  be  suitable  for  de¬ 

vice  applications.  We  grew  and  fabricated  InAlAs/InGaAs/InP  HFETs  in  which  the 
InAlAs  pseudo-insulator  was  AlAs-rich,  with  =  2.0.  The  device  results  were 

excellent  and  showed  no  signs  of  degradation  due  to  lattice  mismatch.  Others  have 
fabricated  similar  HFETs  with  strained  InGaAs  channels  exceeding  <c,AfB  by  a  factor 
of  2-4  without  degradation.  These  results  suggest  that  the  thickness-mismatch  limit 
for  InP-based  HFETs  may  be  approximated  by  the  boundary  between  regions  II  and 
III  in  fig.  6.2.  The  limits  may  be  different  for  other  devices,  but  it  is  clear  that  the 
Matthews-Blakeslee  critical  layer  thickness  should  not  be  taken  as  an  absolute  upper 
limit  in  all  cases. 
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7.2  Suggestions  for  Future  Work 


Strained  layer  epitaxy  is  likely  to  remain  an  important  field  for  many  years  to  come. 
Its  challenges  range  from  fundamental  problems  of  materials  science  to  practical  ap¬ 
plications  for  electronic  and  optical  devices.  Based  upon  the  research  described  here, 
we  have  identified  three  promising  areas  which  could  benefit  from  additional  work. 

In  chapter  6,  we  proposed  a  model  relating  lattice  relaxation,  material  quality,  and 
growth  modes.  A  few  key  experiments  could  confirm  or  refute  this  model.  First,  the 
thermal  stability  of  samples  from  region  III  of  fig.  6.2  (surface  roughening  and  crys- 
taJhne  degradation  without  substantial  relaxation  during  growth)  should  be  tested. 
Our  model  predicts  that  relaxation  will  occur  during  annealing  because  of  the  de¬ 
fects  associated  with  surface  roughness.  Second,  additional  sources  of  heterogeneous 
nucleation  sites  should  be  introduced  into  samples  from  region  II  (high-crystalline 
quality  with  t  >  <c,mb)-  Damage  could  be  created  by  ion  implantation  after  growth 
and  before  anneahng.  Sample  surfaces  could  be  scratched  before  growth  or  between 
growth  and  anneahng.  Other  possibiHties  include  the  use  of  substrates  with  much 
higher  densities  of  threading  dislocations  (p  >  10®/cm^),  and  the  growth  of  epitaxial 
layers  with  extremely  high  doping  densities  (n  or  p  >  10^®/cm®).  We  would  expect 
some  if  not  aU  of  these  procedures  to  result  in  unstable  epitaxial  layers.  HRXRD 
could  be  used  as  the  primary  characterization  tool  to  assess  thermal  stabihty.  TEM 
could  be  apphed  to  selected  samples  to  determine  relaxation  mechanisms. 

Additional  studies  of  the  origin  of  optical  anisotropy  could  help  to  make  VAAE  a 
standard  characterization  tool.  Correlation  of  VAAE  results  with  plan- view  and  cross- 
sectional  TEM  as  weU  as  optical  microscopy  should  lead  to  a  better  understanding  of 
the  growth  morphologies  which  cause  optical  anisotropy.  Thus  far,  we  have  only  used 
VAAE  to  characterize  samples  after  growth.  VAAE  could  also  be  applied  as  an  in  situ 
characterization  technique.  This  would  allow  measurements  on  each  layer  in  complex 
heterostructures  and  might  provide  enough  warning  so  that  strained  layer  growth 
can  be  halted  before  severe  degradation  occurs.  VAAE  could  be  applied  in  either  an 
MBE  or  MOCVD  system.  It  might  be  especiaUy  useful  in  MOCVD  since  RHEED  is 
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not  available.  Some  groups  have  already  installed  ellipsometers  on  MOCVD  systems; 
VAAE  measurements  could  be  made  with  little  additional  effort. 

We  have  investigated  asymmetries  in  x-ray  and  optical  measurements.  Oth¬ 
ers  have  observed  electrical  anisotropies  in  strained  III-V  layers  as  well.  It  is  not 
yet  clear  whether  the  electrical  anisotropies  result  from  asymmetric  arrays  of  mis¬ 
fit  dislocations,  surface  roughness,  3-D  growth,  or  some  combination.  Comparison 
to  TEM  measurements  should  settle  this  question.  Mobilities  must  be  measured 
with  lithographically-defined  Hall  patterns  to  separate  geometry  effects  from  true 
anisotropies.  Such  measurements  may  be  very  sensitive  to  the  onset  of  dislocations. 
They  could  be  combined  with  annealing  to  investigate  thermal  stability.  Actual  de¬ 
vice  measurements  are,  of  course,  more  time-consuming,  but  can  provide  a  wealth  of 
additional  information.  For  example,  our  findings  of  thermal  stability  based  upon  the 
mobility  of  strained-channel  MODFET  heterostructures  should  be  confirmed  by  fab¬ 
ricating  and  testing  MODFETs  on  samples  with  and  without  post-growth  annealing. 
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Appendix  A 

Molecular  Beam  Epitaxy  of 
InGaAs  and  InAlAs 


The  InGaAs  and  InAlAs  layers  characterized  in  this  study  were  grown  by  molec¬ 
ular  beam  epitaxy  (MBE).  In  general,  we  used  established  growth  and  calibration 
procedures.  The  procedures  are  described  in  this  appendix. 

InP  substrates  from  Sumitomo  Electric  were  used  for  aU  growths.  They  were 
5  cm  in  diameter,  350  /im  thick,  and  nominally  within  0.5°  of  the  (001)  orientation. 
The  InP  was  doped  with  Fe  to  make  it  semi-insulating,  with  a  resistivity  of  at  least 
10®  ohm-cm.  The  wafers  were  pre-etched  and  shipped  in  nitrogen.  Wafers  were 
scribed  and  mounted  on  molybdenum  blocks  with  indium.  We  did  no  wafer  cleaning 
or  etching  except  for  blowing  the  samples  with  nitrogren  after  mounting. 

Epitaxial  layers  were  grown  in  a  Riber  2300  solid  source  MBE  system.  A  schematic 
of  the  evaporation  chamber  is  shown  in  fig.  A-1.  The  base  pressure  of  this  chamber 
is  typically  1-2  xl0“^°  Torr  (after  cooling  with  liquid  nitrogen).  The  system  also  in¬ 
cludes  preparation  and  load-lock  chambers.  The  cell  temperatures,  shutter  positions, 
and  substrate  temperature  are  computer-controlled.^^® 

Samples  were  baked  in  the  preparation  chamber  at  200-250° C  for  45-60  min¬ 
utes  prior  to  growth.  They  were  then  transferred  into  the  evaporation  chamber  and 
ramped  up  to  growth  temperature  under  an  As  overpressure.  Samples  were  rotated 
during  layer  growth.  The  ratio  of  the  group  V  (As)  to  group  III  (In,  Al,  and  Ga) 
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Figure  A-1:  Schematic  of  Riber  2300  Molecular  Beam  Epitaxy  system  used  in  this 
study  (courtesy  of  J.C.  Vlcek). 

beam-equivalent-pressures  was  between  15:1  and  25:1  during  growth,  with  background 
pressures  in  the  lO"®  Torr  range. 

Sample  temperature  was  monitored  by  a  thermocouple  and  an  optical  pyrometer 
during  growth.  The  two  measurements  often  differed  by  as  much  as  50“C.  Surface 
reconstruction  patterns  were  monitored  by  reflection  high-energy  electron  diffraction 
(RHEED)  during  ramp-up.  Layers  were  always  grown  at  temperatures  above  the 
2x4  InP  reconstruction  and  typically  5-25°C  below  the  4x2  InP  reconstruction. 
We  estimate  that  the  growth  temperature  was  always  between  460  and  blO^C. 

In  MBE,  it  is  rare  for  a  source  molecule  to  strike  a  gas  molecule  because  of  the 
extremely  low  pressure.  Hence,  the  process  can  be  modelled  as  line-of-sight  evapora¬ 
tion  into  a  perfect  vacuum.  The  flux  of  a  species  D  evaporated  into  a  perfect  vacuum 
is  given  by:^^^ 


V2Tr  MvkTD 

where  Pr,  is  the  equilibrium  vapor  pressure  of  D,  Md  is  the  mass  of  species  D,  Td  is 
the  absolute  temperature,  and  k  is  the  Boltzmann  constant.  The  vapor  pressure  is 
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an  exponential  function  of  temperature: 


Pd  = 


(A.2) 


where  Ea.d  is  the  activation  energy  and  P o,p  is  a  constant  independent  of  tempera¬ 
ture.  Combining  (A.l)  and  (A.2)  for  Ga  evaporation: 


Foa  = 


(A.3) 


y/^-KMaJiTca 

We  grew  a  series  of  In*Gai_a,As  layers  with  a  constant  In  cell  temperature  and 
varying  Ga  cell  temperature.  For  this  case,  assuming  that  the  mole  fractions  of  GaAs 
and  InAs  are  proportional  to  their  fluxes  (no  re-evaporation),  we  have: 


1  -  jg  ^  F^  _  B  ^-E^.gJkToa  (A. 4) 

X  Fin  s/Tca 

To  a  good  approximation,  we  can  neglect  the  \/f  term  in 
1  ~  ^  S  (A.5) 

X 

where  A  is  a  constant. 

A  plot  of  the  logarithm  of  (l-x)/x  versus  inverse  absolute  temperature  should 
approximately  yield  a  straight  line  with  a  slope  of  {-EaI^)-  In  Ag-  A-2,  we  show  the 
data  for  13  samples  of  In*Gai_*As  grown  from  the  same  charges  of  Ga  and  In.  The 
composition  of  each  InGaAs  layer  was  measured  by  HRXRD,  as  described  in  chapter 
2.  A  least-squares  fit  gives  an  activation  energy  of  2.19  eV  or  25,340  K  for  Ga.*  We 
assume  that  the  activation  energy  (slope)  remains  constant  after  cell  recharging.  The 
intercept  of  the  line  (or  the  lattice-matching  temperature)  typicaUy  changes  when  the 
cells  are  recharged. 

Given  the  Ga  cell  temperature,  the  growth  rate  can  also  be  estimated.  We  again 

tThe  activation  energy  of  a  cell  can  also  be  determined  from  flux  measurements.  Such  measure¬ 
ments  gave  a  value  of  2.12  eV  for  the  Ga  cell. 


where  5  is  a  constant. 
(A. 4)  and  obtain: 
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1000/T^(K-') 


Figure  A-2:  Calibration  of  the  composition  of  InxGai_*As:  plot  of  ratio  of  GaAs 
mole  fraction  to  InAs  mole  fraction  versus  inverse  absolute  Ga  cell  temperature  (In 
cell  temperature  was  constant)  for  13  samples.  The  slope,  2.19  eV,  is  the  activation 
energy  for  the  Ga  cell  (see  eq.  A. 5).  Composition  was  measured  by  HRXRD. 


assume  that  all  the  Ga  and  In  will  be  absorbed  on  the  substrate  surface.  The  growth 
rate,  G,  wiU  then  be  proportional  to  the  sum  of  the  Ga  and  In  fluxes: 


G  =  C{Fin  +  FGa)  (A. 6) 

where  G  is  a  constant.  We  assume  that  the  growth  rate,  Go,  and  composition,  Xo, 
are  known  for  the  Ga  temperature,  To’  (To  is  often  the  temperature  which  yields 
lattice- matched  InGaAs  on  InP,  but  it  can  be  any  temperature.)  Combining  (A. 4) 
and  (A. 6): 


Go  =  CFg.{To) 


«0 


1  -  ®0 


+  1 


(A.7) 


We  now  wish  to  determine  the  growth  rate,  Gi,  at  the  Ga  temperature  T\.  We 
obtain  the  composition,  xi,  from  (A. 5)  or  fig.  A-2.  Using  (A. 3),  we  can  approximate 
the  ratio  of  Ga  fluxes  at  Ti  and  To 
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(A.8) 


Fgc{Ti)  ^  J-E,^/k)(llTi-l/To) 

FoaiTo) 

We  can  now  express  the  growth  rate  as: 

G.  =  C1F,„  +  Fo,(r,)l  =  CFo.(r.)  f-^  +  ll  (A.9) 

Li  —  Xi 

Combining  (A.8)  and  (A.9): 


1  —  ®i 


(A.IO) 


Given  the  activation  energy  for  the  Ga  cell  and  a  cahbration  sample  for  which  the 
composition  and  thickness  have  been  measured,  the  composition  and  growth  rate  can 
be  estimated  for  any  Ga  cell  temperature  from  (A. 5)  and  (A.IO).  We  measured  the 
thickness  of  calibration  samples  by  one  of  three  techniques:  profilometry  after  selective 
etching,  HRXRD  fringe  spacing  (see  section  2.1),  and  infrared  reflectance.^"  Using 
the  above  procedures,  we  were  usually  able  to  grow  InxGai_aAs  with  x  within  0.01 
of  the  desired  value.  Thickness  was  typically  within  5-10%  of  the  targeted  value. 

Procedures  analogous  to  those  above  were  used  for  In^Alj—yAs  calibration,  but  the 
error  in  y  was  sometimes  as  large  as  0.02.  The  reason  for  the  larger  errors  appears  to 
be  fluctuations  in  the  activation  energy  of  the  A1  ceU.  Our  HFET  heterostructures 
are  an  ideal  way  to  measure  the  A1  activation  energy  because  they  typically  have 
well-defined  HRXRD  peaks  from  both  the  buffer  and  insulator  layers  (see  fig.  5.2). 
Values  of  Ea,ai  determined  from  HFET  structures  ranged  from  2.4  to  3.4  eV. 

A  key  parameter  in  our  HFETs  is  the  electron  concentration  in  the  channel.  We 
used  Hall/ van  der  Pauw  measurements  of  electron  concentration  to  calibrate  Si  dop¬ 
ing  of  InGaAs.  A  plot  of  electron  concentration  versus  inverse  Si  cell  temperature 
yielded  a  Si  cell  activation  energy  of  4.8  eV  for  n  =  3  x  10^®  to  3  x  10^®/cm®.  In  this 
doping  range,  however,  we  were  only  able  to  consistently  obtain  the  desired  electron 
concentration  within  a  factor  of  two.  We  suspect  the  variations  are  due  to  the  ampho¬ 
teric  nature  of  Si  in  III-V’s.  The  incorporation  of  Si  is  known  to  be  sensitive  to  the 
growth  temperature  which,  as  discussed  above,  is  not  constant  for  all  growths. 
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Appendix  B 

Material  Parameters  of  InP, 
GaAs,  InAs,  and  AlAs* 


Property 

Parameter  (units) 

InP 

GaAs 

InAs 

AlAs 

Lattice  constant 

ao(A) 

5.8688 

5.6532 

6.0584 

5.6623 

Density 

P  (g/cm^) 

4.787 

5.307 

5.667 

3.81 

Thermal  expansion 

atfc(xlO-V'’C') 

4.56 

6.63 

5.16 

Elastic  constant 

Cii(xlO^  ^dyn/cm*) 

10.22 

11.88 

8.33 

15.1 

Elastic  constant 

Ci2(  X  lO^^dyn/cm*) 

5.76 

5.37 

4.53 

5.7 

Elastic  constant 

C44(  X 10^  ^  dyn  jcm^) 

4.6 

5.94 

3.96 

5.89 

Young’s  modulus 

Y  (xlO^^dyn/cm^) 

9.49 

12.40 

7.93 

16.50 

Shear  modulus 

G  ( X 10^^  dyn/cm*) 

3.02 

4.15 

2.59 

5.1 

Poisson’s  ratio 

u 

0.36 

0.311 

0.352 

0.274 

Band  gap 

E,  (=V) 

1.35 

1.42 

0.36 

2.15 

Electron  mass 

m*/mo 

0.08 

0.067 

0.023 

0.14 

Light-hole  mass 

0.089 

0.074 

0.027 

0.15 

Heavy-hole  mass 

ia.hh/”^o 

0.85 

0.62 

0.6 

0.76 

Electron  mobility 

fie{cm^  lY -sec) 

4500 

8500 

28000 

Hole  mobility 

fth(cm*/V-8ec) 

130 

400 

450 

Static  dielectric  constant 

e./eo 

12.4 

13.1 

14.6 

8.2 

Optical  dielectric  constant 

^oo/  ^0 

9.55 

11.1 

12.3 

10.1 

Refractive  index  (1.96  eV) 

n 

3.54 

3.86 

3.96 

3.11 

Extinction  coef.  (1.96  eV) 

k 

0.30 

0.20 

0.61 

0.0015 

Absorption  coef.  (1.96  eV) 

a{fcm) 

61000 

39000 

120000 

290 

*  Parameters  from  references®®’ 
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Appendix  C 

Matthews-Blakeslee  Critical  Layer 
Thickness 


In  this  appendix,  we  briefly  outline  the  Matthews-Blakeslee  theory  of  critical  layer 
thickness  and  modify  it  to  account  for  anisotropic  elastic  properties.  Matthews  and 
Blakeslee  estimated  the  critical  layer  thickness,  tc,  by  equating  the  force  exerted  by 
misfit  strain  with  the  tension  in  a  dislocation  line.^  The  result  was  a  transcendental 
equation:* 


6(1  —  1/  cos^  a)  A  tc  \ 
^  87r/(l  +  i/)cosA  (  V 


(C.l) 


where  6  is  the  magnitude  of  the  Burger’s  vector  of  the  dislocation,  /  is  the  lattice 
mismatch,  a  is  the  angle  between  the  dislocation  line  and  its  Burger’s  vector,  A  is 
the  angle  between  the  slip  direction  and  that  direction  in  the  film  plane  which  is 
orthogonal  to  the  line  of  intersection  of  the  slip  plane  and  the  interface,  and  u  is 
Poisson’s  ratio  which  is  related  to  the  elastic  stiffness  coefficients  by: 


1/  = 


'12 


(C.2) 


Cn  +  C'i2 

For  III-V  semiconductors,  experimental  work  has  shown  that  misfit  dislocations 


‘Matthews  and  Blakeslee  originally  considered  a  special  type  of  strained-layer  superlattice  for 
which  the  right-hand-side  of  (C.l)  has  an  additional  factor  of  4  in  the  numerator. 
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are  primarily  of  the  60®  type  on  (001)  substrates. For  this  case,  cosA  =  cosa  =  1/2, 
b  =  ae/-\/2,  and  (C.l)  becomes: 


t  = 

4\/2rr/ (1  + 1/) 


tte 


(C.3) 


where  a*  is  the  natural  (relaxed)  lattice  constant  of  the  epilayer. 

To  simphfy  calculations,  Matthews“®-“^  assumed  that  the  epitaxial  layer  was 
elastically  isotropic  and  that  the  shear  modulus  of  the  substrate  was  equal  to  the  shear 
modulus  of  the  epilayer.  Hirth^®®  and  Fitzgerald^®®  modified  the  Matthews  theory  to 
consider  elastic  anisotropy,  a  general  feature  of  crystalline  materials.  Fitzgerald  also 
considered  unequal  shear  moduh,  and  obtained: 


tc  = 


O'e  (1  ~  *^/4)  GeGs 


v'2»/(l-.)y(G^  +  G.)>“— -"M  (C.4) 

where  Y  is  Young’s  modulus  of  the  epilayer,  and  Gs  and  Ge  are  the  shear  moduli 
of  the  substrate  and  epilayer,  respectively.  In  the  case  of  biaxial  stress  applied  to  an 
anisotropic  material,  Y  and  G  are  given 


r  =  Cxi  +  ^  (C.5) 

and 

G  =  Cu  —  -(2C44  +  C12  —  Cii)  (C.6) 

We  plot  the  critical  thickness  from  eq.  (C.3)  (isotropic)  and  eqs.  (C.4)-(C.6) 
(anisotropic)  in  figs.  C.l  (InGaAs)  and  C.2  (InAlAs).  Throughout  this  thesis,  we  use 
the  anisotropic  critical  thickness  defined  by  (C.4)  and  refer  to  it  as  tc,MB- 
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Figure  C-1:  Matthews-Blakeslee  critical  layer  thickness  for  InGaAs  layers  on  InP. 


Figure  C-2;  Matthews-Blakeslee  critical  layer  thickness  for  InAlAs  layers  on  InP. 
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